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DELOCALIZATION FOR A CLASS OF RANDOM BLOCK BAND MATRICES 


ZHIGANG BAG AND LASZLO ERDOS 


Abstract. We consider N x N Hermitian random matrices H consisting of blocks of size M > 

The matrix elements are i.i.d. within the blocks, close to a Gaussian in the four moment matching sense, 
but their distribution varies from block to block to form a block-band structure, with an essential band 
width M. We show that the entries of the Green’s function G{z) = {H — z)~^ satisfy the local semicircle 
law with spectral parameter z = E + ir) down to the real axis for any 77 3> N~^, using a combination 
of the supersymmetry method inspired by [19] and the Green’s function comparison strategy. Previous 
estimates were valid only for rj ;§> M~^. The new estimate also Implies that the eigenvectors in the 
middle of the spectrum are fully delocalized. 


1. Introduction 

The Hamiltonian of quantum systems on a graph T is a self-adjoint matrix H = {hab)a,ber, H = H*. 
The matrix elements hab represent the quantum transition rates from vertex a to b. Disordered quantum 
systems have random matrix elements. We assume they are centered, l^hab = 0, and independent subject 
to the basic symmetry constraint hab = hba- The variance := represents the strength of the 

transition from a to b and we use a scaling where the norm ||R|| is typically order 1. The simplest case is 
the mean field model, where hab are identically distributed; this is the standard Wigner matrix ensemble 
[29]. The other prominent example is the Anderson model [2] or random Schrddinger operator, H = A-l-H, 
where the kinetic energy A is the (deterministic) graph Laplacian and the potential V = {Vx)xer is an 
on-site multiplication operator with random multipliers. If T is a discrete d-dimensional torus, then only 
few matrix elements hab are nonzero and they connect nearest neighbor points in the torus, dist(a, b) < 1. 
This is in sharp contrast to the mean field character of the Wigner matrices. 

Random band matrices naturally interpolate between the mean held Wigner matrices and the short 
range Anderson model. They are characterized by a parameter M, called the band width, such that the 
matrix elements hab for dist(a, b) > M are zero or negligible. If M is comparable with the diameter N of 
the system then we are in the mean held regime, while M ~ 1 corresponds to the short range model. 

The Anderson model exhibits a metal-insulator phase transition: at high disorder the system is in 
the localized (insulator) regime, while at small disorder it is in the delocalized (metallic) regime, at 
least in d > 3 dimensions and away from the spectral edges. The localized regime is characterized 
by exponentially decaying eigenfunctions and off diagonal decay of the Green’s function, while in the 
complementary regime the eigenfunctions are supported in the whole physical space. In terms of the 
localization length i, the characteristic length scale of the decay, the localized regime corresponds to 
£ N, while in the delocalized regime £ ^ N. Starting from the basic papers [1, 15], the localized regime 
is well understood, but the delocalized regime is still an open mathematical problem for the d-dimensional 
torus. 
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Since the eigenvectors of the mean field Wigner matrices are always delocalized [13, 14], while the 
short range models are localized, by varying the parameter M in the random band matrix, one expects 
a (de)localization phase transition. Indeed, for d = 1 it is conjectured (and supported by non rigorous 
supersymmetric calculations [16]) that the system is delocalized for broad bands, M and localized 

for M <C The optimal power 1/2 has not yet been achieved from either sides. Localization has 

been shown for M <g; in [22], while delocalization in a certain sense was proven for M ^ in 

[11]. Interestingly, for a special Gaussian model even the sine kernel behavior of the 2-point correlation 
function of the characteristic polynomials could be proven down to the optimal band width M ^ 
see [18, 20]. Note that the sine kernel is consistent with the delocalization but does not imply it. We 
remark that our discussion concerns the bulk of the spectrum; the transition at the spectral edge is much 
better understood. In [24] it was shown that the edge spectrum follows the Tracy Widom distribution, 
characteristic to mean field model, for M ^ but it yields a different distribution for narrow bands, 

M < 

Delocalization is closely related to estimates on the diagonal elements of the resolvent G{z) = {H—z)~^ 
at spectral parameters with small imaginary part rj = \mz. Indeed, if Ga^E + iij) is bounded for all i 
and all Li € M, then each £^-normalized eigenvector u of is delocalized on scale rj~^ in a sense that 
maxi \ui\'^ < r], i.e. u is supported on at least r]~^ sites. In particular, if Ga can be controlled down to 
the scale rj ^ 1/N, then the system is in the complete delocalized regime. Moreover, boundedness of Gu 
also implies that the local semicircle law holds for the same regime of ry. 

For band matrices with band width M, or even under the more general condition < M~^, the 
boundedness of Gu was shown down to scale ry ^ M~^ in [14] (see also [12]). If M ^ it is 

expected that Gu remains bounded even down to jy 3> N~^ which is the typical eigenvalue spacing, the 
smallest relevant scale in the model. However, with the standard approach [14, 12] via the self-consistent 
equations for the Green’s function does not seem to work for ry < 1 /M ; the fluctuation is hard to control. 
The more subtle approach using the self-consistent matrix equation in [11] could prove delocalization 
and the off-diagonal Green’s function profile that are consistent with the conventional quantum diffusion 
picture, but it was valid only for relatively large ry, far from M~^. Moment methods, even with a delicate 
renormalization scheme [23] could not break the barrier rj ^ M~^ either. 

In this paper we attack the problem differently; with supersymmetric (SUSY) techniques. Our main 
result is that Gu{z) is bounded, and the local semicircle law holds for any ry i.e. down to the 

optimal scale, if the band width is not too small, M ^ but under two technical assumptions. 

First, we consider a generalization of Wegner’s n-orbital model [21, 28], namely, we assume that the band 
matrix has a block structure, i.e. it consists oi M x M blocks and the matrix elements within each block 
have the same distribution. This assumption is essential to reduce the number of integration variables in 
the supersymmetric representation, since, roughly speaking, each M x M block will be represented by 
a single supermatrix with 16 supersymmetric variables. Second, we assume that the distribution of the 
matrix elements matches a Gaussian up to four moments in the spirit of [27]. Supersymmetry heavily 
uses Gaussian integrations, in fact all mathematically rigorous works on random band matrices with 
supersymmetric method assume that the matrix elements are Gaussian, see [4, 5, 6, 18, 19, 20, 25, 26]. 
The Green’s function comparison method [14] allows one to compare Green’s functions of two matrix 
ensembles provided that the distributions match up to four moments and provided that Gu are bounded. 
This was an important motivation to reach the optimal scale ry ^ N~^. 

In the next subsections we introduce the model precisely and state our main results. While SUSY 
approach is ubiquitous in physics, see e.g. the basic monograph by Efetov [7], its application in rigorous 
proofs is notoriously difficult. Initiated by T. Spencer (see [25] for a summary) and starting with the 
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paper [4] by Disertori, Pinsker and Spencer, only a handful of mathematical papers have succeeded in 
exploiting this powerful tool. Our supersymmetric analysis was inspired by [19], but our observable. Gab, 
requires a partly different formalism, in particular we use the singular version of the superbosonization 
formula [3]. Moreover, our analysis is considerably more involved since we consider relatively narrow 
bands. In Section 1.3, we explain our novelties compared with [19]. 

1.1. Matrix model. Let 


Hn = (hab) 

he a,n N X N random Hermitian matrix, in which the entries are independent (up to symmetry), centered, 
complex variables. In this paper, we are concerned with Hn possessing a block band structure. To define 
this structure explicitly, we set the additional parameters M = M(N) and W = W (N) satisfying 

W = N/M. 

For simplicity, we assume that both M and W are integers. Let S = (Sjk) he a W xW symmetric matrix, 
which will be chosen as a weighted Laplacian of a connected graph on W vertices. Now, we decompose 
Hn into W x W blocks of size M x M, and relabel 

hjk,a^ ■— hab, j,h — I 5 ■ ■ • ; hV, cx, (3 — 1, . . . , ilL, 

where {j,k) = ([a/M], [ 6 /M]) is the index of the block containing hab, and 

{a, 13) = (a - (j - 1)M, b - {k - 1)M^ 

describes the location of the entry in the block. Moreover, we assume 

^hji^a/3hj'k',a'^' — jlz> 6j'6 {5jk 3~Bjk)- (^■^) 

That means, the variance profile of the random matrix \/MHn is given by 

S = (Sj/c) := I + S, (1-2) 

in which each entry represents the common variance of the entries in the corresponding block of \/MHn- 
Moreover, if hjk,ai 3 ^s are Gaussian, (1.1) also implies that for each off-diagonal entry hjk^ap, its real part 
and imaginary part are i.i.d. N{0,Sjk/2M) variables. 

1.2. Assumptions and main results. In the sequel, for some matrix A = (oy) and some index sets 
I and J, we introduce the notation to denote the submatrix obtained by deleting the i-th row and 

j-th column of A for all* G I and j G J. We will adopt the abbreviation 

ylbb) ^ ;= yl(«IW). ( 1 . 3 ) 

In addition, we use UAjjmax := rnaxi^- ja^ j to denote the max norm of A. Throughout the paper, we need 
some assumptions on S. 

Assumption 1.1 (On S). Let Q = (V,£l) he a connected simple graph with V = {1,... ,W}. Assume 
that S is a W xW symmetric matrix satisfying the following four conditions. 

(i) S is a weighted Laplacian on Q, i.e. for i j, we have Sij > 0 if {i,j} G S and Sij = 0 if 
{i,j} ^ £, and for the diagonal entries, we have 

^ii — ^ ^ ^ij, 


\/i = l,...,W. 
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(ii) S defined in (1.2) is strictly diagonally dominant, i.e., there exists some constant cq > 0 such 
that 

l + 2Sii>co, M i = l,...,W. 

(iii) For the discrete Green’s functions, we assume that there exist some positive constants C and 7 
such that 

max ||(5«)-i|Uax<C'iy^. 

(iv) There exists a spanning tree Go = (V,fo) C G, on which the weights are bounded below, i.e. for 
some constant c > 0 , we have 

5ij>c, if {i,j}G£o 

Remark 1.2. From Assumption 1.1 (ii), we easily see that 

S > col (1.4) 

for the same positive constant cq. In addition, the lower bound c in (iv) can be weakened to N~^ for some 
sufficiently small constant £ > 0. But for simplicity, we will not try to optimize this bound in this paper. 

Example 1.1. Let A he the standard discrete Laplacian on the d-dimensional torus fl with 

periodic boundary condition, where tr = W/d. Here by standard we mean the weights on the edges of the 
box are all 1. Now let S = aA for some positive constant a < l/4d. It is then easy to check Assumption 
1.1 (i), (ii) and (iv) are satisfied. In addition, if d = 1, it is well known that we can choose 7 = 1 m 
Assumption 1.1 (iii). For d > 3, one can choose 7 = 0 . For d = 2, one can choose 7 = £ for arbitrarily 
small constant e. For instance, one can refer to [8] for more details. 

For simplicity, we also introduce the notation 

cr^- := r := (cTyOAT.Ar = ImI'm. = (1.5) 

where 1m is the M-dimensional vector whose components are all 1 and S is the variance matrix in (1.2). 
It is elementary that 

Spec(r) = Spec(.S)U{0}c [0,1]. (1.6) 

Our assumption on M depends on the constant 7 in Assumption 1.1 (iii). 

Assumption 1.3 (On M). We assume that there exists a (small) positive constant £1 such that 

M (1.7) 

Remark 1.4. A direct consequence of (1.7) and N = MW is 

4+2T+ei 

M > N . (1.8) 

Especially, when 7 = 1 , one has M Actually, through a more involved analysis, (1.7) (or (1.8)) 

can he further improved. At least, for 7 < 1, we expect that M is enough. However, we will not 

pursue this direction here. 

Besides Assumption 1.1 on the variance profile of H, we need to impose some additional assumption on 

the distribution of its entries. To this end, we temporarily employ the notation iJ® = (/i®j) to represent a 

random block band matrix with Gaussian entries, satisfying (1.1), Assumption 1.1 and Assumption 1.3. 
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Assumption 1.5 (On distribution). We assume that for each a,b G {1,..., N}, the moments of the 
entry hab match those of h^j^ up to the ^th order, i.e. 

E{Rehab)'"i\mhabY ='E{RehlbfO^KbY^ Vfc,£eN, s.t. k + £ < 4. (1.9) 

In addition, we assume the distribution of hab possesses a subexponential tail, namely, there exist positive 
constants ci and C 2 such that for any 7 > 0 , 

p{\hab\ > rnnhabY)^) < e^e-^ ( 1 . 10 ) 

holds uniformly for all a,b = 1,..., N. 

The four moment condition (1.9) in the context of random matrices first appeared in Tao and Vu’s 
work [27]. 

To state our results, we will need the following notion on the comparison of two random sequences, 
which was introduced in [9] and [12]. 

Definition 1.6 (Stochastic domination). For some possibly N-dependent parameter set UjV; o-nd two 
families of random variables X = (Xjv(u) : IV e N,it G Uw) and Y = {Yn(u) : N G N,u € Uat), we say 
that X is stochastically dominated by Y, if for all e > 0 and D > 0 we have 

sup p(x 7 v(m) > N^Yn{u)) < N-^ (1.11) 

«6Ujv ^ 

for all sufficiently large N > iVo(e,D). In this case we write 

X ^ Y. 


For example, by (1.1) and Assumption 1.5, we have 

\hab\-<^, Va, 6 =l,...,A. (1.12) 

Note that S is doubly stochastic. It is known that the empirical eigenvalue distribution of Hn converges 
to the semicircle law, whose density function is given by 

Qsc[x) := ^ - a ;2 • l(jx| < 2 ). 

We denote the Green’s function oi Hn by 

G{z) = Gn{z) := {Hm — z)~^, z = E + ip G C'*' := {w G C : Imw > 0} 


and its (a, 6 ) matrix element is Gabiz). Throughout the paper, we will always use E and p to denote 
the real and imaginary part of z without further mention. In addition, for simplicity, we suppress the 
subscript N from the notation of the matrices here and there. The Stieltjes transform of Qscix) is 

-z + - 4 


msc{z) = j 




£1 

-2 X- z 


da; = 


where we chose the branch of the square root with positive imaginary part for z G C“'". Note that msc(z) 
is a solution to the following self-consistent equation 


msc{z) = 


Yz)' 


(1.13) 


The semicircle law also holds in a local sense, see Theorem 2.3 in [12]. For simplicity, we cite this 
result with a slight modification adjusted to our assumption. 
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Proposition 1.7 (Erdos, Knowles, Yau, Yin, [12]). Let H be a random block band matrix satisfying 
Assumptions 1.1, 1.3 and 1.5. Then 


max \Gab{z) - Sabmsc{z)\ ^ 

a,b 


1 

y/Mr]' 


if 


for any fixed small positive constants k and e. 


E € [—2 “t- R, 2 — k\ 


and M < ?? < 10 


(1.14) 


Remark 1.8. We remark that Theorem 2.3 in [12] was established under a more general assumption 
^k^jk ~ 1 '^jk — G/M. Especially, the block structure on the variance profile is not needed. In 
addition, Theorem 2.3 in [12] also eovers the edges of the spectrum, which will not be discussed in this 
paper. We also refer to [I]] for a previous result, see Theorem 2.1 therein. 

Our aim in this paper is to extend the local semicircle law to the regime rj ^ N~^ and replace M with 
N in (1.14). More specifically, we will work in the following set, defined for arbitrarily small constant 
K > 0 and any sufficiently small positive constant £2 '.= e 2 (£i), 

D(iV, K, e2 ):=\z = EEip€£-.\E\<y/2- n, < p < (1.15) 

Throughout the paper, we will assume that £2 is much smaller than £ 1 , see (1.7) for the latter. Specifically, 
there exists some large enough constant C such that £2 < SijG. 


Theorem 1.9 (Local semicircle law). Suppose that H is a random block band matrix satisfying Assump¬ 
tions 1.1, 1.3 and 1.5. Let k be an arbitrarily small positive constant and £2 be any sufficiently small 
positive eonstant. Then 

max \Gab{z) - Sabmsc{z) \ -< 1— (1-16) 

a,b \/Nr] 

holds uniformly on D(7V, ^,£ 2 ). 

Remark 1.10. The restriction \E\ < •\/2 — k in (1.15) is technical. We believe the result can be extended 
to the whole bulk regime of the spectrum, i.e., \E\ <2 — k, see Section 12 for further comment. The upper 
bound of r] in (1.15) is also technical. However, for rj > one can control the Green’s function 

by ( 1 . 14 ) directly. 


Let Ai,..., Atv be the eigenvalues oi H^. We denote by := {un ,..., UiN) the normalized eigenvector 
of H]\[ corresponding to Ai. From Theorem 1.9, we can also get the following delocalization property for 
the eigenvectors. 


Theorem 1.11 (Complete delocalization). Let H be a random block band matrix satisfying Assumptions 
1.1, 1.3 and 1.5. We have 

max ||ui||oo^lV“ 2 . (1-17) 

i-.\\i\<^/2-K 

Remark 1.12. We remark that delocalization in a certain weak sense was proven in [11] for an even 
more general class of random band matrices if M However, Theorem 1.11 asserts delocalization 

for all eigenvectors in a very strong sense (supremum norm), while Proposition 7.1 of [11] stated that 
most eigenvectors are delocalized in a sense that their substantial support cannot be too small. 


1.3. Outline of the proof strategy and novelties. In this section, we briefly outline the strategy for 
the proof of Theorem 1.9. 

The first step, which is the main task of the whole proof, is to establish the following Theorem 1.14, 
namely, a prior estimate of the Green’s function in the Gaussian case. For technical reason, we need the 
following slight modification of Assumption 1.3, to state the result. 
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Assumption 1.13 (On M). Let ei be the small positive constant in Assumption 1.3. We assume 

N{\0gN)-^° >M> ll/4+27+ei^ 

In the regime M > N{\ogN)~^^, we see that (1.16) anyway follows from (1.14) directly. 

Theorem 1.14. Assume that H is a Gaussian block band matrix, satisfying Assumptions 1.1 and 1.13. 
Let n be any fixed positive integer. Let k he an arbitrarily small positive constant and £2 be any sufficiently 
small positive constant. There is Nq = NQ(n), such that for all N > Nq and all z € D(iV, ^,£ 2 ), we have 

E\Gab{z)\^^ <N^°{Sab + j^), ya,b=l,...,N (1.19) 

for some positive constant Co independent of n and z. 

Remark 1.15. Much more delicate analysis can show that the prefactor N^° can be improved to some 
n-dependent constant Cn. We refer to Section 12 for further comment on this issue. 

Using the definition of stochastic domination in Definition 1.6, a simple Markov inequality shows that 
(1.19) implies 

\Gab{z)\ ^ Sab - 7 x 7 =) V a, 6 = 1 ,..., A. ( 1 - 20 ) 

VNr) 

The proof of Theorem 1.14 is the main task of our paper. We will use the supersymmetry method. We 
partially rely on the arguments from Shcherbina’s work [19] concerning universality of the local 2-point 
function and we develop new techniques to treat our observable, the high moment of the entries of G{z), 
under a more general setting. We will comment on the novelties later in this subsection. 

The second step is to generalize Theorem 1.14 from the Gaussian case to more general distribution 
satisfying Assumption 1.5, via a Green’s function comparison strategy initiated in [14], see Lemma 2.1 
below. 

The last step is to use Lemma 2.1 and its Gorollary 2.2 to prove our main theorems. Using (1.20) below 
to bound the error term in the self-consistent equation for the Green’s function, we can prove Theorem 
1.9 by a continuity argument in z, with the aid of the initial estimate for large ry provided in Proposition 
1.7. Theorem 1.11 will then easily follow from Theorem 1.9. 

The second and the last steps are carried out in Section 2. The main body of this paper. Sections 3-11 
is devoted to the proof of Theorem 1.14. 

One of the main novelty of this work is to combine the supersymmetry method and the Green’s 
function comparison strategy to go beyond the Gaussian ensemble, which was so far the only random 
band matrix ensemble amenable to the supersymmetry method, as mentioned at the beginning. The 
comparison strategy requires an apriori control on the individual matrix elements of the Green’s function 
with high probability (see (1.20)), this is one of our main motivations behind Theorem 1.14. 

Although we consider a different observable than [19], many technical aspects of the supersymmetric 
analysis overlaps with [19]. For the convenience of the reader, we now briefly introduce the strategy of 
[19], and highlight the main novelties of our work. 

In [19], the author considers the 2-point correlation function of the trace of the resolvent of the 
Gaussian block band matrix H, with the variance profile S' = 1 -I- oA, under the assumption M ^ N 
(note that we use M instead of W in [19] for the size of the blocks). The 2-point correlation function 
can be expressed in terms of a superintegral of a superfunction F({Si})^]^) with a collection of 4 x 4 
supermatrices Si := Z*Zi. Here for each i, Zi = 'I' 2 y, 4 > 2 ,i) is an M x 4 matrix and Z* is 

its conjugate transpose, where and 4 ' 2 ,i are Grassmann M-vectors whilst <i)i ^ and $ 2 ,i are complex 
M-vectors. Then, by using the superbosonization formula in the nonsingular case (M > 4) from [17], 
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one can transform the superintegral of to a superintegral of where each Si is a 

supermatrix akin to Si, but only consists of 16 independent variables (either complex or Grassmann). 
We will call the integral representation of the observable after using the superbosonization formula as the 
final integral representation. Schematically it has the form 



for some functions g(-), fc(-) and fg(-), where we used the abbreviation S := {Si}J^i and Sc and Sg 
represents the collection of all complex variables and Grassmann variables in S, respectively. Here, 
g{Sc) and fc(>5o) are some complex functions and ^g{Sg,Sc) will be mostly regarded as a function of the 
Grassmann variables with complex variables as its parameters. The number of variables (either complex 
or Grassmann) in the final integral representation then turns out to be of order W, which is much smaller 
than the original order N. In fact, in [19] it is assumed that W = 0(1) although the author also mentions 
the possibility to deal with the case W ~ for some small positive e, see the remark below Theorem 1 
therein. 

Performing a saddle point analysis for the complex measure exp{Mfc(iSc)}, one can restrict the integral 
in a small vicinity of some saddle point, say, Sc = Sco- It turns out that fc(5co) = 0 and fc(5c) decays 
quadratically away from Sco- Consequently, by plugging in the saddle point Sco, one can estimate g(5c) 
by g(Sco) directly. However, for exp{Mfc(5o)} and exp{fg(5g, 5c)}, one shall expand them around the 
saddle point. Roughly speaking, in some vicinity of Sco, one will find that the expansions read 

^Mfc(Sc) = exp{—u'Au + ec(u)}, gfal^g.^o) = exp{—p'IHIx}p(p, r, u), (1-22) 

where u is a complex vector of dimension 0{W), which is essentially a vectorization of '/M{Sc — Sco)', 
ec(u) = o(l) is some error term; p and r are two Grassmann vectors of dimension 0(W); A is a complex 
matrix with positive-definite Hermitian part and H is a complex matrix; p(p, t,u) is the expansion of 
exp{fg(5g, 5c) — fg(5g,5co)}, which possesses the form 

0{W) 

PiP,r,u)= ^ (1.23) 

£=0 

where pe{p,T,u) is a polynomial of the components of p and t with degree 2£, regarding u as fixed 
parameters. Now, keeping the leading order term of p(p, t,u), and discarding the remainder terms, one 
can get the final estimate of the integral by taking the Gaussian integral over u, p and t. This completes 
the summary of [19]. 

Similarly to [19], we also use the superbosonization formula to reduce the number of variables and 
perform the saddle point analysis on the resulting integral. However, owing to the following three main 
aspects, our analysis is significantly different from [19]. 


• (Different observable) Our objective is to compute high moments of the single entry of the Green’s 
function. By using Wick’s formula (see Proposition 3.1), we express E|Ggfc|^” in terms of a superintegral 
of some superfunction of the form 


F 


(f 


VI/ ilr* . 

^ ^ a,j 


a=1.2; ) := {4>l,q,p4>l,p,a4>2,p,a4’‘2,q,p) 


for some p,q £ {1,..., IT} and a, f G {1,..., M}, where is the o-th coordinate of and the 

others are defined analogously. Unlike the case in [19], F is not a function of {Si}^i only. Hence, using 
the superbosonization formula to change Si to Si directly is not feasible in our case. In order to handle 
the factor ((^i,q,/ 3 (()i,p,a<A 2 .p.a<(> 2 .(j,/ 3 )", the main idea is to split off certain rank-one supermatrices from Sp 
and Sq such that this factor can be expressed in terms of the entries of these rank-one supermatrices. 
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Then we use the superbosonization formula not only in the nonsingular case from [17] but also in the 
singular case from [3] to change and reduce the variables, resulting the final integral representation of 
E|C?jfep”. Though this final integral representation, very schematically, is still of the form (1.21), due to 
the decomposition of the supermatrices Sp and Sg, it is considerably more complicated than its counter¬ 
part in [19]. Especially, the function g(5c) differs from its counterpart in [19], and its estimate at the 
saddle point follows from a different argument. 

• (Small band width) In [19], the author considers the case that the band width M is comparable with TV, 
i.e. the number of blocks W is finite. Though the derivation of the 2-point correlation function is highly 
nontrivial even with such a large band width, our objective, the local semicircle law and delocalization 
of the eigenvectors, however, can be proved for the case M ^ in a similar manner as for the Wigner 
matrix {M = N), see [12, 14]. In our work, we will work with much smaller band width to go beyond 
the results in [12, 14], see Assumption 1.3. Several main difficulties stemming from a narrow band width 
can be heuristically explained as follows. 

At first, let us focus on the integral over the small vicinity of the saddle point, in which the exponential 
functions in the integrand in (1.21) approximately look like (1.22). 

We regard the first term in (1.22) as a complex Gaussian measure, of dimension 0{W). When IT ~ 1, 
one can discard the error term ec(u) directly and perform the Gaussian integral over u, due to the fact 
f duexp{—u'Re(A)u}|ec(u)| = o(l). However, such an estimate is not allowed when W ^ (say), 
because the normalization of the measure exp{—u'Re(A)u} might be exponentially larger than that of 
exp{—u'Au}. In order to handle this issue, we shall do a second deformation of the contours of the 
complex variables in the vicinity of the saddle, following the steepest descent paths exactly, whereby we 
can transform the complex Gaussian measure to a real one, thus the error term of the integral can be 
controlled. 

Now, we turn to the second term in (1.22). When IT ^ 1, there are only hnitely many Grassmann 
variables. Hence, the complex coefficient of each term in the polynomial p(p, t,u), which is of order 
for some £ £ N (see (1.23)), actually controls the magnitude of the integral of this term against 
the Gaussian measure exp{—p'HIr}. Consequently, in case of IT ^ 1, it suffices to keep the leading order 
term (according to one may discard the others trivially, and compute the Gaussian integral over 

p and T explicitly. However, when IT ~ N'^ (say), in light of the Wick’s formula (3.2) and the fact 
that the coefficients are of order the order of the integral of each term of p(/9, x,u) against the 

Gaussian measure reads detH*^'!-') for some index sets I and J and some £ £ N. Due to the fact 

IT - N^, detHClJ) is typically exponential in IT. Hence, it is much more complicated to determine and 
compare the orders of the integrals of all terms. In our discussion, we perform a unified estimate 

for the integrals of all the terms, rather than simply compare them by 

In addition, the analysis for the integral away from the vicinity of the saddle point in our work is 
also quite different from [19]. Actually, the integral over the complement of the vicinity can be trivially 
ignored in [19], since each factor in the integrand of (1.21) is of order 1, thus gaining any o(l) factor for 
the integrand outside the vicinity is enough for the estimate. However, in our case, either exp{Mfc(iSc)} 
or f dSg exp{fg(Sg, Sc)} is essentially exponential in IT. This fact forces us to provide an apriori bound 
for f dSg exp{fg(Sg,Sc)} in the full domain of Sc rather than in the vicinity of the saddle point only. In 
addition, an analysis of the tail behavior of the measure exp{Mfc(5c)} needs also to be performed. 
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• (General variance profile S) In [19], the authors considered the special case S = aA with a < l/4d. We 
generalize the discussion to more general weighted Laplacians S satisfying Assumption 1.1, which, as a 
special case, includes the standard Laplacian A for any fixed dimension d. 

1.4. Notation and organization. Throughout the paper, we will need some notation. At first, we 
conventionally use U{r) to denote the unitary group of degree r, as well, U{1, 1) represents the t/(l, 1) 
group. Furthermore, we denote 

U{r) = U{r)/U{1Y, U{1, 1) = t/(l, 1)/C/(1)2. 

Recalling the real part E oi z, we will frequently need the following two parameters 

iE + V4:-E^ iE - V4- A2 

a+ =-^-, a_ =-^-. 

Correspondingly, we define the following four matrices 

= diag(a+, a_), Hip = diag(a_, a+), Z1+= diag(a+, a+), £)_ = diag(a_, a_). (1-24) 

We remark here D± does not mean “£>+ or For simplicity, we introduce the following notation for 

some domains used throughout the paper. 

I:=[0,1], L:=[0,27r), S : unit circle, M+:=[0, oo), R_ := —R+, F := a+R+. 

For some i x i Hermitian matrix A, we use Ai(A) < ... < Xi{A) to represent its ordered eigenvalues. 
For some possibly iV-dependent parameter set Uat, and two families of complex functions {uNiu) : 
N G N,u G and {6jv(m) ■ N G N,u G Uat}, if there exists a positive constant C > 1 such that 
C~^\b]s[(u)\ < |a 7 v(M)| < C\b]s[(u)\ holds uniformly in N and u, we write 

ajv(u) ^ bjv(u). 

Conventionally, we use {e^ : i = !,...,£} to denote the standard basis of in which the dimension £ 
has been suppressed for simplicity. For some real quantities a and b, we use a Ab and a V & to represent 
min{a, b} and max{a, b}, respectively. 

Throughout the paper, c, c', ci, C 2 , C, C, Ci, C 2 represent some generic positive constants that are 
possibly n-dependent and may differ from line to line. In contrast, we use Co to denote some generic 
positive constant independent of n. 

The paper will be organized in the following way. In Section 2, we prove Theorem 1.9 and Theorem 1.11, 
with Theorem 1.14. The proof of Theorem 1.14 will be done in Section 3-Section 11. More specifically, 
in Section 3, we use the supersymmetric formalism to represent EjGyp” in terms of a superintegral, 
in which the integrand can be factorized into several functions; Section 4 is devoted to a preliminary 
analysis on these functions; Section 5-Section 10 are responsible for different steps of the saddle point 
analysis, whose organization will be further clarified at the end of Section 5; Section 11 is devoted to the 
final proof of Theorem 1.14, by summing up the discussions in 3-Section 10. Finally, in Section 12, we 
make some further comments on possible improvements. 

2. Proofs of Theorem 1.9 and Theorem 1.11 

At first, (1.19) can be generalized to the generally distributed matrix with the four moment matching 
condition via the Green’s function comparison strategy. 
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Lemma 2.1. Assume that H is a random block band matrix, satisfying Assumptions 1.1, 1.5 and 1.13. 
Let K be an arbitrarily small positive constant and £2 be any sufficiently small positive constant. There 
is Nq = iVo(n), such that for all N > Nq and all z € 'D(N, k, £ 2 ), we have 

+ Va,&=l,...,iV (2.1) 

for some positive constant Co uniform in n and z. 

By the definition of stochastic domination in Definition 1.6, we can get the following corollary from 
Lemma 2.1 immediately. 

Corollary 2.2. Under the assumptions of Lemma 2.1, we have 

\Gab{,z)\ -< 6ab H—7 ^f=^ ^ a,b = 1 ,..., (2-2) 

y/Np 

uniformly on T){N, k,£ 2 ). 

In the sequel, at first, we prove Lemma 2.1 from Theorem 1.14 via the Green’s function comparison 
strategy. Then we prove Theorem 1.9, using Lemma 2.1. Finally, we will show that Theorem 1.11 follows 
from Theorem 1.9 simply. 

2.1. Green’s function comparison: Proof of Lemma 2.1. To show (2.1), we use Lindeberg’s re¬ 
placement strategy to compare the Green’s functions of the Gaussian case and the general case. That 
means, we will replace the entries of if® by those of H one by one, and compare the Green’s functions 
step by step. Ghoose and fix a bijective ordering map 

: {(*, j) : 1 < * < j < {l,-■ •,‘?(^)}, c(iV) := (2.3) 

Then we use Hk to represent the NxN random Hermitian matrix whose (i, j)-th entry is hij if w(i, j) < k, 
and is hL otherwise. Especially, we have Hq = and = H. Correspondingly, we define the 

Green’s functions by 

Gk[z)-.= (Hk-zy\ k=l,...,g{N). 

Fix /c and denote 

vj~^{k) = {a,b). (2.4) 

Then, we write 

Hk-i = Hi + V,6, := (1 - {hl,e,el + hle.el ), 

Hk = HI +\Nab, ^ab ■= - {hab^a^b + , 

where Hi is obtained via replacing hab and hba by 0 in Hk (or replacing and hf^ by 0 in Hk-i). In 
addition, we denote 

GUz) = {Hl-zr\ 

Set £3 = £ 3 ( 7 , £ 1 ) to be a sufficiently small positive constant, satisfying (say) 

< — gi 
- 100 ' 5 -b 27 -h £1 ’ 


(2.5) 
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where 7 is from Assumption 1.1 (iii) and ei is from (1.7). For simplicity, we introduce the following 
parameters for £ = 1,..., <?(iV) and i,j = 1,..., iV, 




( 2 . 6 ) 


where C is a positive constant. Here we used the notation 5|j = 1 if two index sets I and J are the same 
and (5|j = 0 otherwise. It is easy to see that for rj < we have 


07*1 < 200, V£=l,...,^(iV), ^,j = l,...,iV, 


(2.7) 


by using ( 1 . 8 ). 

Now, we compare Gk-i{z) and Gk{z). We will prove the following lemma. 


Lemma 2.3. Suppose that the assumptions in Lemma 2.1 hold. Additionally, we assume that for some 
sufficiently small positive constant 63 satisfying (2.5), 

\{Gi)i,{z)\<N^f |(G°),,(z)HiV-^ V£=l,...,c(iV), Vi,j = l,...,iV (2.8) 

uniformly for z G 'D{N, k, 62 ). Let n gN be any given integer. Then, if 

E|(Gfe_i)i,(z)|2"<0fe_i.,,(<5,, + ^^), V*,j = l,...,iV, (2.9) 

we also have 

E|(Gfc)y(z)|2-<0fe,,,(j,, + ^^), y^,J = l,...,N (2.10) 

for any k = 1,..., <;{N). 


Proof of Lemma 2.3. Fix k and omit the argument z from now on. At first, under the conditions (2.8) 
and (2.9), we show that 

E|(G°).,r <30o(% + ^^), V*,j = l,...,7V. (2.11) 

To see this, we use the expansion with (2.4) 

(G°)., = [Gk-ih + (Gfc_iW,,G°).„ 


which implies that for a sufficiently large constant D > 0 

_/V2e3 ,2n 


in ^ / \ 

n{Glf,?- < E |(Gfc_i).,| + < 30o(5., + j^), 


1 

{N^y 


where the first step follows from (1.12), (2.8), Definition 1.6 and the trivial bound 7 ^ for the Green’s 
functions, and the second step follows from (2.9), (2.7) and the fact /\/M <C 1/y/TVrj for z G 

D(N, K,e2)- 

Now, recall (2.4) again and expand Gfe_i(z) and Gk{z) around G^{z), namely 


Gfc_i = G° + ^(-1)'(G°V,6)'G° + (-l)'"+i(G°V„,)™+iGfe_i, 

1^1 

m 

Gk = Gl + ^(-l)^(G°W,fc)^G° + {-ir+yGl\N,tr+^Gk. (2.12) 

1=1 


We always choose m to be sufficiently large, depending on £3 but independent of N. Then, we can write 

171 

(Gfc_i)ij = {G^)ij + + Rm+l,ij, 

e=i 
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where 


m 

i=l 


(2.13) 


Ra.. := (-1)'((G°V,,)'G°)^^, S,.., := (-l)'((G°W,b)^G°)^^, £ = 1,... ,m, 

R^+i,„ := (-!)'"+! ((G°V,b)™+iGfe_i)^^, S^+i.i, := (-!)'"+'((G°W„fc)'"+iGfc)^^. (2.14) 

At first, by taking m sufficiently large, from (2.8) and (1.12), we have the trivial bound 

|R„+i.,,|, ^ M-^7V(™+2)e3 « (2.15) 

For Re^ij and Se^ij., we split the discussion into off-diagonal case and diagonal case. In the case off ^ j, we 
keep the first and the last factors of the terms in the expansions of {{G^\/abYGk)ij iiGk'^abYG^)ij, 
namely, {G^)ij' and [GY)i'j for some i',j' = a or b, and bound the factors in between by using (1.12) and 
(2.8), resulting the bound 

^ |(G0).,.(G0).,,|, e = l,...,m. (2.16) 

i' ^j'—a,b 

For i = j, we only keep the first factor of the terms in the expansions of {{G°\/abYGk)ii and {{G°\NabYGk)ii, 
and bound the others by using (1.12) and (2.8), resulting the bound 

\R,MASeM<M-^N^^^(^\{GlU + \{Glu), i = (2.17) 

Observe that, in case i ^ j, if {i,j} ^ {a,b}, at least one of {G^)ij' and (G°)i'j is an off-diagonal entry 
of G^ for i',j' = a or b. 

Now we compare the 2n-th moment of \{Gk-i)ij \ and \iGk)ij\- At first, we write 


E|(Gd)„f”=E((Gd),,r((Gd)„r, d = k-l,k (2.18) 

By substituting the expansion (2.13) into (2.18), we can write 

E|(Gd)yf" = A(^,j) + Rd(f,J), d=k-l,k, (2.19) 


where A(f,j) is the sum of the terms which depend only on and the first four moments of hab, and 
'R.d{i,j) is the sum of all the other terms. We claim that Rd(f, j) satisfies the bound 




d = fc — 1, fc. 


( 2 . 20 ) 


(A? 7 )"-§ ' {NrtY 

for some positive constant G. Now, we verify (2.20). According to (2.11) and the fact that the sequence 
Ri^ij,..., Rm,ij, Rm+i,ij^ as well as Si^^,..., Sm,ij, Sm-i-i.ii, decreases by a factor /^YM in magnitude, 
it is not difficult to check the leading order terms of Rfc_i(f, j) are of the form 


E((G?)..)"((Gg).,) 

and those of Rfe(f, j) are of the form 

E((G°)y 

with some p,qi,q[ € N such that 


2n-p-ELl(«+9^) 

J. X 


-i 0 \ \P 


Glh 


2n-p-J2Yiile+qe) 


n cge cl'e 


£=1 


'^iiqe + q'i) =5, 0 < p < 2n -'^{qe + q'^). 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


e=i 


e=i 
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Every other term has at least 6 factors of hab or or their conjugates, thus their sizes are typically 
controlled by M~^{Nr])~'^, i.e. they are subleading. Hence, it suffices to bound (2.21) and (2.22). In the 
sequel, we only estimate (2.21) in details, (2.22) can be handled in the same manner. 

Now, the five factors of hab or hba within the R^yj’s in (2.21) are independent of the rest and estimated 
by For the remaining factors from G°., we use (2.11) to bound 2n of them and use (2.8) to 

bound the rest. In the case that i ^ j and {i,j} ^ {a, f*}, by the discussion above, we must have an 
off-diagonal entry of in the product {G^)ij' {G^)i'j for any choice of i',/ = a or b. Then, in the bound 
for Ri^ij in (2.16), for each {G^)ij'{G^)i'j, we keep the off-diagonal entry and bound the other by 
from assumption (2.8). Hence, by using (2.16) and (2.23), we see that for some v,jV S {i,j, a, 6} with 
ir ^ jr, r = 1 ,..., ^(< 7 ^ + the following bound holds 


(2.21) < f 


N^3 . 




E 


(^0).^. |2n-ELl(«+9^) 


-liQe+Qe) 

n 


,|)^3( 


00 

^/M) {NijY' 


(2.24) 


where the last step follows from (2.11) and Holder’s inequality. In case oi i Y j but {i,j} = {a,b}, we 
keep an entry in the product (GY)ij' i bound the other by . We remark here in this case 
the entry being kept can be either diagonal or off-diagonal. Consequently, for some ir,jr G {b j) a, 6}, r = 
1 ,... ,^{qi + q'f)^ we have the bound 


/ \ 5 

^ ( tm ) ^ 




n 


') ^ 


0 n 


T-rOr \ j — 




(2.25) 


by using (2.11) and Holder’s inequality again. Hence, we have shown (2.20) in the case oi i Y j- For 
i = j, it is analogous to show 


Ktm) 


(2.26) 


by using (2.11), (2.17) and Holder’s inequality. Hence, we verified (2.20). Consequently, by Assumption 
1.5, (2.19) and (2.20) we have 


E|(Gfe-i),,f" 


E|(Gfc).,f" 



{Nr]Y-i 


1 

{Nrj)'^ 




which together with the assumption (2.9) for E|(Gfc_i)ijp" and the definition of Qe^ij’s in (2.6), we can 
get 

Hence, we completed the proof of Lemma 2.3. □ 


To show (2.1), we also need the following lemma. 

Lemma 2.4. Suppose that the assumptions in Lemma 2.1 hold. Fix the indices a,b ^ {1,...A^}. Let 
be a matrix obtained from H with its (a, b)-th entry replaced by 0. Then, if for some tjq > 1/N there 
exists 

\Gu{z)\^l, |(G°)i,( 2 ;)H 1 for r]>rio, 'di = l,...,N, 

then we also have 

\G,j{z)\<^, \{G°)ij{z)\for ^<77 < 770 , W i,j = 1,..., N. 


(2.27) 
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Proof of Lemma 2.4- The proof is almost the same as the discussion on pages 2311-2312 in [10]. For the 
convenience of the reader, we sketch it below. At first, according to the discussion below (4.28) in [10], 
for any i,j = 1,..., iV, we have 

\Gij{E + ir;)] < Cmax^ \ mGu{E + 12^17). 

^ k>0 

Now, we set 


ki := max{fc : 2 ^7] < rjo}, ^2 := max{fc : 2^77 < 1 }. 

According to our assumption, both ki and ^2 are of the order log N. Now, we have 

All ^2 00 

\mGu{E + i2'^T]) = Y lmG«(A + 12^) + Y '<^Gu{E + 12^) + Y '<^Gu{E + i2'^'n) 

A :>0 A ;—0 A:—Aii A:—A12 + I 

^ — Y] 7^^'‘^Gu{E + irjo) + (^2 ~ ^i) + 1 ^ — 

^ ko ^ ^ 

where in the second step, we used the fact that the function y y\mGu{E + iy) is monotonically 
increasing, the condition (2.27) and the fact rj < tjq. Hence, we conclude the proof of Lemma 2.4. □ 

Now, with Theorem 1.14, Lemma 2.3 and Lemma 2.4, we can prove Lemma 2.1. 


Proof for Lemma 2.1. The proof relies on the following bootstrap argument, namely, we show that once 

l(G,), 7 l^l, V£ = l,...,<iV), Vf,i = l,...,iV (2.28) 

holds for r] > rjo with tjq G it also holds for rj > r]QN~^^ for any £3 satisfying (2.5). 

Assuming (2.28) holds for y > rjo, we see that 


maxl(G°)yl = maxl(G^)y + {{Gi)\N abG'l)i 


^ max l(G^)jj 


^=maxl(G^)yl •maxl(G^)yl ^ 1 + maxKG?)^^]. 

VM hj vM 


Consequently, for rj > rjo, we also have 

1(G°),,1^1, V£ = l,...,^(iV), Vz,j = l,...,iV. 


(2.29) 


Therefore, (2.27) holds. Then, by Lemma 2.4, we see that (2.8) holds for 77 > tjqN . Furthermore, by 
Lemma 2.3 and Theorem 1.14 for Gq, i.e. the Gaussian case, one can get that for any given n, 

E] (G^),, [2- < (<5,, + < 200 (<5y + , for M-^N^->y> , 

V£ = l,...,<7(iV), V7,j = l,...,iV. (2.30) 

Note that since (2.30) holds for any given n, we get (2.28) for >y> 

Now we start from 770 = By Proposition 1.7 we see that (2.28) holds for all 77 > 770 . Then 

we can use the bootstrap argument above finitely many times to show (2.28) holds for all 77 > 
Consequently, we have (2.8) for all 77 > Then, Lemma 2.1 follows from Lemma 2.3 and Theorem 

1.14 immediately. □ 
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2.2. Proof of Theorem 1.9. Without loss of generality we can assume that M < N{\ogN)~^^, other¬ 
wise, Proposition 1.7 implies (1.16) immediately. Now, recalling the notation defined in (1.3), we denote 
the Green’s function of as 

G«(z) := (i7« -z)-\ 

with a little abuse of notation. We only need to consider the diagonal entries Gu below, since the bound 
for the off-diagonal entires of G{z) is implied by (2.1) directly. Set 

/S.i =/S.i{z) := — - \- z + a^iGag. (2-31) 

^ii 

a 

We introduce the notation 

Ad = Ad{z) := max \Gii{z) - msc{z)\. 

i 

We have the following lemma. 

Lemma 2.5. Suppose that H satisfies Assumptions 1.1, 1.5 and 1.13. We have 

|A( 2 )H^, V*=l,...,iV, (2.32) 

uniformly for z € D(Af, k, £ 2 ). 

The proof of Lemma 2.5 will be postponed. Using Lemma 2.5, we see that, with high probability, 
(2.31) is a small perturbation of the self-consistent equation of nisc, i-C. (1-13), considering = 1- 

To control we use a continuity argument from [12]. 

We remind here that in the sequel, the parameter set of the stochastic dominance is always D(A^, r, £ 2 ), 
without further mention. We need to show that 

A, ^ (2.33) 

and first we claim that it suffices to show that 

l{Ad<N-^)Ad^^. (2.34) 

Indeed, if (2.34) were proven, we see that with high probability either A^ > N~~r or A^ ^ < 

N~^ for z G D(Af, k,£ 2 ). That means, there is a gap in the possible range of Ad. Now, choosing £ in 
(1.14) to be sufficiently small, we are able to get for 77 = 

Ad^N-^, y Eg[-2 + k,2-k\, Vi=l,...,W (2.35) 

By the fact that Ad is continuous in z, we see that with high probability, A^ can only stay in one side of 
the range, namely, (2.33) holds. The rigorous details of this argument involve considering a hne discrete 
grid of the 2 -parameter and using that G{z) is Lipschitz continuous (albeit with a large Lipschitz constant 
1 / 77 ). The details are found in Section 5.3 of [12]. 

Hence, what remains is to verify (2.34). The proof of (2.34) is almost the same as that for Lemma 3.5 
in [14]. For the convenience of the reader, we sketch it below without reproducing the details. We set 

1 ^ 

m = m{z) :=—'^G^i{z), Ui = 0 ^( 2 ) := - m, i = l,...,N. 

' i=l 

We also denote u:=(ui,...,UAr). By the assumption Ad < N~~i', we have 

u, = 0 (Af-T). 


(2.36) 
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Now we rewrite (2.31) as 

0 = Gii + 


=--Gu 


+ Gli. 


^ + Ea """ ' Z + m{z) 

By using (2.32), Lemma 5.1 in [14], and the assumption < N~~r ^ we can show that 


12 ,; = — 


+ 0(||u||^) +0(max|Z\i|). 


(2.37) 


(2.38) 


{z + m(z))2 

One can refer to the derivation of (5.14) in [14] for more details. Averaging over i for (2.37) and (2.38) 
leads to 

1 


i{z) 


z + m{z) 


= -n. 


and 


1 ^ 

12 := = 0 (||ull^) + 0 (maxlAl) 


(2.39) 


(2.40) 


i=l 


Plugging (2.36) and (2.32) into (2.40) yields 

\n\^N-^. (2.41) 

Using (2.41), the fact \rh(z) — msc{z)\ < Ad < N~~i', and Lemma 5.2 in [14], to (2.39), we have 

l7Ti(2;) - msc{z)\ < |12| = 0(1 Ju]]^) + O(maxlziil), (2.42) 

where in the first step we have used the fact that z G D(iV, ^,£ 2 ) thus away from the edges of the 
semicircle law. Now, we combine (2.37), (2.38) and (2.39), resulting 


" (&SSe ^ ^ ^ ^ mfx 


i = l,...,iV. (2.43) 


We just take the above identity as the definition of w^. Analogously, we set w := (wi,..., w^r)'. Then 

(2.40) and (2.43) imply 

||w||oo = O(llull^) +0(maxlz\il) +O(llulloo • |(z + 7 ti(z))“^ - (z + msc(z))~^l) 

< O(llull^) +0(maxlZ\il) +O(llulloo • |m(z) - msc(z)l) 

< 0(||u 11^)+0(maxlA|) (2.44) 

where the second step follows from the fact \z + msc{z)\ > 1 inD(A^,K,e 2 ) (see (5.1) in [14] for instance), 

(2.41) and (2.42), and in the last step we used (2.42) again. 

Now, using the fact m1^{z) = (msc{z) + z)“^ (see (1.13)), we rewrite (2.43) in terms of the matrix T 
introduced in (1.5) as 

, -1. 


u = {I - ml^{z)T) w. 


Consequently, we have 


< 


(1 - ml^{z)r) 


-1 




|w||oo := T(z)llwl 


(2.45) 


Then for z G D(A^, k, £ 2 ), using (1.6) and Proposition A.2 (ii) in [12] (with i5_ = 1 and 9 > c), we can get 

r{z) = 0{\ogN). (2.46) 

Plugging (2.46) and (2.44) into (2.45) yields 

||u||oo ^ Odlu]]^ +0(maxlA|)) ^ ||u||L + 
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where the second step follows from (2.32). Then (2.36) further implies that 

1 

U DO ^ 


y/NT]' 


which together with (2.42) and (2.32) also implies 


\m{z) - msc{z)\ < 


y/Nrj 


Hence 


Ad < ||u||oo + \m{z) - msciz)\ -< 


y/mj' 


Therefore, we completed the proof of Theorem 1.9. 


Proof of Lemma 2.5. For simplicity, we omit the variable z from the notation below. At first, we recall 
the elementary identity by Schur’s complement, namely, 

1 


Gu = 


- (hf )*G«hf 


(2.47) 


where we used the notation to denote the i-th column of H, with the z-th component deleted. Now, 


we use the identity for a,b ^ i (see Lemma 4.5 in [12] for instance), 

Gal = Gat - Ga^GUGu)-^ = Gat " GayG,b[hu -Z- (hf ^ )*G(*) hf ^. 


(2.48) 


By using (1.10) and the large deviation estimate for the quadratic form (see Theorem C.l of [12] for 
instance), we have 


(h«).GWh«-^a2,.GW|^ J 

a^i 


maxlGi*2l2+ maxlG^ji , 

IVl a a^b 


y{i)l 


which implies that 

(h«)*GWhf> 


a^i 


-< max JG^ I + \/t 7 max JgS + max < 3 max JG^^^ ], 


M af^i 


a^b 


a.,b^i 


(2.49) 


(2.50) 


where we have used the fact that = 1 in the hrst inequality above. Plugging (1.20) and (2.50) 

into (2.48) and using Corollary 2.2 we obtain 


|GSHl + ^(l + 3majlG«l) 


which implies 


max I 

a,b^i 


maxlGi;^] ^ 1, 

a,b^i 


(hd>)*GWhf> 


^ 1 . 


In addition, (1.20), (2.48) and (2.51) lead to the fact that 


\Gat{z) - G^^l{z)\ -< —, \G^ah^)\ ^ ^at + 


Nr]' ' ' y/mj' 

Now, using (2.31), (2.47), (2.49) and (2.52), we can see that 

lAI = \-hu + (h«)*GWh« - ^aAG,,| ^ 

a 

Therefore, we completed the proof of Lemma 2.5. 


y a, h i. 


y/mj' 


(2.51) 

(2.52) 

(2.53) 

□ 
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2.3. Proof of Theorem 1.11. With Theorem 1.9, we can prove Theorem 1.11 routinely. At first, due 
to Definition 1.6 and the fact that Gab{z) and msc{z) are Lipschitz functions of z with Lipschitz constant 
it is easy to strengthen (1.16) to 

max sup \Gab{z) - 5ab'msc{z)\ -< 

zGD(]V,K.e2) yAl? 

which implies that 


max sup \Gaa{z)\ -< G (2-54) 

“ zGD(Af.K,e2) 

for some positive constant G due to the fact that msc{z) ^ 1. Recalling the normalized eigenvector 
Ui = {uii ,..., UiAf) corresponding to Aj, and using the spectral decomposition, we have 


maxlmG'aa(-2:) 


= max 


N 

E' 


^ |A,-i?| 2+^2 


N 

E 

2=1 


lloV 




(2.55) 


For any |Ai| < v^— k, we set E = \i on the r.h.s. of (2.55) and use (2.54) to bound the l.h.s. of it. Then 
we obtain 


||2 
11 oo 


V 


-< 1 . 


Choosing r] = N above and using the fact that £2 can be arbitrarily small, we can get (1.17). Hence, 
we completed the proof of Theorem 1.11. 


3. Supersymmetric formalism and integral representation for the Green’s function 

In this section, we will represent E\Gij{z)\'^'^ for the Gaussian case by a superintegral. The final 
representation is stated in (3.30). We make the convention here, for any real argument in an integral 
below, its region of the integral is always R, unless specified otherwise. 


3.1. Gaussian integrals and superbosonization formulas. Let 4> = ((/>i,..., 0^)' be a vector of 
complex components, xj} = {tpi,... ,ipk)' be a vector of Grassmann components. In addition, let </>* and 
1 / 1 * be the conjugate transposes of (p and ip, respectively. We recall the following well-known formulas 
for Gaussian integrals. 


Proposition 3.1 (Gaussian integrals or Wick’s formulas). 

(i) Let A be a k X k complex matrix with positive-definite Hermitian part, i.e. ReA > 0. Then for 
any £ £ N, and ii,... ,i£, ji,..., ji £ {1,..., k}, we have 



dRec^gdlm^g 

TT 


e 

exp{-</.*A</.} 

b=l 


1 

det A 


E )a.v( 6 ) 


creP(^) b=l 


(3.1) 


where P(£) is the permutation group of degree t. 

(ii) Let B he any k x k matrix. Then for any I £ {0,..., kf, any £ distinct integers ii,... ,i(, and 
another £ distinct integers ji, ..., £ {1,..., k}, we have 


„ fc i 

/ Wdiffadipa exp{-'j/)*B'0} = (-l)^+2:i=ih<.+J<») detB('l-'\ 

a=l 6=1 


(3.2) 


where I = {ii,.. .,ie}, and J = {ji,.. 
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Now, we introduce the superbosonization formula for superintegrals. Let x = (Xij) be an £ x r matrix 
with Grassmann entries, f = (/y) be an £ x r matrix with complex entries. In addition, we denote their 
conjugate transposes by %* and f* respectively. Let be a function of the entries of the matrix 


5(f,r;x,x*) 


f*x f*f / 


Let A{x,X*) be the Grassmann algebra generated by x^’s and Xy’s. Then we can regard F as a 
function defined on a complex vector space, taking values in A{x^X*)- Hence, we can and do view 
F(5(f, f*; X, X*)) as a polynomial in Xy’s and Xy’s, hr which the coefficients are functions of /^’s and 
/ij’s. Under this viewpoint, we state the assumption on F as follows. 


Assumption 3.2. Suppose that F(5(f, f* ; Xj X*)) ® holomorphic function of fij’s and fij’s if they 

are regarded as independent variables, and F is a Schwarz function of Re fij’s and \rr\ fij’s, by those we 
mean that all of the coefficients o/F(iS(f, f* ; Xi X*))? as functions of fij’s and fij’s, possess the above 
properties. 


Proposition 3.3 (Superbosonization formula for the nonsingular case, [17]). Suppose that F satisfies 
Assumption 3.2. For £ > r, we have 


I Px f*f dA(x)dz>(y)du;d| ^ ^ ) 


det y 


der (x — wy-i^) 


(3.3) 


where x = [xij) is a unitary matrix; y = {yij) is a positive-definite Hermitian matrix; io and ^ are two 
Grassmann matrices, and all of them are r x r. Here 

df = n = n dXii, 


db(y) = l(y > 0) J^dy** dRey^fedlmy^fc, dwd| = duijdfij, 

i—1 j>k 


and d/i(-) is defined by 


r(r-l)/2 ’’ J 

d/i(x) = —— • n ■ dy(u). 


n; 1 ^! 27ri 

iii=i 

under the parametrization induced by the eigendecomposition, namely, 


X = I/*xU, X = diag(a;i,..., Xr-), VGU{r). 

Here d/i(U) is the Haar measure on U{r), and A(-) is the Vandermonde determinant. In addition, the 
integral w.r.p.t. x ranges over U(2), that w.r.p.t. y ranges over all positive-definite matrices. 


For the singular case, i.e. r > i, we only state the formula for the case of r = 2 and i = 1, which is 
enough for our purpose. We can refer to formula (11) in [3] for the result under more general setting on 
r and i. 


Proposition 3.4 (Superbosonization formula for the singular case, [3]). Suppose that F satisfies Assump¬ 
tion 3.2. If r = 2 and (■ = !, we have 


/ 



X*X 

f*x 




dwd/i(x) ■ l{y > 0)dy • dwd^ F 




ww* \ y(y-|x 
yww* J det^ x 


(3.4) 
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where y is a positive variable; x is a 2-dimensional unitary matrix; us = and ^ = (^ 1 ,^ 2 ) o,re 

two vectors with Grassmann components. In addition, vs is a unit vector, which can be parameterized by 

w = (^l/\/l + |wp,w/\/l + |uip^ , wGC. 

Moreover, the differentials are defined as 

dw = ] ,„.„ dRewdlmw, dwdi = TT dw^d^^. 

(1 + kr) 


In addition, the integral w.r.p.t. x ranges over U{2). 

In our discussion, for w, we will adopt the parametrization 

V = l/\/l + IwP, e^^=w/\w\, n € I, 0 gL 
for convenience. Accordingly, we can get 


dw = vdvd0. 


3.2. Initial representation. For a = 1, 2 and j = 1,. .., W, we set 

*^a,j — (0a,j,l5 ■ ■ ■ : j ■ - ■ ; 

= «,i, . . . , K.w)'^ = «.l, ■ • ■ > K.wY- 

For each j and each a, ^ is a vector with complex components, and d!aj is a vector with Grassmann 
components. In addition, we use d>*^ and to represent the conjugate transposes of $a,j and d/aj 
respectively. Analogously, we adopt the notation $* and di* to represent the conjugate transposes of 
and respectively. We have the following integral representation for the moments of the Green’s 
function. 


Lemma 3.5. For any p, q = 1,. .. ,W and a, (3 = I, ..., M, we have 

I ^P(J,ck/3 ('2^) I — j ddiddl 0(4^2,p,a4^2,q,p') 

X exp |i«'i(z - + i$*(z - - idfUz - H)'^2 - - Ff)$2}, (3.5) 

where 


M 

«= n nn 

a—1,2 j=l oi' — l 


dRe(j)a,j,cy.'d\m4>a 


W M 

d^' = nnn dliSa,j,a'dfia,j,G 

0—1,2 j—1 = l 


Proof. By using Proposition 3.1 (i) with £ = n and Proposition 3.1 (ii) with .£ = 0, we can get (3.5) 
immediately. □ 


3.3. Averaging over the Gaussian random matrix. Recall the variance profile S in (1.2). Now, we 
take expectation of the Green’s function, i.e average over the random matrix. By elementary Gaussian 
integral, we get 

If ^ 

]E|Gpq,a/3(-2:)P" = d^ddl {fil,q,p(t>l^p^a(l>2,p,a4>2,q,pY QXp ^l'^{TrXjJZ + TrYjJZ)^ 

X exp '^Y^SjkTrXjJXkJ — SjkTrYjJYkj'^ 

j,k j,k 

X exp I - ^^SjfcTrfijJSfcjj, 

j,k 


(3.6) 
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where 

J = diag(l, -1), Z = diag( 2 :, z), 

and for each j = 1,, W, the matrices Xj, Yj, Clj and Sj are 2x2 blocks of a supermatrix, namely, 


4 

Clj \ 



1112.,"JllD 

^4^2., 


b J- 






5,= 


Remark 3.6. The derivation of (3.6) from (3.5) is quite standard. We refer to the proof of (2.1)) in 
[19] for more details and will not reproduce it here. 

3.4. Decomposition of the supermatrices. From now on, we split the discussion into the following 
three cases 

• (Case 1): Entries in the off-diagonal blocks, i.e. p Y Qj 

• (Case 2): Off-diagonal entries in the diagonal blocks, i.e. p = q, ct ^ (3, 

• (Case 3): Diagonal entries, i.e. p = q, a = [3. 

For each case, we will perform a decomposition for the supermatrix Sj (j = p ov q). For a vector v and 
some index set I, we use to denote the subvector obtained by deleting the i-th component of v for all 
i G I. Then, we adopt the notation 


sj" = 


4 

% 


4 I) 




4 


^3 

y}\ 


Here, for A = Xj,Yj, Qj or :^j, the notation is defined via replacing ^a,j, ^a,j, and j by 


respectively, for a = 1,2, in the definition of A. In addition, the notation A1®1 
is defined via replacing ^a,j, and by ([a.j.i, Ya,j,i, and ifa.j.i respectively, for a = 1, 2, 

in the definition of A. Moreover, for A = Sj, Xj, Yj, Clj or Sj, we will simply abbreviate and 

AH“}> by /\{a,b) aY\ respectively. Note that is of rank-one. 

For Case 1, due to symmetry, we can assume a = /3 = 1. Then we extract two rank-one supermatrices 
from Sp and Sq such that the quantities (j) 2 ,p,i 4 >i,p,i and 4 >i,q,i 4 > 2 ,q,i can be expressed in terms of the 
entries of these supermatrices. More specifically, we decompose the supermatrices 


Consequently, we can write 




— (Yq^^) 12(Y^'^^) 2 


(3.7) 


(3.8) 


q 

For Case 2, due to symmetry, we can assume that a = 1, /3 = 2. Then we extract two rank-one 
supermatrices from Sp, namely, 

4 = 41’2) ^5[2]_ (3 9) 

Consequently, we can write 

(>1,P,24>1,p,i4>2,p,1<['2,p,2 = ( 4 ^^ ) 12 ( 4 ^^ )21 • (3.10) 

Finally, for Case 3, due to symmetry, we can assume that a = 1. Then we extract only one rank-one 
supermatrix from Sp, namely. 


4 = 411 + 54 . 


( 3 . 11 ) 
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Consequently, we can write 

<^1,P,1'/'1,P.1^2,p,1?!>2,p,1 = ) 11 )22 • 

Since the discussion for all three cases are similar, we will only present the details for Case 1. More 
specifically, in the remaining part of this section and Section 4 to Section 10, we will only treat Case 1. 
In Section 11, we will sum up the discussions in the previous sections and explain how to adapt them to 
Case 2 and Case 3, resulting a final proof of Theorem 1.14. 


3.5. Variable reduction by superbosonization formulae. We will work with Case 1. Recall the 
decomposition (3.7). We use the superbosonization formulae to reduce the number of variables. We shall 
treat Sk {k ^ p, q) and 5^ (j = p, q) on an equal footing and use the formula (3.3) with r = 2,i = M 
for the former and r = 2,i = M — 1 for the latter, while we separate the terms {j = p,q) and use 
the formula (3.4). For simplicity, we introduce the notation 

{ Sj, if jy^P,q, 

(3.12) 

if j=p,q. 

Accordingly, we will use Xj, Clj, and Yj to denote four blocks of Sj. With this notation, we can 
rewrite (3.6) with a = /3 = 1 as 

If n ^ 

E|Gpq,ii(z)p" = y d$d4' (^^i,,,i^f)i,p,i^2,p,i</>2,9,i) expjiy^ (TrXjJZ + TrYjJZ^'^ 

^ ^ ^ TVfe J - T rYj Jffc j) I exp I - ^ ^ SjfeT rClj JSfc J | 

j,k j,k 

X n exp{iTrXWjZ + iTrFWjz} exp {- jl™ J - Tr^-JF™ j) } 

k=p,q k=p,q j=l 

X n exp{|^(rrlWjxWj-TrFWjFWj)} exp { - JS'V} 

k,i—p,q k,i—p^q 

^ n exp I — ^ y] Sjfe ^ J + Trf2|. ^ I (3.13) 

k=p,q j 

Now, we use the superbosonization formulae, i.e., (3.3) and (3.4), to change to the reduced variables 
as 




[ 1 ] 


->■ X, 


Yj 


[ 1 ] 

k ’ 


n 


dl ^ 








* ^ 


—^ ^j, j — 1,..., IF, 

W ^ wWjfW fW ^ F™ — nWwWfwWF k-n n 


(3.14) 


Here, for j = 1,..., IF, Vj is a 2 x 2 unitary matrix; F,- is a 2 x 2 positive-definite matrix; ilj = {ujj^ap) 
and Sj = {^j^ap) are 2x2 Grassmann matrices. For k = p oi q, X^^^ is a 2 x 2 unitary matrix; is a 
positive variable; = (w[,^j,is a column vector with Grassmann components; = (^j;^i,Cfe 2 ) 
is a row vector with Grassmann components. In addition, for k = p,q, 


[1] - 


wL' = 




du 





el, 




e L. 


Then by using superbosonization formulae, we arrive at the representation 


(3.15) 


E|Gp,,n(z)p" 
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(n 


lyllZ+4 I dXWdyWdwWda;Wd|WdXdFdf]dS (yWy[i](wW(wW)*)^^(wW(wW)*)^^y 


w 

X exp Ii ^ (^TrXj JZ + TrYj JZ^ | ^ X! J - TrY^ JYu j) | 


j=i 


X exp < — 


j,k j 


j,k 

det“ Yi 


n expji 


! iTrXl^bz + iTtY;^ 

^ n exp{|^(TrXWjxWj-rryWjrWj)} n exp{-^Tru;W(wW)VwW|Wj| 

k^i—p^q k,£—p,q 

W 1 ^ 

n n exp { - — ^%feTra;W(wW)*JS, j} 


di] 


det^{X,-n,Y-%) 

1 ^ 


n 


det (Xfe — ^kY)^ ^Sfe) 


det Yfc 


k=p,q 


Jz] n exp{-^s,fc(rrX,JX^ 


[ 1 ] 


, j-Trr,jrWj 


k—p,q 




)} 


k=p,q 

X n 

k=p,q 


4 = 1 


k=p,q 


4 = 1 


det^(xW) 


(3.16) 


where we used the notation yPl := w^l := (wp^,Wq^^). The differentials in (3.16) are defined 

by 

dXPl := H dA(xW)d/i(xW), dybl := [] l(yW > 0)d2/W, 


j=p,q 


j=p,q 


dwP^ := , 


j=p,q 

w 


du.wd«w := n n difji 

a=l,2j=p,q 


W 


W 


dX := d/i(Xj), dr := dP(Y,), dnd.=, := dwj^a^gd^j,^;?. 

j—1 j—1 a,j8—l,2j—l 

The regions of the integral of Xj’s and are all U(2), and those of r, ’s are the set of all positive- 

definite matrices. The integral of ranges over I and that of ranges over L, for k = p,q. 

Now we change the variables as XjJ —>■ Xj, YjJ —>■ Bj, VljJ —>■ Hj,E.jJ —>■ Ej and perform the scaling 
Xj — >■ —MXj, Bj — ^ MBj, rij —>■ and Ej —>■ y/MEj. Consequently, we can write 

E|Gp,.ii(z)p"= J dXPldyWdwWdu;Wd|WdXdBdfldS exp { - M(if (X)-f T(B))} 

xr{n, S, X, B) ■ Q(0, S, u;P 1 , |d], xPl, yPl, wdl) • X(X, B, XPl, yPI, wPl), (3.17) 

where the functions in the integrand are defined as 

K{X) -=-1^2 SjkTrXjXk + iEj2 TrX, + ^ log det X,, 

4,fe 4 4 

L{B) := ^'^SjkTrBjBk - 'iE'^TtBj - ^ log det 


j,k 


V{n,E,X,B) :=exp{-^s,feTrfljSfc} 


1 


j,k J detail+ M-^X-^n,B-%) 

det (Xfc -|- M ^VlkB^ ^^k) 


n 

k=p,q 


det Bk 


X 
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X n exp I - 

k=p,q 

X n exp I - 

k^£—p,q 


,wW):= n 


J-(X,i?,xW,yW,wW) := /(X,^™) g(i3,yW,wW), 


(3.18) 


with 


w 

/(X,xW):=exp{M, 7 ^rrX, j} J] 


i=i 


k,^—p,q 


1 


,n, det^(xW) 


■ exp {iTrX™ JZ - ^ s,feTrX,xW j}, 


(3.19) 


g(B,yW,wW) :=exp|-M?7^rrBj j| exp | - I^TrF^ Jf/’ j| 

j = l k,l—p,q 

X ((wW(wW)*)^^(wW(wW)*)^y n (yr’)”+"exp{irryWjZ-^5,,TrB,rWj}. 


(3.20) 


k—p,q 


In (3.17), the regions of Xj^s and are all 17(2), and those of Bj’s are the set of the matrices A satisfy¬ 
ing AJ > 0. Roughly speaking, here we collected the terms containing the Grassmann variables from ’s 
and l^’s, resulting the factor Q(n,S,a;W,TfW,yW,wW). Then, we put the terms containing the 
Grassmann variables in llj’s and S^-’s together (except those in 2(11, S, u) , wdl)), resulting 

the factor 7^(0, S, X, B). Finally, we sorted out the terms containing rj, and the variables from 7f|,^Ps 
and w^^s (except those in 2(11, S, a;d], X^l, ydl, wdl))^ resulting the factor X{X, B, Tfdl^ ybl ^ w^l). 
This separation indicates the order in which we will perform the integrations. 


3.6. Parametrization for X, B. Similarly to the discussion in [19], we start with some preliminary 
parameterization. At first, we do the eigendecomposition 

Aj = p;AjP„ Bj = Q-^BjQj, Pj e 17(2), Q, e 17(1,1), (3.21) 

where 

Xj = diag(xj,i,a::j,2), Bj = diag(&j,i,-6j-,2), Xj^i,Xj^2 £ S, l>j,i,&j,2 G K-i-- (3.22) 

Further, we introduce 

Vj=PjP: €11(2), T,=Q,Q^^ €U{1A), j = 1,...,W. (3.23) 


Especially, we have Fi = Ti = I. Now, we parameterize Pi, Qi, Vj and Tj for all j = 2,. 



u = — Uj = ^1 — V, Vj G I, 6 ,9j € L, 


W as follows 


Qi — 


s 


te-“" 
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— \/ 1 “t" 1 -h ij ^ ^ ^+5 ^ 


(3.24) 


Under the parametrization above, we can express the corresponding differentials as follows. 

w w 


dXdB = d^,{Pi)dv{Qi) ■ nw)d-(T,).n dbj^idbj^2 ■ 


dxj,i dxj^2 
27 ri 27 ri 


w 


x2^(7r/2)^^ Y[{xj^i - Xj-2)^(6 j.i + bj^2f, 
t=i 


(3.25) 


where 


d^(Pi) = 2'ydw • d/r(Vj) = 2z;,di;i • diy{Qi) = 2tdt ■ di^(T,) = 2UdU • 

27r 27r 27r 27r 

In addition, for simplicity, we do the change of variables 

Uj ^ Pi OjQi, Sj ^ QC^SjPi. (3.26) 

Note that the Berezinian of such a change is 1. After this change, P(U, S, A, P, yl^l, wl^l) turns out to 
be independent of Pi and Qi. 

To adapt to the new parametrization, we change the notation 

KiX)^KiX,V), LiB)^L{B,T), P(U, S, A, P) ^ P(U, S, A, P, U, T), 

Q(U,S,a;W,, A™,y™, wW) ^ Q(r!,S,A, Qr,yW, w™), 

P(A, P, All], y[11, will) ^ P(A, P, U, T, Pi , Qi , A™ , yW , w™), 

/(A,aW) ^ /(Pi,U,A,aW), 5(iJ,yW,wW) ^ g(Qi,r, Ay''',wW). (3.27) 

We recall here that K{X) does irot depend on Pi, as well, L{B) does irot depend on Qi. Moreover, 
according to the change (3.26), we have 

1 


P(U, S, A, P, U, T) = exp { - ^ SjkTrn.Ek] ■ [] 


, det^ {l + M-W;X-^V,n,T-^B-^T,E,) 


X n 

k=p,q 


det (U^AfeUfc + M-iUfeT-ip-iPfeSfc) 


det Pfc 


(3.28) 


and 


Q(U,S,u;[iU[il,Pi,Qi,AW,yW,wW) 


n exp{ - ^^s,.(TrPi*U,QiwW|Wj + Tru;[il(wW)VQr'S,T^i)} 
1 


k—p,q 


X n (l - , ,3.29) 


k,i—p^q k—p^q 

Consequently, using (3.25), from (3.17) we can write 


E|G„.u(^)|""= „ / WMVAMTA 

W 

X exp {-M(A(A, V) + L{B, T ))} • [^(a:,-1 - a.,.2)"(&,.i + 6, .2)^ • A(A, P, U, T). (3.30) 
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where we introduced the notation 

A{X, B, V, T) := J dxWdyWdwWdwWd|Wdl]dSd/i(Pi)dj2(Qi) V{n, E, X, B, V, T) 

X J-(X,B,y,T,Pi,Qi,xW,yW,wW). (3.31) 


In (3.30), the regions of Vj^s are all 17(2), and those of Tj’s are all 17(1,1). Observe that all Grassmann 
variables are inside the integrand of the integral A(X, B, V, T). Hence, (3.30) separates the saddle point 
calculation from the observable A(X, B, V, T). 

To facilitate the discussions in the remaining part, we introduce some additional terms and notation 
here. Henceforth, we will employ the notation 

(A'I'l)-' = {(AW)-‘,(Am)-‘}, = {(!/fy,(sW)-'} 

for the collection of inverse matrices and reciprocals, respectively. For a matrix or a vector A under 
discussion, we will use the term A-variables to refer to all the variables parametrizing it. For example, 
Tfj-variables means Xj^i and Xj^ 2 , and 7f-variables refer to the collection of all Tfj-variables. Analogously, 
we can define the terms T-variables, y^^l-variables , fl-variables and so on. We use another term A-entries 
to refer to the non-zero entries of A. Note that A"j-variables are just Tfj-entries. However, for Tj, they 
are different, namely. 


Tj-variables : vs. Tj-entries : Sj, ‘'^^7 


Analogously, we will also use the term T-entries to refer to the collection of all Tj-entries. Then H-entries, 
wl^l-entries, etc. are defined in the same manner. It is easy to check that Q^^-entries are the same as 
Qi-entries, up to a sign, as well, T^“^-entries are the same as Tj-entries, for all j = 2,..., W. 

Moreover, to simplify the notation, we make the convention here that we will frequently use a dot to 
represent all the arguments of a function. That means, for instance, we will write P(H, S, X, B, V, T) as 
P(-) if there is no confusion. Analogously, we will also use the abbreviation Q(-), A(-), and so on. 

Let a := {oi,..., ae} be a set of variables, we will adopt the notation 


i3(a; Ki,/t2,K3) 

to denote the class of all multivariate polynomials p(a) in the arguments oi,..., such that the following 
three conditions are satisfied: (i) The total number of the monomials in p(a) is bounded by ki; (ii) the 
coefficients of all monomials in p(a) are bounded by K 2 in magnitude; (iii) the power of each in each 
monomial is bounded by /ca, for all* = !,...,£. For example, 

5b-j + 3bj,it] + 1 G 3,5,2). (3.32) 

In addition, we define the subset of 0(a; ki, ^ 2 , w^a), namely, 

0deg(a;Ki,«:2,«:a) 

consisting of those polynomials in 0(a; ki, K 2 , ^a) such that the degree is bounded by Ka, i.e. the total 
degree of each monomial is bounded by K 3 . For example 


bbj l + Sbj^it'^ + 1 G 0deg({&j_i, 3, 5, 3 ). 
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4. Preliminary discussion on the integrand 


In this section, we perform a preliminary analysis on the factors of the integrand in (3.17). For 
convenience, we introduce the matrix 


a = 




(4.1) 


4.1. exp{—M(if (X, V)-\-L{B^ 2^))}- Recall the parametrization of Bj, Xj, Tj and Vj in (3.22) and (3.24), 
as well as the matrices defined in (1.24). According to the discussion in [19], there are three types of 
saddle points of this function, namely, 

• Type I : For each j, {Bj,Tj,Xj) = {D±,I,D±) or 

dj € L, Vj = 0 if Xj = Xi , and Vj = 1 if Xj ^ Xi. 

• Type II ; For each j, {Bj,Tj,Xj) = {D±,I,D+) and Vj G U{2). 

• Type III ; For each j, {Bj,Tj, Xj) = (i4±,/,£>_) and Vj £ 17(2). 

(Actually, since 0j and Vj vary on continuous sets, it would be more appropriate to use the term saddle 
manifolds.) Note that at each type of saddle points, we have {Bj,Tj) = (D±,I) for all j. We will see that 
the main contribution to the integral (3.17) comes from some small vicinities of the Type I saddle points. 
Furthermore, the contributions from all the Type I saddle points are the same, which can be explained 
as follows. At first, by the definition in (3.23), we have Vi = I. If we regard 9j’s in the parametrization 
of Vj's as fixed parameters, it is easy to see that there are totally 2^ choices of Type I saddle points. 
Moreover, if Vj =1, we can do the transform 

Xj^3Xj3 = Xi, Vj^3Vj=I in 11 (1,1). 


Consequently, it suffices to consider two saddle points 

(B,,T„Xj,Vj) = (D±,I,D±,I), or (D±,I,D^,I), (4.2) 

corresponding to Xi = D± or Dip, respectively. Furthermore, the contributions to the integral (3.17) 
from the vicinities of these two saddle points are also the same. To see this, we recall the fact that the 
original integrand in (3.17) is a function of the entries of Xj = Pf^XjPj. Now we do the transform 
Pj —>■ 3Pj and Xj 3Xj3 for all j = 1,..., IF to change one saddle in (4.2) to the other. Now, since the 
Haar measure on U(2) is invariant under the shift Pi —>■ 3Pi , the integral over Pi-variables is unchanged. 
That means, for Type I saddle points, it suffices to consider 

• Type!’: For each j, (Bj,Tj,Xj,Vj) = (D±,I,D±,I). 

In summary, the total contribution to the integral (3.17) from all Type I saddle points is 2^^ times that 
from the Type T saddle point. 

Following the discussion in [19], we will show in Section 5 that both K(X, V) — K(D±, I) and L(B, T) — 
L(D±,I) have positive real parts, bounded by some positive quadratic forms from below, which allows 
us to perform the saddle point analysis. In addition, it will be seen that in a vicinity of Type T saddle 
point, exp{—M(Ar(X, F) + L(B,T))} is approximately Gaussian. 


4.2. Q(n, S, Fi, Qi, y[^l, wl^l). The function Q(-) contains both the 11, S-variables from 

F(-), and the Pi, Qi, X'3\yV\ , -variables from P(-). In addition, note that in the integrand in (3.17), 
Q(-) is the only factor containing the u}7] and variables. Hence, we can compute the integral 


Q(H,S,Pi,gi,AW,yW,wW) := J Q(0, S, Ri, Qi, w™) (4.3) 


at first. The explicit formula for Q(-) is complicated and irrelevant for us. From (3.29) and the definition 
of the Grassmann integral, it is not difficult to see that Q(-) is a polynomial of the (AI^I)"^, (yt^^)”^, 
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w^^l, Pi, Qi, ri and S-entries. In principle, for each monomial in the polynomial Q(-), we can combine the 
Grassmann variables with P(-), then perform the integral over and S-variables, whilst we combine the 
complex variables with P(-), and perform the integral over yl^l, w^^l. Pi and Qi-variables. A formal 
discussion on Q(-) will be given in Section 6.1. However, the terms from Q(-) turn out to be irrelevant 
in our proof. Therefore, in the arguments with Q(-) involved, a typical strategy that we will adopt is as 
follows: we usually neglect Q(-) at first, and perform the discussion on P(-) and P(-) separately, at the 
end, we make necessary comments on how to slightly modify the discussions to take Q(-) into account. 


4.3. P(H,S,A,P,y,T). We will mainly regard P(-) as a function of the fl and S-variables. As mentioned 
above, we also have some H and S-variables from the irrelevant term Q(-). But we temporarily ignore 
them and regard as if the integral over H and S-variables reads 

P(A, P, V, T) := J dfldE r{fl, S, A, P, V, T). (4.4) 

We shall estimate P(-) in three different regions: (1) the complement of the vicinities of the saddle points; 
(2) the vicinity of Type I saddle point; (3) the vicinities of Type II and III saddle points, which will 
be done in Sections 6.2, 9.1 and 10.1, respectively. (Definition 5.5 gives the precise definition of the 
vicinities.) In each case we will decompose the function P(-) as a product of a Gaussian measure and 
a multivariate polynomial of Grassmann variables. Consequently, we can employ (3.2) to perform the 
integral of this polynomial against the Gaussian measure, whereby P(-) can be estimated. Especially, it 
turns out that in the vicinity of Type I saddle points, P(-) is approximately the normalizing constant of 
the Gaussian measure obtained from exp{—M(Ar(A, V) -I- P(P,T))}. 


4.4. P(X,P,l/,T,Pi,Qi,AW,yW , wl^l). Observe that T is the only term containing the energy scale rj. 
As in the previous discussion of P(-)i here we also ignore the Pi, Qi, A[^l,y[^l, wl^l-yariables from the 
irrelevant term Q(-) temporarily, and investigate the integral 


F(A, P, V,T) = J dAWdy[ildwWd^(Pi)di.(Qi) P(A, P, V, T, Pi, Qi, A^, yW, w™) 

= J dAWd;i(Pi)/(A,F,Pi). J dyWdwWdi.(Qi)g(P,T,Qi,yW,wW). (4.5) 

We shall also estimate F(-) in three different regions: (1) the complement of the vicinities of the saddle 
points; (2) the vicinity of Type I saddle point; (3) the vicinities of Type II and III saddle points, which 
will be done in Sections 6.3, 9.2 and 10.2, respectively. 

Especially, when we restrict the A, P, V and T-variables to the vicinity of the Type I saddle points, 
the above integral can be performed approximately, resulting our main term, a factor of order 1 / (N. 
This step will be done in Section 9. It is instructive to give a heuristic sketch of this calculation. At first, 
we plug the Type I saddle points into (4.5). We will show that the integral of /(•) approximately reads 

g-(a+-a_)JV,y f(D±,I,P,) ^ 

which is the easy part. Then, recalling the definition of g{-) in (3.20) and the parameterization (3.15), 
we will show that the integral of g{-) approximately reads 


g(o+ a^)Nfj J dy[i]dw[^l(izz(Qi) 5(P±,/,(5i,y[^l,w[^l) 

pOO p [1] [1] 

- / 2 tdt / dffWdaW * 

Jo 

pOO 1 

- / 2 tdt ■ - 


0 




(4.6) 
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where in the second step above we used the fact that 

[ da ■ ~ T. 

Jh 

We notice that the factor ' in (4.6) actually comes from the term 

((w™ (wW)*) ,2 (w™ (w™)*)2l) 

in (3.20). This factor brings a strong oscillation to the integrand in the integral (4.6). In Case 2, an 
analogous factor will appear, resulting the same estimate as (4.6). However, in Case 3, such an oscillating 
factor is absent, then the estimate for the counterpart of the integral in (4.6) is of order 1/Nr] instead of 
l/(iV 77 )”+^. The detailed analysis will be presented in Sections 10 and 11. 


5. Saddle points and vicinities 


In this section, we study the saddle points of K(X,V) and L{B,T) and deform the contours of the 
H-variables to pass through the saddle points. Then we introduce and classify some small vicinities of 
these saddle points. The derivation of the saddle points of K{X,V) and L{B,T) in Section 5.1 and 5.2 
below is essentially the same as the counterpart in [19], the only difference is that we are working under 
a more general setting on S. Hence, in Section 5.1 and 5.2, we just sketch the discussion, list the results, 
and make necessary modifications to adapt to our setting. In the sequel, we employ the notation 


ba ■— (^1,05 - ^ ■ — (^l,a : • • • : ^ W , a ') : a — 1 5 2, 

t := V := { v2 ,..., vw ), a := { a2 , ■ ■ ■ , aw ), 6 := { 62 ,. ■. , 0w )- (5.1) 


As mentioned above, later we also need to deform the contours, and discuss the integral over some 
vicinities of the saddle points, thus it is convenient to introduce a notation for the integral over specific 
domains. To this end, for a = 1, 2, we use l[] and 1“ to denote generic domains of and respectively. 
Analogously, we use I* and I” to represent generic domains of t and v, respectively. These domains will 
be specified later. Now, for some collection of domains, we introduce the notation 



X exp {-M(K(A, V) + LiB, T ))} • f[(x,,i - x,, 2 )^( 6,.1 + 6 ,, 2 )" • A(l, B, H, T). (5.2) 

j=i 

For example, we can write (3.30) as 


E|Gp,.ii(z)p" 










(5.3) 


which is the integral over the full domain. 

5.1. Saddle points of L{B,T). We introduce the function 

—iAa —logo, a G C. (5.4) 

Recalling the definition of i(-) in (3.18), the decomposition of Hj’s in (3.21) and the definition of T^’s in 
(3.23), we can write 

= -\Y.^rkTr{T-^B,T, - T-^B.T,)^ + ^ (^TrB^ - iETrB, - logdet^) 

j,k j 
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= ~ ~ - logdet 4 ) 

j,k j 

+ 2 ^)i2p(&j,i + bj^2){bk,i + bk,2) 

j,k 

=:£{h,)+i{-h 2 )+£s{B,T), (5.5) 

where we used the notation introduced in (5.1), and the functions £(■) and £si') are defined as 

£(a) := -^'^Sjk{aj - Qkf + '^k{aj), a= (ai,...,av^) G C’^, 
j,k j 

£s{B,T) := i^s,fe|(TfcT-i)i2p(6,,i + 6,.2 )(Vi +^fc,2)- (5-6) 

i.fe 

Following the discussion in [19] with slight modification (see Section 3 therein), we see that for |i?| < 
■\/2 — K, the saddle point of L{B, T) is 

(4,T,) = (Z1±,/), Vj = 1,...,VF, (5.7) 

where D± is defined in (1.24). For simplicity, we will write (5.7) as {B,T) = {D±,I) in the sequel. 
Observe that 


L{D±J) = £{a+)+£{a.), £s(77±,/)=0. (5.8) 

We introduce the notation 

i++(a) := ^(a) - £(a+), i+-(a) := £(a) - £{a-) L_(a) := £(-a) - £(a_), (5.9) 

where £{a+) represents the value of £{sl) at the point a = (a+,..., a+), and £{a-) is defined analogously. 
Correspondingly, we adopt the notation 

LiB,T) :=LiB,T)-L{D±,I)=£++(hi)+L.ih 2 )+£siB,T), (5.10) 

which is implied by (5.5), (5.8) and (5.9). Now, for each j = 1,..., W, we deform the contours of 
and bj ^2 to 


bj^i G F := {ra+jr G K+}, &j ,2 G F = {—ra_|r G M+} 


(5.11) 


to pass through the saddle points of i3-variables, based on the following lemma which will be proved in 
Section 7. 


w toW/ =E|Gpq,ii(z)|^ 


Lemma 5.1. With the notation introduced in (5.2), we have 

'Tfy^W f^W _ 7" 

X (^i , i , X , Zj , K_|_ , ii ) ~ ’ ■^'^+ ’ 

We introduce the notation 

rj,i = \bj,i\, rj^2 = \bj,2\, j = I,... ,W. (5.12) 

Along the new contours, we have the following lemma. 

Lemma 5.2. Suppose that \E\ < \/2 — k. Let bi G F'^, b 2 G f^. We have 

w w 

ReL(B,T) >cJ2E " 1)' + - 1)) + Re£s(B,T) > c ^ ^(r,., - 1)^ 

a—l,2j—l a—l,2j—l 

(5.13) 


for some positive constant c. 
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Proof. Since \E\ < \f2 — k, we have Re( 6 j 4 + bj^ 2 ){bk,i + bk, 2 ) > 0 for bi G and b 2 € thus 
Re£s{B,T) > 0, in light of the dehnition in (5.6). Consequently, according to (5.10), it suffices to prove 

w 

Re£++(bi) + Rei—(b 2 ) > c ^ ^ - 1)^ + (j’yo - logr^.a - 1)) (5.14) 

0 — 1,2 j = l 

for some positive constant c. To see this, we observe the following identities obtained via elementary 
calculation, 

Re£++(bi) = —-—(2 “ 1)^) + “ logr^, 1 - l) 

j 

Rei__(b2) = ——(2 - rk,2f - - 1 )^J + {rj,2 - logrj-2 - l), 

i.fc j 1=1 

which together with \E\ < \f2 — n and (1.4) implies (5.14) immediately. Hence, we completed the proof 
of Lemma 5.2. □ 

5.2. Saddle points of K{X,V). Analogously, recalling the definition in (5.6), we can write 

K{X, H) = i ^ [\TrX^ - iETrX, - logdet A,) 

i.fc j 

= ^ - E - logdet A,) 

i.fc 3 

+ ^Y^jk\iVkV*)i2\'^{Xj,i - Xj,2){xk,l - Xk,2) 
j,k 

= -£(xi)-^(x 2 ) + £s(A,H), (5.15) 

where £{■) is defined in the first line of (5.6) and €s(A, V) is defined as 

is{X,V) = - 2 ^j. 2 )( 2 :fc,i -a:fc. 2 ). (5.16) 

3,k 

Analogously to the notation L{D±,I), we will use K{D±,I) to represent the value of K{X,V) at 
{Xj,Vj) = {D±,I) for all j = 1,...,IT. In addition, K{D+,I) and K{D-,I) are dehned in the same 
manner. Observing that 

es{D±,I)=is{D+,I)=isiD_,I) = 0, (5.17) 

we have 

A(Z)±,/) =-f(a+)-£(a_), A(Z)+,/) =-2f(a+), A(Z)_, J) =-2£(a_). (5.18) 

Moreover, we employ the notation 

A(A, H) = K{X, H) - K{D±,I) = -i++(xi) - £+_(x 2 ) + £s{X, H). (5.19) 

We will need the following elementary observations that are easy to check from (5.18) and (5.6) 

A:(i:'±,/) + L(i:>±,/) =0, ReA(L>+,/) = ReA'(D_,/) = ReA(Il±,/). (5.20) 

In addition, we introduce the W x W matrix 

S^ = {s]k), s], := s,k\{VkV;U\ 


(5.21) 




and the 2W x 2W matrices 
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s = ses, 



(5.22) 


where S" depends on the F-variables according to (5.21). Here we regard H-variables as fixed parameters. 
Due to the fact \iykV*)i 2 \ G I, it is easy to see that S’' is a weighted Laplacian of a graph with 2W 
vertices. In particular, S" < 0. By the definition (5.21), one can see that 5’) = 0 for alH = 1,..., ID. 
Consequently, we can obtain 




Similarly to (1.4), we get 


-Bjj, if j = l,...W 

-Sj-w,j-w, if j = W + l,...,2W 


/ + S" > Co/, 


(5.23) 


where cq is the constant in Assumption 1.1 (ii). Moreover, it is not difficult to see from the definitions in 
(5.16), (5.21) and (5.22) that 

i^s,ferr(A,-Xfe)2+4(A,H) = -ix'S’'x, (5.24) 

j,k 


where we used the notation x := (x'j^jX^)'. 
Now let 


dj = argxyi, -dw+j = argXj, 2 , V j = 1,..., VF. 

Then, recalling the parametrization of V}’s in (3.24), we have the following lemma. 


(5.25) 


Lemma 5.3. Assume that Xj^i,Xj ^2 G D for all j = 1,..., W. We have 

2W p 2 

Rek{X, V)>^ El®') jfc(cosdj — cosdfc) 2+c^ (sinr^j --) (5.26) 

j,k=l j=l 

for some positive constant c. In addition, ReK{X,V) attains its minimum 0 at the following three types 
of saddle points 

• Type I: For each j, Xj = D± or Dip, 

9j G L Vj = 0 if Xj = Xi, and Vj = 1 if Xj Xi, 

• Type II: For each j, Xj = D+, Vj G f/(2), 

• Type III: For each j, Xj = /!_, Vj G U{2), 

which are the restrictions of three types of saddle points in Section 4-1, on X and V-variables. 

Remark 5.4. The Type I saddle points of {X,V) are exactly those points satisfying 

V*X,Vj=D±, Vj = l,...,IT, or V;XjVj = D^, Vj = l,...,IT. 

In Lemma 5.3, we wrote them in terms of Xj, Vj and 9j in order to evoke the parameterization in (3.22) 
and ( 3 . 24 ). 


Proof. By (5.15), (5.24), the definitions of the functions /(•) in (5.7) and k(-) in (5.4), we can write 
K(X,V) = --x*S’'x-^ (ik(xj-i) +k(xy2) -lk(a+) -k(a-))- 
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By using (5.25) and the fact \xj^a\ = 1 for all j = 1,..., VF and a = 1,2, we can obtain via elementary 
calculation 

2W 

ReK{X,V) = 4 ^ (S"),fc(cos^?j - cosAf 

j,k=l 

2W 2W p 2 

“4 (§")jfc(sini?j -sin7?fc)2 + ^ (^siniJj - —) . (5.27) 

j,k=l j=l 

In light of the fact §’' < 0 and (5.23), we have 

/+!§’'>/ +S’'> Co/. (5.28) 

Applying (5.28) to the last two terms on the r.h.s. of (5.27) yields (5.26). 

Now, we show that ReK{X,V) attains its minimum 0 at three types of points listed in Lemma 5.3. 
It is elementary to check that these points are minimum points along the contour by plugging them into 
the definition of ReK{X, V). In the sequel, we show that they are the only solutions to the equation 

ReK{X,V)=0. (5.29) 


At first, by the second term on the r.h.s. of (5.26), we see that for any solution to (5.29), 


sin -dj 


E 
"2 ’ 


Vj = l,...,2VL, 


which implies that Xj^a = or a_ for all j = 1,..., IL and a = 1, 2, by recalling the definition (5.25) 
and the definitions of a+ and a_ in Section 1.4. Consequently, for each j, Xj can only be one of D±, 
Hzp, and /!_. 

Suppose that Xi = £>+, we claim that Xj = /!+ for all j. Otherwise, owing to the fact that the graph 
0 is connected, there exists {i,j} G £ such that Sij > 0 and Xi = D+ but Xj = D±, Dzp or £)_. Without 
loss of generality, we assume Xj = D±. In this case, we use the fact 


Re/C(X,IG) > -(S’')y(cosi?i — cosiJj)^ + -(S’')i+wj(cosi?i+vy — cos’djY 


1 


1 


+ -(S’’)jj+iv(cosi?i - cos-dj+w) + ^(S’')j+wj+w(cosi?i+iv - cos-dj+w) , 

which follows from (5.26) directly. Now, by the assumption Xi = D+ while Xj = D±, we have 


cosdi = cosi?i+w = cos'dj = Re(a+) = 


V4- /;2 


COS1? 


3+W 


= Re(a_) = — 


V4- A2 


(5.30) 


which together with (5.30) implies that 


ReAr(i:, V) > 


A-E^ 




3+W 


(§") 


i+W,3+W 


4-^2 

5ij > 0, 


(5.31) 


contradicting to (5.29). Analogously, we can show that Xj can not be Hip or H_. Consequently, for a 
solution to (5.29), if Xi = /!+, we have shown that Xj = 11+ for all j. Similarly, we can show that if 
Xi = D-, then Xj = D- for all j. These two kinds of solutions are collected as the Type II and Type 
III saddle points, respectively. 

What remains is to show that if Xi = D± or Hzp, the solution to (5.29) must be one of the Type I 
saddle points. We only show the case of Xi = D±. Assume that {!,*} G f in the graph Q, i.e. Si^ > 0. 
At first, similarly to the discussion from (5.30) to (5.31), we can show that Xi can only be D± or Dzp. If 
Xi = D±, then by using (5.30) with j = 1, we have 

ReK{X,V) > + (S’').+rv,i) = > 0, 
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and the equality holds if and only if Vi = I, according to the assumption Vi = I and the definition in 
(5.21). The discussion on the case of Xi = Dzp is analogous. Consequently, we have 

V;x,Vi = Vi*XiVi = D±. (5.32) 

Since the graph Q is connected, we can show that (5.32) holds for all i = 1,...,1T. Analogously, if 
Xi = £>zp, we can show that V*XjVj = Hip for all j = 1,..., W. Therefore, we completed the proof of 
Lemma 5.3. □ 

5.3. Vicinities of the saddle points. Having studied the saddle points of L{B,T) and K{X,V), we 
then introduce some small vicinities of them. To this end, we introduce the quantity 

0 = 0(iV, eo) := (5.33) 

for small positive constant Eq which will be chosen later. Let a = (ai,... ,aw) £ be any complex 
vector. In the sequel, we adopt the notation 

a + d := (ai + d,..., aw + d), da := (doi,..., daw) j V d £ C, 

arg(a) := (arg(ai),..., arg(aw)) • (5.34) 

Now, we define the following domains . 

= T’)_{N,eo) := {a € T'^ : ||a- a+H^ < ^}, 

T'_ = T^(iV,Eo) := {ae : ||a + a _||2 < ^}, 

= T^(7V,Eo) := {as : 11 arg(a;ia)| |i < ^}, 

(iV,Eo) := {ae : 11 arg(alia)||i < ^}, 

Ts = TsiN,eo) := {a e : -a'S^^^a < ^}, (5.35) 

where the superscripts b and x indicate that these will be domains of the corresponding variables. In 
order to define the vicinities of the Type I saddle points properly, we introduce the permutation of 
{1, 2}, for each triple Vj)- Specifically, recalling the fact of Uj = from (3.24), we define 

:= = ei) ^ ^i). 

Denoting by e = (ei,..., ew) and e(a) = (ei(a),..., evi/(a)) for a = 1, 2, we set 

^e{a) (^l,ei (a): ■ • ■ ; ^ ^ ^ (j32,e2 : ■ * ■ ; ). (5.36) 

With this notation, we now define the Type I, II, and III vicinities, parameterized by (bp, b 2 , Xp, X 2 , t, v) 
of the corresponding saddle point types. We also define the special case of the Type I vicinity, namely. 
Type r vicinity, corresponding to the Type T saddle point defined in Section 4.1. 

Definition 5.5. 

• Type I vicinity : (bp,b 2 ,x,(p),x,( 2 ),t,Ve) £ X X T- X T- X T 5 X T 5 for some e. 

• Type r vicinity : (bp, b 2 , Xp, X 2 , t, v) G x T5_ x x Tf x Tg x Tg. 

• Type II vicinity : (bp,b 2 ,xp,X 2 ,t, v) G T(j_ X x x x Tg x 

• Type III vicinity : (bp, b 2 , xp, X 2 , t, v) G T(j_ X x Tf x Tf x Tg x 1^“^. 
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In the following discussion, the parameter eg in 0 is allowed to be different from line to line. However, 
given £i in (1.7), we shall always choose £2 in (1.15) and £q in (5.33) according to the rule 

C£2 < ^0 < ^ (5.37) 

for some sufficiently large C > 0. Consequently, by Assumption 1.13 we have 

N{\0gN)-^° >M = > Vl/(4+27+ei)(l-4eo)jV/4eo ^ ^r{i+2J+4so )^ 1^27 04^ (5 gg) 

To prove Theorem 1.14, we split the task into three steps. The first step is to exclude the integral outside 
the vicinities. Specifically, we will show the following lemma. 

Lemma 5.6. Under Assumptions 1.1 and 1.13, we have, 

l(r'^,f^,E'^,E^,R5^-i,l'^-i) = 2 ^I(T';,T'L,T:;,T^,Ts,Ts) 

+ I(T^, T'L, T^, , Ts, I^-i) 

+ l(T^,T'L,T^,T^,Ts,I^-i) +0(e-®). (5.39) 

Remark 5.7. The first three terms on the r.h.s. of (5.39) correspond to the integrals over vicinities of 
the Type I, II, and III saddle points, respectively. Note that for the first term, we have used the argument 
in Section ).!, namely, the total contribution of the integral over the Type I vicinity is 2'^ times that 
over the Type V vicinity. 


The second step, is to estimate the integral over the Type I vicinity. We have the following lemma. 

Lemma 5.8. Under Assumptions 1.1 and 1.13, there exists some positive constant Cq uniform in n and 
some positive number Nq = Nq^u) such that for all N > Nq, 

AfC^o 

2^\l{TX,Tt,ri,T-_,Ts,Ts)\ < (5.40) 

The last step is to show that the integral over the Type II and III vicinities are also negligible. 


Lemma 5.9. Under Assumptions 1.1 and 1.13, there exists some positive constant c such that. 


= 0(e-=^), l(T^, T^, T^, , Ts,I^-i) = ©(e'^'^). 


Therefore, the remaining task is to prove Lemmas 5.1, 5.6, 5.8 and 5.9. For the convenience of the 
reader, we outline the organization of the subsequent part as follows. 

At first, the proofs of Lemmas 5.1 and 5.6 require a discussion on the bound of the integrand, especially 
on the term A(-), which contains the integral over all the Grassmann variables. To this end, we perform 
a crude analysis for the function A(-) in Section 6 in advance, with which we are able to prove Lemmas 
5.1 and 5.6 in Section 7. 

Then, we can restrict ourselves to the integral over the vicinities, i.e., prove Lemmas 5.8 and 5.9. It will 
be shown that in the vicinity (bi,b2,t) € x x T 5 , the factor exp{—ML(i3,T)} is approximately 
the product of a complex Gaussian measure of the H-variables and a real Gaussian measure of the t- 
variables. Here, by “complex Gaussian measure” we mean a function of the form exp{— u'Au}, where u 
is a real vector, while A is a complex matrix with positive-definite Hermitian part. In order to estimate 
the integral against this Gaussian measure (in an approximate sense), we shall get rid of the o(l) term 
in the integral of the form 

J du exp { — u'Au -I- 0 ( 1 )} f(u), (5-41) 

for some function f, which however cannot be done directly, owing to the fact that A is complex. In 
our case, this problem can be solved by further deforming the contours of the H-variables, following the 
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steepest descent paths exactly in the vicinity. By doing this, we can get a real Gaussian measure, thus 
the remainder terms can be easily controlled when integrate against this measure. The situation for 
exp{— V)} is a little bit more complicated due to different types of the saddle points. However, 
in the Type I vicinity, we can do the same thing. Hence, in Section 8 , we will analyze the Gaussian 
measure (in an approximate sense) exp{—M(/^(X, V) + L{B, T))}, especially, we will further deform the 
contours of X and H-variables in the vicinities, whereby we can prove Lemmas 5.8 in Section 9. In the 
Type H and HI vicinities, we will bound exp{—M/L(X, H)} by its absolute value directly. It turns out 
to be enough for our proof of Lemma 5.9, which is given in Section 10. 


6 . Crude bound on A{X,B,V,T) 

In this section, we provide a bound on the function A(-) in terms of the H, T-variables, which holds 
on all the domains under discussion in the sequel. Here, by crude bound we mean a bound of order 
exp{0(ITiV'^^)}), which will be specified in Lemma 6.1 below. By the definition in (3.31), we see that 
A(-) is an integral of the product of Q(-), 'Pi') and j'^(-)- We will mainly treat Q(-) as a function of 
variables, treat V{-) as a function of H, S- variables, and treat J^(-) as a function of y W,wW,Pi,Qi- 
variables. However, in the function Q(-), we actually have every argument mentioned above. Hence, 
we perform the integral over -variables and ^^^^-variables for Q(-) at first. The resulting function 
Q(-) turns out to be a polynomial of the remaining arguments. As mentioned in Section 4.2, a typical 
procedure we will adopt is to ignore Q(-) at first, then estimate the integrals of V{-) and P'(-)j which are 
denoted by P(-) and F(-), respectively (see (4.4) and (4.5)), finally, we make necessary comment on how 
to modify the bounding scheme to take Q(-) into account, whereby we can get the desired bound for A(-). 

For the sake of simplicity, from now on, we will use the notation 


‘^ 1,1 ■— >^j,3 ■— ■— ^j,22, 

0,1 ■= O.IU 0,2 := 0,211 0.3 := 0,12 > 0.4 := 0.22 ■ 


( 6 . 1 ) 


Moreover, we introduce the domains 


S := ■! re‘’^ : |r — 1| < 




K=K{E) := 


jw G C : 0 < argw < -f f|, if P > 0, 
{w G C : 2^2+ - I < argw < o}, if E < 0. 


( 6 . 2 ) 


By the assumption that \E\ < \/2 — k in (1.15), it is easy to see that | arga;| < ir/i — c for all w G KUK, 
where c is some positive constant depending on k. Our aim is to show the following lemma. 

Lemma 6.1. Suppose that bi,b 2 ,Xi,X 2 G x x x . Under the assumption of Theorem 
1 . 14 , we have 

w ^ 

IA(X,B,V,T)I < J] (^-1 +^-i + i) := 

J=i 

Remark 6.2. Obviously, using the terminology introduced at the end of Section 3.6, we have 

p(r"0t) G o({ri) 0 i'y 200 }i 0 i;Ki,K 2 ,K 3 ), ^ 2,^3 = 0(1). (6.3) 
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6.1. Integral of Q. In this section, we investigate the function 

Q(r!,S,Pi,Qi,xW,yW,wW) := J da;Wd|W S, Qi, Y™, w™). (6.4) 

Recall 0deg(a-; K 2 , ^ 3 ) defined at the end of Section 3.6, the parameterization in (3.15) and (3.24) and 
the notation introduced in (6.1). We shall show the following lemma. 

Lemma 6.3. If we regard a, v]p \ v]l\ Pi and -entries as fixed parameters, we have 

Q(-) e ndeg(6; K1,K2,R3), Ki = R2,R3=0(1), (6.5) 

where © is the set of variables defined by 

e:= .. 

^ ^ k—p Q'd 6=1 4 

Proof. Note that Q(-) can be regarded as a function of the Grassmann variables in and Hence, 
by the definition in (3.29), it is a polynomial of these variables with bounded degree. To have a closer 
look on this polynomial, we start with 

1 


exp{ - k=p,q. 


( 6 . 6 ) 


Observe that, in the exponent of ( 6 . 6 ), one ^[i'^^-variable must be accompanied by one fl-variable. In 

tb] 

uj o. J-iic-iij uy vviijii xoo^o'../U uc* 


addition, we combine the factor \I\[M with flj’s. Then, by Taylor expansion with respect to 
variables, it is easy to find 

exp|-^ G ndeg(©i,fc; ni, K 2 ,/ts), m = K2,K3 = 0{1), 


(6.7) 


where 


[ 1 ] ^i,aKk,h ' 


©i,fc.-{t,s,e , I 


y/M 


2 = 1 ,..., 4 : 6 = 1,2 


k=p,q. 


Analogously, we can show that for k = p,q 
r 1 


exp< - 


where 


y/M 


y^ SkjTru:^^\w^^^)*JQ^ G 0 deg( 62 .fc; ^2, ^ 3 ), ki = W°'^^\ K2,«^3 = 0(1), 

j 

( 6 . 8 ) 


P. 

— it s Sj.g k,b ) jL _ „ „ 

,fc.—-vC, s,e , > , fe—P, 9- 

I y/M J a=l,...,4;b=l,2 


& 2 ,k ■= s, e 

Here we have used the fact that Q^^-entries are the same as Qi'Sntries, up to a sign. In a similar manner, 
one can show that for k,i = p, q, 

1 


and 


exp{ - —Sfc«rra;|3'(w|3')*J(viff')^^^'j| G Odeg(© 3 ,fc,f; ^i, ^ 2 , K 3 ), ki,K 2 ,K 3 = 0(1), (6.9) 

(1 - G 0deg(©4,fc;Kl,K2,K3), Ki,K2,K3 = 0 (1), (6.10) 


where 


© 


W := {e'"!' 


W 




14 : 1 }.,. s,,:={toH)-', 4 :ui:i}, 

J a, 0 = 1,2 ’ ’ J a, 0 = 1 ,. 
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Hence, by (6.7)-(6.10) and (3.29), we see that 


Q(-) e C!deg(65;Kl,/t2,K3), Kl = K2,K,3 = 0(1), 


( 6 . 11 ) 


where 



^3,r^k,b [ 1 ] ^[ 1 ] ) 

Vm ’ Vm 


By the definition in (6.4), Q(-) is the integral of Q(-) over the and ^^^^-variables. Now, we regard all 
the other variables in G 5 , except and ^^^^-variables, as parameters. By the definition of Grassmann 
integral, we know that only the coefficient of the highest order term nfe=p g 00=1 2i 2 (') 
survives after integrating and -variables out. Then, it is easy to see (6.5) from (6.11), completing 

the proof. □ 


6.2. Integral of V. In this subsection, we temporarily ignore the O and S-variables from Q(-), and 
estimate P{X,B,V,T) defined in (4.4). Recalling and rj ^2 defined in (5.12), we can formulate our 
estimate as follows. 


Lemma 6.4. Suppose that the assumptions in Lemma 6.1 hold. We have 

w 

\P{X,B,V,T)\ < e°^^^ n (6-12) 

i=i 

Proof. We start with one factor from 'P{-) (see (3.28)), namely 

1 

det"" (1 + M-^v;x-W,n,T-^B-%E,) 

= exp { - Mlogdet (l -h Xr^V*Xj^V^VljT-^Bj^'E.j) | 

= 1 + T ^^ 3 ) ■ (6.13) 

i=\ ^ 

Here p£(-) is a polynomial in X~^, Sj-entries with bounded degree and bounded 

coefficients. Here we used the fact that V* and r^“^-entries are the same as Vj and Tj-entries, respectively, 
up to a sign. Moreover, if we regard p^(-) as a polynomial of VLj and S^-entries, it is homogeneous, with 
degree 2£, and the total degree for flj-variables is I, thus that for Sj-entries is also ^. More specifically, 
we can write 

4 t 

pi{Xj , Bj , Vj , Tj , Llj,^j) = Pe,a.i3(Xj , Bj , Vj , Tj ) ^j,Pi, 

ai,...,ae, i=l 

pi,...,Pi = l 

where we used the notation in (6.1) and denoted a = (ai,... ,a^) and /3 = (/3i,... ,/3|). It is easy to 
verify that vjj is of the form (6.13) by taking Taylor expansion with respect to the Grassmann variables. 
The expansion in (6.13) terminates at f = 4, owing to the fact that there are totally 8 Grassmann 
variables from Llj and Sj. In addition, it is also easy to check that pe,a,i 3 {') is a polynomial of X~^, B~^, 
Vj., Tj-entries with bounded degree and bounded coefficients, which implies that there exist two positive 
constants Ci and C 2 , such that 

\pi,a,i3{-)\ + r Ptj+ 1 ) ^ (6.14) 

uniformly in £, a and /3. Here we used the fact that X~^ and 1^-entries are all bounded and T^-entries 
are bounded by 1 -I- tj . 
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Now, we go back to the definition of P(-) in (3.28) and study the last factor. Similarly to the discussion 
above, it is easy to see that ioi k = p or g, 

det {V*XuVk + M-^nuT-^B-^Tu'^k) 

'^k • — 

det Bk 

4 I 

= -Bfc) + EE Bk^ Vfc, T]^') ^k,ak^i^k,j3k^i: (6.15) 

i—1 ol,I 3 i—1 

where po(’) = delX^/det^fc and p£,a,) 3 (*)’s some polynomials of Xu, B'^^, 14 , T^-entries with 
bounded degree and bounded coefficients. Similarly, we have 


|Po(-)li \Pl.,a,p{-)\ < + ’’fci + + 1)*^^ 

for some positive constants Ci and C 2 . 

According to the definitions in (6.13) and (6.15), we can rewrite (3.28) as 

w 

7’(11,E;, A:,i 3 ,y,r) = exp I - p P vjk- 


(6.16) 


(6.17) 




In light of the discussion above, nj!li '^j nfe=p q is a polynomial of A B ^,V,T,Vt and S-entries, 
in which each monomial is of the form 


1=1 k=p,q 
^-1 n -1 


W ij 


n n 

1=1 i=l 


(6.18) 


where we used the notation 


£ (.^ 1 ,..., ), oi ( 0 ^ 2 ,..., oi \^), f3 (/^i; • ■ •; ’ 

^1 “ ('^ 1 . 11 ■ ■ ■ (3j = {13jp, ■ .. ,/3j^i^), 

and cif 5 ^(’) is a polynomial of X,X~^, B~^, V and T-entries. Moreover, all the entries of i, a and (3 
are bounded by 4. By (6.14) and (6.16), we have 


w 


< eO(^) n +^,”2 +^4 + 1 ) 

1=1 


c 


(6.19) 


In addition, it is easy to see that the number of the summands of the form (6.18) in ri]=i Ilfe^p q 
is bounded by . 

Define the vectors 


^ ■ (^1 ; ^2; ^3; ^4); ^ ■ (^15^25^35^4); 

where 

iVfy = {uJi a: ... 5 U’W.ck); ~ (Ci,ck; • ■ ■ ; ^W.a); H ~ 1; ^5 3:. 

Here we used the notation (6.1). In addition, we introduce the matrix 

m = s®s®s®s. 


( 6 . 20 ) 


SjkTrVlj^ik = H1HI■=.^ 


It is easy to check 
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By using the Gaussian integral formula for the Grassmann variables (3.2), we see that for each <5 and 
/3, we have 

^ w 

/ dndE ■ exp ^-^SjkTrnjEk^ ■ ^ | detH^'l-'^l, (6.21) 

j,k j^l 

for some index sets I and J with |l| = |J|. By Assumption 1.1 (i) and (ii), we see that the 2-norm of each 
row of S is 0(1). Consequently, by using Hadamard’s inequality, we have 

|det]Hl('l-’)| = ( 6 . 22 ) 

Therefore, (6.17)-(6.22) and the bound for the total number of summands of the form (6.18) in 


U.Y=i '^3 Ilk=p,q '^k imply that 


w 


|P(A, B, fo, T)| < n +*3 + 1) 

3 = 1 

Thus we completed the proof. 


0 ( 1 ) 


□ 


6.3. Integral of B. In this subsection, we also temporarily ignore the Pi, Qi-variables 

from Q(-), and estimate F{X,B,V,T) defined in (4.5). We have the following lemma. 

Lemma 6.5. Suppose that the assumptions in Lemma 6.1 hold. We have 

k—p,q 

Proof. Recalling the decomposition of P(-) in (3.18) together with the parameterization in (3.27), we will 
study the integrals 


G{B,T) := J dKQi)dyWdwW g(gi, T, ^ y™, w™), 
F(1,R) := J d^l{Pl)dX+ /(Pi,R,l,XW) 


(6.23) 

(6.24) 


separately. Recalling the convention at the end of Section 3, we use /(•) and g{-) to represent the 
integrands above. One can refer to (3.19) and (3.20) for the dehnition. 

From the assumption 77 < it is easy to see 

\¥{X,V)\<e°^'^^"+, (6.25) 

since Pi, fo, A, Adl-variables are all bounded and | det a],^^|, | det Ag^^| ~ 1 when xi,X 2 G E defined in 

( 6 . 2 ). 

For G(P,r), we use the facts 

Re(rrPgyWj) > 0, Re{iTrY+JZ) =-riTrY+ <0, TrY^JY^JPO, k,£=p,q, 


|(wW(wW)*)^2l ^ 1’ I W’(wi^A)2il ^ 

to estimate trivially several terms, whereby we can get the bound 

w 

\g{-)\<exp{-MpJ2TrRe{B,)j} { - hkTrRe{Bk)Y+j}. 


(6.26) 


(6.27) 


3=1 


k—p,q 


Here Re{Bj) = ^T, '^Re{Bj)TjQi. Hence, integrating yY and yq^ out yields 




/R^ 


d2/Wcl2/W|g(Qi,T,p,yW,wW)|<C 


exp I - M 77 X;]IiT?’Re(Pj)j| 

Uk=p,q ((wW)VRe(P,)wW)^^ 


(6.28) 
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for some positive constants C and Ci depending on n, where we used the elementary facts that Skk > c 
for some positive constant c and 

TrRe(B,)yW J = yW(wW)VRe(B,)wW, k=p,q, j = (6.29) 

Now, note that 

JRei?jw[,^' > Xi{JReBj), k=p,q, j = I,... ,W. (6.30) 

In addition, it is also easy to see Xi{Tj) = Sj — tj and Ai((5i) = s — t, according to the definitions in 
(3.24). Now, by the fact JA~^ = AJ for any A € U{1, 1), we have 

JReBj = QiTjdia.g(^Rebj^i,Rebj^ 2 )TjQi. (6.31) 


Consequently, we can get 


Ai(JRei?j) > {sj — tjY{s — t)^min{Re&j 4 , Re 6 y 2 } = 


min{Re 6 yi, R&bj 2 \ 

(Sj +t/) 2 (s + t )2 


(6.32) 


by recalling the facts = 1 and s^ — t^ = 1- Therefore, combining (6.28), (6.30) and (6.32), we have 


w 

d4^'dyW|5(-)l < C'(s + t)'‘^iexp|-Mr7^rr(ReBj)j| 


(sfe + 


i=i 


k^p g (min{Re 6 fe,i, Re 6 fc, 2 }) 


Cl ' 


(6.33) 


Now, what remains is to estimate the exponential function in (6.33). By elementary calculation from 
(6.31) we obtain 


Tr(ReBj)T ^ (Refo^^i 4- Rel>j 2)(('Sj -t-tj)(s -t-t ) — dstsjtj'^. 


Observe that 


(s^ + + t^) — dstSjtj = 


(sj + ^ + 4t^ + 4t‘f + 4t^ + 1 ^ 1 + 2t 


(sj + + f^) + 4stsjtj 


> 


3 ' “j 


> 


2 (l + 2t2)(l + 2t2) -2(l + 2t2)- 


It implies that 


i=i 


exp { - Mi^Y.'^r{ReB,)j] < exp { - 2 M, 7 ^ + 2 ^^)} 

1=1 -1 

< exp { - cMiq ''W 2 ^^ (1 + 2t^)}, 


(6.34) 


1=1 ' -1 

for some positive constant c, where in the last step we used the fact 


R^bj^a > cTj^a, V j = 1,..., IT, a = 1,2, (6.35) 

for some positive constant c, in light of the assumption |iil| < \f2 — k and the definition of K in ( 6 . 2 ). 
Plugging (6.34) into (6.33), estimating (s + t)^ < 2{l + 2t‘^), and integrating t out, we can crudely bound 




1 v,.!+r,-. 2 \-C '2 


, rr_ 

^ l + 2 t 2 y ((min{Re&fc.i, Re&fe, 2 })‘^' 


(1 + 2 t 2 )Ci 


Now, we use the trivial bounds 

w , ^ 

I’i.i + r,,2 \ -C 2 




and 


l + 2tl 


min{Re5fc,i, Re&fe,2} 


< C(1 + 2t/.)(rf.. + ,). 


. (6.36) 


(6.37) 


(6.38) 
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Inserting (6.37) and (6.38) into (6.36) and integrating out the remaining variables yields 

\G{B.T)\<c(^y^ n + (6-39) 

k=p,q 

Combining (6.25) and (6.39) we can get the bound 

k—p,q 

Hence, we completed the proof of Lemma 6.5. □ 

6.4. Summing up: Proof of Lemma 6.1. In the discussions in Sections 6.2 and 6.3, we ignored the 
irrelevant factor Q(-). However, it is easy to modify the discussion slightly to take this factor into account, 
whereby we can prove Lemma 6.1. 

Proof of Lemma 6.1. At first, by the definition in (6.4), we can rewrite (3.31) as 

A(-) = J dAWdyWdwWdHdSd^(Pi)dj^(Qi) P(-)Q(-)-^(-)- 

Now, by the conclusion ki = in Lemma 6.3, it suffices to consider one term in Q(-), which is a 

monomial of the form 

s, (j/W)-\ (yW)-i)q(H,S), 

regarding cr, rip v^q \ Pi-variables, A [^1-variables and wl^l-variables as bounded parameters. Here p(-) is a 
monomial of t, s, (?/p^^)“^, and q(-) is a monomial of H, S-variables, both with bounded coefficients 

and bounded degrees, according to the fact K 2 ,K 3 = 0(1) in Lemma 6.3. Now we define 

Pq{X,B,V,T) := J dndEVin,E,X,B,V,T)-qin,E), 

Fp(A,H,t/,T) := J dAWdyWdwWd/4(Pi)dz/(gi) 

xP(A, B, V, T, Pi, Qi, aW, y W, w™) • p{t, s, (y™)"!, (y™)"!). (6.40) 

By repeating the discussions in Sections 6.2 and 6.3 with slight modification, we can easily see that 

w 

|Pp(A, B, V,T)\< + tj + 

4=1 

|Fp(A,P,H,T)|<eO(^^'^) n + 

k=p,q 

hold as well. Therefore, we completed the proof of Lemma 6.1. □ 

7. Proofs of Lemmas 5.1 and 5.6 

In this section, with the aid of Lemma 6.1, we prove Lemmas 5.1 and 5.6. According to Lemmas 
5.2 and 5.3, one can see that away from the saddles, ReL{B,T) and ReA(A,H) increase quadratically 
in P-variables and A-variables, respectively. Hence, it is easy to control the integral (5.2) over these 
variables outside the vicinities. However, like P-variables, the domain of t-variables is also not compact. 
This forces us to analyze the exponential function 

M(t) :=exp{-MRe(£s(-B,T))| (7.1) 

carefully for any fixed P-variables. 
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Recall the definition of the sector K in (6.2). For bi G and b 2 G we have 


minRe( 6 j 4 + + &fc. 2 ) > cinin 

j,k i,k 


E 


rj,ark,b > cininr^^ 
j,a 


:= 2t(B), 


(7.2) 


for some positive constant c depending on k from (1.15). From now on, we regard M(t) as a measure of 
the t-variables and study it in the following two regions separately: 


(f) : t G (ii) : t G \ 

Roughly speaking, when t G 1^“^, we will see that M(t) can be bounded pointwisely by a Gaussian 
measure. More specifically, we have the following lemma. 


Lemma 7.1. With the notation above, we have 

M(t) < exp I - ^^{B)'^Sjk{tk - tjf\, Vt G 


However, the behavior of M(t) for t G \ is much more sophisticated. We will not try to 

provide a pointwise control of M(t) in this region. Instead, we will bound the integral of q(t) against 
M(t) over this region, for any given monomial q(-) of interest. More specifically, recalling the definition 
of 0 in (5.33) and the spanning tree Go = (V,£’o) in Assumption 1 . 1 , and additionally setting 

M 

^ ■= l“2l(-B) min s^-, (7.3) 

we have the following lemma. 


Lemma 7.2. Let q(t) = 0^=2 ® monomial of t-variables, with powers nj = 0(1) for all j = 

2,..., W. We have 


r JfL / ,sO(W=) r ^ ^ 

/ TTdtj M(t)q(t) < (l + £-5 exp<^-022l(R)+O(W2logAf)^ 

V y I j 


(7.4) 


Remark 7.3. Roughly speaking, by Lemma 1.2 we see that the integral of q{t)-variables against the 
measure M(t) over the region \ is exponentially small, owing to the fact 0^ ^ W'^ logiV. 


We will postpone the proofs of Lemmas 7.1 and 7.2 to the end of this section. In the sequel, at first, 
we prove Lemmas 5.1 and 5.6 with the aid of Lemmas 6.1, 7.1 and 7.2. Before commencing the formal 
proofs, we mention two basic facts which are formulated as the following lemma. 


Lemma 7.4. Under Assumption 1.1, we have the following two facts. 

• For the smallest eigenvalue of there exists some positive constant c such that 


(7,5) 


• Let Q = {p 2 , ■.., pw)' be a real vector and pi = 0. If there is at least one a G {2, ... ,W} such 
that Qa > then we have 


Proof. Let g 
Then 


^jk{Pj Qk) — ("^' 6 ) 

3,k 

{p 2 ,. ■., pwY be a real vector and pi = 0. Now, we assume \pa\ = max^= 2 ....,w \Pfi\- 
-p'S^^'lp ^Y.j,k^lk{Pj - Pkf ^ c{pa- Pl)^ C 

\\p\\l ~ - w^pi 
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where the second step follows from Assumption 1.1 (iv) and Cauchy-Schwarz inequality. Analogously, we 
have 

according to the definition of 0 in (5.33). Hence, we completed the proof. □ 


Recalling the notation defined in (5.2) and the facts \xj^a\ = 1 and \bj^a\ = »"j.o for all j = 1,..., VF 
and a = 1 , 2 , for any sequence of domains, we have 

|i(i?,i^,iT,i^,i*,r)| 

^ w ^ w ^ w ^ w ^ w ^ w 

< gO(Wlog7V) 


I Y[^h 2 n 2vjdvj 


1 


2 j = l 


*1 j = l 


2 i=i 


j=2 


3=2 


W 

X exp{-M(ReK(A,H) + Rei(H,T))} • |A(A, H, H, r)| • . (7.7) 

1=1 

In addition, according to Lemma 6.1, we have 

WWW 

\A{X,B,V,T)\ l[2t, ■l[2v, ■ n(^4.i +^.. 2 )" < eO(^^'^)p(r,r-i,t), (7.8) 

1=2 1=2 1=1 

for some polynomial p(r,r“^,t) with positive coefficients, and 

p{r,r~^,t) & Q({rj,i,rj,2,rjl,rj2,tj}^^-Ki,K2,K3'^, ki = K 2 ,K 3 = 0(1). (7.9) 

7.1. Proof of Lemma 5.1. At first, since throughout the whole proof, the domains of Xi, X 2 , and 
v-variables, namely, E^, E'^ and will not be involved, we just use *’s to represent them, in order 

to simplify the notation. 

Now, we introduce the following contours with the parameter D G K+, 
r© := {ra+|r G [0,2)]} C L, R© = [0, (Rea+)2)] C R+, £© := {(Rea+)2) + i(lma+)r-|r G [0,2)]}. 

In addition, we recall the sector K defined in (6.2). Then, trivially, we have 


M-i-, r, /l© c IK, IR_|_, r, >C© € IK, V 2) € . 

We claim that the integrand in (5.2) is an analytic function of the H-variables. To see this, we can 
go back to the integral representation (3.17) and the definitions of L{B) and X, B) in (3.18). 

Note that since exp{Mlogdeti?i} = (detSi)'^, actually the logarithmic terms in L[B) do not produce 
any singularity in the integrand in (3.17). In addition, according to the fact that the = 0 for any 
Grassmann variable y and f > 2, the factors det^(I + M“^A“^Hji?“^Sj) is actually a polynomial of 
Sj, and i?~^-entries with degree 16. The other factors containing R-variables V{-) can be checked 
analogously. Hence, it is easy to see that exp{—ML(-)}P(-) is analytic in H-variables. Consequently, we 
have 

i((r© u £©)’^, (f© u £©)’^, *, = i((R©)^, (R©)’^, *,*,*). 

Hence, to prove Lemma 5.1, it suffices to prove the following lemma. 


Lemma 7.5. Suppose that \E\ < \f2 — k. As 2) —>■ 00 , the following convergence hold, 

(z) : l((r©U£©)^,(f©U£©)^,*,*,R)^'-\*) - i((r©)^, (f ©)'^, ^ 0, 

{it) : i((r)^,(f)^,*,*,R5!'-\*) -i((r©)^,(f©)^,*,*,R)^-\*) ^0, 
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-X 


(m) : 

Proof. For simplicity, we use the notation 






0 . 


1^1 := (Fs, X {T^UC^r \ X (fs,)^, 

rb,2 


t&,3 


;; 

\VE 


= T^ xr'^\{r^r X (Fs,)'^, 

= X Rf \R^ X Kg'. 


Now, recall the definition of the function £(a) in (5.6) and the representation of L{B,T) in (5.5). Hence, 
in light of the definition of M(t) in (7.1), we have 

exp { — MReL{B, T)} = exp { — M(Ref(bi) + Re€(—b 2 )) }M(t). (^-lO) 

By the assumption |F1| < — k, we see that Rebj^abk,a > 0 for all bj^a, bk,a £ IKUlK. Consequently, when 

bj^i £ K and bj ^2 £ K for all j = 1,..., W, we have for some positive constant c dependent of k in (1.15), 

Re^(bi) + Re€(-b2) > X! X! (^(^ + + (-l)“"^^-E'lm6j,Q - logr^.a) 

a=l,2 j 

^ ^ J2 J2da - J2 (7.11) 

a— 1,2 j a— 1,2 j 

where we used Assumption 1.1 (ii) and the fact that (—l)°'~^^E\mbj^a > 0. 

Now, when (bi ,b 2 ) £ lg‘ for i = 1 , 2 ,3, we have X]a=i 2 some positive constant c, 

which implies the trivial fact 

EEd.a)EEd. + i®*- E-12) 

a— 1,2 j 0 — 1,2 j 

Consequently, we can get from (7.10), (7.11) and (7.12) that for some positive constant c, 

w 

exp { - MReA(.B, T)} < Y[ a . M(t) 

0=1,2 j = l 

holds in ig* for i = 1, 2, 3. In addition, by the boundedness of V and A-variables, we can get the trivial 
bound MK{X,V) = 0{N). Hence, from (7.7) and (7.8) we see that the quantities in Lemma 7.5 (i), (ii) 
and (iii) can be bounded by the following integral with i = 1,2,3, respectively. 


w 


w 


w 


w 


w 


^0(W log N) 


w 

X exp ^-M(ReK(X,V) + ReL(B,T))j ■ |A(A,H,H,T)| • +C, 2 )" 


2vjdvj 


1=1 


< e 


0{N) -cMD'^ 


VE 

I® 1=1 


EE 


EE 


n "i>-M(t).p(r,r-gt). (7.13) 

1=2 0 = 1 , 21=1 

According to the facts ki = and K 2 = 0(1) in (7.9), it is suffices to consider one monomial in 

p(r,r“^,t) with bounded coefficient. That means, it suffices to estimate the integral 

» EE „ EE 

/. . d^l n * = 1,2,3, (7.14) 




0 = 1,2 1=1 


1=1 

for some monomial 

EE EE 

q(i’,i’"gt) = |£i,o|,ni = 0(1), for j = l,...,IF, a = l,2, (7.15) 

0 = 1 , 21=1 1=2 
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where the bound on ij's and n^’s follows from the fact that K 3 = 0(1) in (7.9). 

Bounding t^-’s by 1 trivially in the region t G 1 ^ 7-1 and using Lemma 7.2 in the region t G 
we can get for * = 1,2, 3, 


„ w w 

J \ ^ A — 1 — 10 — 1 


i = l 


a=l,2 


By the definition of 2l(i?) in (7.2) and the assumption M ^ we see that 


w 




w 


M(l+o(l)) 

j,a 


a=l,2 j=l a=l,2 j=l 

Consequently, by using elementary Gaussian integral, we can get the trivial bound 

(7.14) < i = 1,2,3, 

and then we have for i = 1, 2, 3, 

(7.13) < ^0, as 00 . 

Thus we completed the proof. 


□ 


7.2. Proof of Lemma 5.6. Plugging the first identity of (5.20) and (7.8) into (7.7), we can write 


w 


w 


w 


w 






Tfe 

b j=l 


jb 

b i=l 


^1 j=l 


^2 j=l 


W „ W 


X I dtj duj exp |—M ^ReiL(X, F) + ReL(.B, T)^ I p(r, r ^,t), (7.16) 

•_Q J 


i=2 "^ j=2 

where p(r,r”^,t) is specified in (7.9). 

Lemma 5.6 immediately follows from the following two lemmas. 


Lemma 7.6. Under Assumptions 1.1 and 1.13, we have 

l(r^,f’^,E'^,E^,R)^-i < e-^\ (7.17) 

Lemma 7.7. Under Assumptions 1.1 and 1.13, we have 

l(r^,f^,E^,E'^,l'^-i,I^-i) = 2'^l(T(;,Ti,T“,T“,Ts,Ts) 

+ I(T';, T'L, , Ts, I^-i) 

+ l(T^,T'L,T“,T^,Ts,l‘^-i) +0(e-®). (7.18) 


In the sequel, we prove Lemmas 7.6 and 7.7. 


Proof of Lemma 7.6. Recall (7.16) with the choice of the integration domains 

To simplify the integral on the r.h.s. of (7.16), we use the fact ReK{X, F) > 0 implied by (5.26), together 
with the facts that the x and v-variables are bounded by 1. Consequently, we can eliminate the integral 
over X and v-variables from the integral on the r.h.s. of (7.16). Moreover, according to (7.9), it suffices 
to prove 


w ^ w . w 

TT d6j 1 / TT d6j 2 / TT dtj exp MReL(i?, T)| • q(r, r“^, t) < e“® (7.19) 

7=1 ’ •^r'^7=i ’ JJ 2 '' ^ 

instead, where q(-) is the monomial defined in (7.15). 
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Now, by the first inequality of (5.13), we have 

w 

exp( - MReL(B,T)} < exp^-M ^ ^ - 1)^ + - loga - 1)) | • M(t). (7.20) 

a=l,2 j=l 

At first, we integrate t-variables out by using Lemma 7.2, namely, 

w 

M(t) • q(r,r"\t) 


gW-l\jvi/-l “ 


w 


< 


n + ^ exp^-0^2l(B)+0(W^ log JV)j. 

a—1,2 j—1 

Recall the definitions of £ and Ql(B) in (7.3) and (7.2), and the fact = o(M^/^). The proof of the 
following fact is an elementary exercise 


(7.21) 


sup exp {-M(r-logr- 1 )}/= 0 ( 1 ), if i=0{Mi). 

rGR+ 


(7.22) 


Hence, we get the bound 
w 


n n ^ “ 1)} ■ ■ (i+") 




a=l,2 j=l 


VK 


< max 

6=1.2 /c=l 


1 + ^k,b) 


“ n n - log -1)} ■ ( 

' ’ ' a=l,2j = l 

Consequently, (7.20)-(7.23) imply that 

w 

TT dtj exp MRei/(.B, T)| • q(r, r~^, t) 
-1\ -^4; I J 


o(rv=) 


= ). 


(7.23) 




< 


w 

exp {-0221(H)+ 0(1^2 log iV)}- Y[ exp {-cM(r,-„ - l)^} 


a=l,2 j = l 


W 


a=l,2 j=l 

for some positive constant c, where in the last step we use the obvious fact 

w 

022t(H) + > 0" » log N 

^ a=l,2 j=l 


(7.24) 


by (5.38) and the definition of 2t(H) in (7.2). Plugging the bound (7.24) into the l.h.s of (7.19) and taking 
the integral over H-variables we can see that (7.19) holds, which further implies (7.17). Therefore, we 
completed the proof of Lemma 7.6. □ 

To prove Lemma 7.7, we split the exponential function into two parts. We use one part to control the 
integral, and the other will be estimated by its magnitude. More specifically, we shall prove the following 
two lemmas. 


Lemma 7.8. Under Assumptions 1.1 and 1.13, we have 

r ^ r ^ r ^ r ^ r ^ 

L, L, 


/rw 


i=i 








Ijw-l 


i=2 




X exp { - ^M{ReK{X, V) + ReL(H, T))} • p(r, t) < . 


(7.25) 
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Lemma 7.9. //(bi, b 2 , xi, X 2 , t, v) £ x x x x ^ x ^, but not in any of the Types 

I, II, III vicinities in Definition 5.5, we have 

exp I - ^M{ReK{X,V) + ReL{B,T))'^ < e"®. (7.26) 

With Lemmas 7.8 and 7.9, we can prove Lemma 7.7. 

Proof of Lemma 7.7. For the sake of simplicity, in this proof, we temporarily use Ifuii to represent the 
l.h.s. of (7.18), i.e. the integral over the full domain, and use X/, Xu and I/// to represent the first three 
terms on the r.h.s. of (7.18). Now, combining (7.16), (7.26) and (7.25), we see that, 

- X/ - Xu - XuA < ■ e-® • < e"®, 

in light of the definition of 0 in (5.33) and the assumption (5.37). Hence, we completed the proof of 
Lemma 7.7. □ 


Proof of Lemma 7.8. At first, again, the polynomial p(-) in the integrand can be replaced by the monomial 
q(-) defined in (7.15) in the discussion below, owing to the fact that ki = exp{0(lF)} in (7.9). Then, 
the proof is similar to that of Lemma 7.6, but much simpler, since t-variables are bounded by 1 now. 
Consequently, we can eliminate X, t and v-variables from the integral directly and use the trivial bounds 

w w 

q(r,r-\t)< ReX(.B,r) > c ^ ^(r^-« - 1)^ (7.27) 

a—l,2j=l a—l,2j—l 


where the latter is from (5.13). Hence, it suffices to show 

WWW 

f Ild&j,! / f[dbj,2 Y[ f[exp{-cM{rj^a-lA}rfA (7.28) 

a=i,2i=i 

Note that (7.28) follows from elementary Gaussian integral immediately. Therefore, we completed the 
proof of Lemma 7.8. □ 


Proof of Lemma 7.9. At first, according to (5.13) and (5.26), we see both MReL{B, T) and MReK{X, V) 
are nonnegative on the full domain. Hence, it suffices to show one of them is larger than 0 outside the 
Type 1, 11, HI vicinities. 

Note that for each type of vicinity, we have 

(bi,b 2 ,t) G X T^_ X Ts. (7.29) 

Now, if (7.29) is violated, we have (bi, b 2 ) G F^ x F^ \ T(j_ x or t G \ Tg. If the former holds, 
by using (5.13) and the definition of and in (5.35), we have 

w 

MReL{B,T) > cM ^ - if = cM\\hi - a +||2 + cA/||b 2 + a _||2 > 0, 

a=l,2 j=l 

which shows (7.26) if (bi,b 2 ) G F'^ x \ T*; x T^. 

Hence, it suffices to consider the case (bi,b 2 ) G T(j_ x T'L, t G 1'^“^ \ Tg. Using Lemma 7.1, we can 
see that 

M(t) < exp | —^2t(H)(—< exp 2l(i?)0| < e“®, (7.30) 

where in the second step we used the definition of Tg in (5.35) and in the last step we used the fact 
21(5) > c if (bi,b 2 ) G T'), X T^. Then (7.20) and (7.30) also imply (7.26). Now, we turn to show 
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MReK{X^ F) > 0 outside the vicinities. Recalling the definition of 'dj^s in (5.25), we split the discussion 
into two cases 

(i) : (siniJj - < ^, V j = 1,..., 21R, (m) : (^siniJ^ - > ^, for some j £ {1,..., 2VF}. 

Using (5.26), we can get MReK{X, U) > 0 in case (ii) immediately, so we can assume case (i) below. 
Then (i) implies that 

|arg(a+^a;ya)P A |arg(aI^Xj,a)p < ^, V j = 1,..., lU; a = 1, 2. (7.31) 

Now, we claim that it suffices to focus on the following three subcases of (7.31), 

(i’) There exists a sequence of permutations of {1, 2}, namely, e = (ei,..., ew), such that 


1 a.Tg{a_^_^Xj, 


0 

< — , 
- M’ 

1 arg(a; 

- 

,0-(2))P 

0 

- M’ 

II 

> 

(ii’) There exists 








1 

arg(a“ 


0 

< — , 

- M’ 

V j 

= 1 ,.., 

.,W;a = 

= 1,2. 

(iii’) There exists 








1 

arg(al 


VI 

V j 

= 1 ,.., 

.,lU;a = 

= 1,2. 


To see this, note that for those N-variables which satisfy (7.31) but do not belong to any of the case 
(i’), (ii’) or (iii’) listed above, there must be a pair {i, j} £ £ such that 


1 arg(a+^a;,,i)p, | arg(a+^a;,,2)P < 

1 arg(a_^a;yi)p, | arg(a_^a;y2)P < 

(7.32) 

or there exists a permutation such that 



1 arg(a+^a;i,,,(l))|^ | arg(al^a;i,,,(2 ))|^ < 

1 arg(a:f^a::j-1)|^ | arg(a:f^xy2)|^ < 

(7.33) 

or 

|arg(a+^a;i_,.(i))|Marg(aI^a;,_,,(2))P < 

|arg(aI^a::j-i)|Marg(aI^a:y2)P < 

(7.34) 


For each of (7.32), (7.33) and (7.34), we can perform a discussion similar to (5.30)-(5.31), to show that 
MReK{X, V) > cM » 0. 

Hence, it suffices to focus on cases (i’), (ii’) and (iii’) in the sequel. Now, we denote the domains of 
N-variables satisfying (i’), (ii’) and (iii’) by TJ, and respectively. In addition, in the remaining 

part of this proof, we temporarily use the term Type A vicinity to represent its restriction on A-variables, 
for A = I, II, III. Obviously, we have 

Type A vicinity C T’^, A = I, II, III. 

If (xi,X 2 ) £ Tf but outside the Type I vicinity, 

II arg(a:^^Xe(i ))||2 + || arg(aI^Xe( 2))||2 > ^, V e, 

which easily implies that 

2W _ 2 

MReA:(A, V)>cMJ2{ sin > ©■ (7-35) 

i=i 

Now, we turn to the case that there exists one sequence of permutations e such that (x^ji),Xg( 2 )) G 
X T^, but Vg ^ Tg. In this case, we just go back to the first line of (5.15) and do the transform 

Vj —^ t jXjJ 
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for those j with e{j) ^ e(l), where 3 is defined in (4.1). Then, it suffices to consider 

(xi,X 2 ) G X T!i, but V ^ T 5 , or (xi,X 2 ) G Tf. x but v ^ T 5 . 

In either case, we can show that MReK{X,V) > 0, analogously to case of (bi,b 2 ) G x but 

t ^ Ts, in (7.30). 

Now, what remains is to show that for those (xi,X 2 ) G but outside the Type A vicinity (A=II, 
III), we have MReK{X,V) > 0. We only discuss the case A = II, the other is analogous. Note that 
outside the Type II vicinity of X variables we have 

||arg(a;^xi)||^ + 11 arg(a:f^X2)|I 2 > (7.36) 

Observe that now we are already in T|j, which means that all Xj^aS are close to a+ and far away from 
a_. That means, we have sin(arg(a;ja)) — E/2 ~ arg(al|l^a;ja). Consequently, (7.36) also implies (7.35). 
Therefore, we completed the proof of Lemma 7.9. □ 

7.3. Proof of Lemma 7.1. Using the definition of is{B,T) in (5.6) and 2l(i3) in (7.2) and the fact 
|(rfcr“^)i 2 | = we have 

Re£s{B,T) > (7.37) 

j,k 

Simple estimate using = 1 + shows that 

\sjtke — SktjB ^1 > —{tk — tj) (2 _|_ 2^2 _|_ 2 ^ 2 ) — ~ ^ 3 ) ■ (7.38) 

Notice that the assumption t G 1'^“^ was used only in the last inequality. By (7.37), (7.38) and the 
definition (7.1), Lemma 7.1 follows immediately. 


7.4. Proof of Lemma 7.2. Let F = R+ \ I. Now we consider the domain sequence J = (J 2 , ... ,Iw) G 
{1,1'^}'^“^. We decompose the integral in Lemma 7.2 as follows 


w 


w 


/RW- 


]^dtj M(t)q(t) = 






dtj M(t)q(t). 

' i=2 


(7.39) 


Note the total number of the choices of such J in the sum above is — 1. It suffices to consider one 
of these sequences J G {I, in which there is at least one i such that Ji = F. 

Recall the spanning tree Qg = (V;^o) in Assumption 1.1. The simplest case is that there exists a linear 
spanning tree (a path) Qg with 


£g = {(^,^ + l)}Z7' C£. 


(7.40) 


We first present the proof in this simplest case. 

Now, we only keep the edges in the path £g, i.e. the terms with fc = j — 1 in (7.37), we also trivially 
discard the term 1/(1 + 2tj) from the sum 1/(1 + ‘2-tZi) + 1/(1 + in the estimate (7.38) (the first 
inequality), and finally we bound all M2l(ii)Sj_ij/4 by £ defined in (7.3) from below. That means, we 
use the bound 

M(t) < [] exp { - S"’ } := n 

3=2 j=2 


(7.41) 
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Consequently, we have 


w 


w 


w 


^03 = 2 j=2 


dtj M[(t)q(t) < 




nj =2 ^3 j=2 j=2 


(7.42) 


Note that, as a function of t, Mj(t) only depends on tj-i and tj. 

Having fixed J, assume that k is the largest index such that Ik = i-e. tfc+i, ■ ■ ■ ,tw S I. Now, we 
claim that 


w 

E 




^E 


3=2 l + 3=2 1 + 

To see (7.43), we use the following elementary facts 


> 


300fc2’ 


if tk€F. 


{tj - tj-if ^ 1 {tj - tj-if 


l + 2t]_. 


> 


ti 


i-i 


= - l)^ if tj-i G H"" 


and 


{tj tj-l)'^ > ^ ('f _/ '(2 -r . 

1-|_ 2^2 — C-i) J C-i Gi 

for all j = 2,..., W. We show (7.43) by contradiction. If (7.43) is violated, we have 

{t,-t,.,r ^ 

^-n j V J 


(7.43) 


(7.44) 


(7.45) 


1 + 


300fc2^ 


which together with (7.44) and (7.45) implies that 

+ + lo¬ 

using (7.46) recursively yields 


/ 1 \fc-i / i \K-i 1 


1 \ ^-1 


(7.46) 


(7.47) 


where in the second step we used the fact ti = 0. Note that (7.47) contradicts tk G F. Hence, we verified 
(7.43). 

Now, we split n ^2 (f) i'2fo parts. We use one to control the integral, and the other will be 
estimated by (7.43). Specifically, substituting (7.43) into (7.42) we have 
» w „ w w 


>n 


dtj M(t)q(t) < 


e 6 oofe^ 


j=2 d'i j—2 




*+ i =2 


(7.48) 


J=2 


Therefore, what remains is to estimate the integral in (7.48), which can be done by elementary Gaussian 
integral step by step. More specifically, using (7.44) and (7.45) and the change of variable tj/tj-i — l^tj 
in case of t 3 _i € F and tj — tj-i —>■ tj in case of tj-i G I, it is elementary to see that for any £ = 0{W)^ 

f dtj t^^{Mj{t))^ <£\\[l + ct~^\^^^'‘ (tjtl + l) < (7.49) 

J K_|_ 

Starting from j = IT, using (7.49) to integrate (7.48) successively, the exponent of tj increases linearly 
{rij = 0(1)), thus we can get 

/ J^dtj M(t)q(t) < e"5o^ + 

■^03=2 b 3=2 

Then (7.4) follows from the definition of £ in (7.3) and (5.38). Hence, we completed the proof for (7.4) 
when the spanning tree is given by (7.40). 
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Now, we consider more general spanning tree Go and regard 1 as its root. We start from the general¬ 
ization of (7.41), namely, 

M(t)< Y[ exp{-£(^^|f} := 

{*,i}6fo * {i,j}e£o 

Here we make the convention that dist(l,i) = dist(l, j) — 1 for all {i,j} G £o, where dist(a, b) represents 
the distance between a and b. Now, if there is k' such that Ik' € we can prove the following analogue 
of (7.43), namely, 

i'tj — uY ^ 1 

^ 1 -k 2t2 - 300fc2 

{*.j} 6 £o * 

by performing the argument in (7.44)-(7.47) on the path connecting k' and the root 1. Consequently, we 
can get the analogue of (7.48) via replacing Mj(t)’s by Mij (t)’s. Finally, integrating tj's out successively, 
from the leaves to the root 1, yields the same conclusion, i.e. (7.4), for general Go- Therefore, we completed 
the proof of Lemma 7.2. 


8. Gaussian measure in the vicinities 

From now on, we can restrict ourselves to the Type I, II and III vicinities. As a preparation of the 
proofs of Lemmas 5.8 and 5.9, we will show in this section that the exponential function 

exp{-M(A:(X,H) + L(H,T))} (8.1) 

is approximately a Gaussian measure (unnormalized). 

8.1. Parametrization and initial approximation in the vicinities. We change the x, b, t, v- 

variables to a new set of variables, namely, x, b, t and v. The precise definition of x differs in the 
different vicinities. To distinguish the parameterization, we set >£■ = ±, -I-, or —, corresponding to Type 
I, II or III vicinity, respectively. Recalling from (1.24). For each j and each >f, we then set 

Xj = H^diag ^exp{ixj_i/V^},exp{ii;j_2/V^}^ , Xj^a/VM G [-7r,7r], 

Bj = D± + D±diag(hj^i/y/M,bj^2/^/^'j , tj (8-2) 

If x: = ±, we also need to parameterize Vj by 

Vj=Vj/'/M- (8.3) 

We set the vectors 

ba ■ — (^l,a5---5 b\^^a) ; • ■ • ; Il'IV.a ); R — 1; 2, 

t := (I 2 , ■. -Jw), V := (v 2 ,.. ■,vw)- 

Accordingly, recalling the quantity 0 from (5.33), we introduce the domains 

t = t(N, eo) := {a e : I |a| 1^ < 0}, tg = ts(iV, Sq) := {a S : -a'^S^^^a < 0}. 

We remind here, as mentioned above, in the sequel, the small constant Sq in T and T g may be different 
from line to line, subject to (5.37). Now, by the definition of the Type F, II and III vicinities in Definition 
5.5 and the parametrization in (8.2) and (8.3), we can redefine the vicinities as follows. 

Definition 8.1. We can redefine three types of vicinities as follows. 

• Type r vicinity : (bi, b 2 , Xi, X 2 , t, v)GTxTxTxTxTgxTg, with x = ±. 

• Type II vicinity : (bi,b 2 ,Xi,X 2 ,t, v) sTxTxTxTx Tg x with }<= +. 

• Type III vicinity : (bi, 1)2, Xi, X 2 , t, v) sTxTxTxTx Tg x with h = 
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We recall from (7.6) the fact 


teTs^||t|U = O( 0 ). 


(8.4) 


Now, we use the representation (5.2). Then, for the Type I vicinity, we change x, b,t, v-variables to 
X, b, t, v-variables according to (8.2) with x: = ±, thus 

w !. w 

2^l(T(;,T^,T“,T(i,Ts,Ts) = 




(n!)^4^7r2^+4 J^2w-2 27r 27r Jf 

\ I -^2 3=2 




^ w w r ^ r ^ °4-° 

. J I l[ 2 v,dv, f[exp{i^^^^} 

I I J Is ,'=2 Is 1=2 1 = 1 \ IVI 


i=i 


i=i 


i =2 


j=2 


4 = 1 


w 

X exp { - M{KiX, y) + L{B, T ))} • [](x,,i - x,, 2 )"( 6 ,.i + 6 ,, 2 )' • A(X, B, F, T). (8.5) 

4=1 

For the Type II or III vicinities, i.e. x: = -|- or —, we change x, b, t-variables to x, b, t-variables. 
Consequently, we have 

w . w 


f—n2 W 

*'& 'yb yx yx y ttW— 1 ^ _ \ ^>c) 


l(T^,Ti,T^,T^,Ts,I^-i) = 


^ d 0 i ^ da 




(n!)^ gVE7^2VE+4 J^2w-2 27r 27r Jf 


’4.2 


4=1 


4=1 


'^ 4=1 


y y y j[2v,dv, 

-'1 j_i '^^Sj—2 j—2 j—1 ^ 

W 

X exp { - M{K{X, V) + L{B,T))} ■ [](a :,-1 - a;,- 2 )'( 6,-1 + 6 ,, 2 )' • A{X,B,V,T). ( 8 . 6 ) 

4=1 

We will also need the following facts 
w 

n 1(2^4 ,1 -2^4,2)"(&4,i +&4,2)"I = \A{X,B,V,T)\ < 

4=1 


(8.7) 


if 


Xi,X 2 eS'^, bj^i = a++ o{l), bj ^2 =-a-+ o{l), tj=o{l), Vj = l,...,iV, 


which always hold in these types of vicinities. The first estimate in (8.7) is trivial, and the second follows 
from Lemma 6.1. 

Now, we approximate (8.1) in the vicinities. For any G L, we introduce the matrices 


£+W ■■= 


0 e*’’ 

3"*’’ 0 


£-{^) ■■= 


0 


,-ii9 


0 


Then, with the parameterization above, expanding Xj in (3.22) and Tj in (3.24) up to the second order, 
we can write 

Xj =D^ + -^D^dia.g{xjf,Xj^2) + ^RJ, x = ±,+,-, 

Tj = I -\ —r^=£’_|_(cTj) H— tjR^- 

y/M M ^ 

For XT = ±, we also expand Vj in (3.24) up to the second order, namely, 


( 8 . 8 ) 










55 


We just take (8.8) and (8.9) as the definition of R* and K". Note that Rj is actually >f-dependent. 
However, this dependence is irrelevant for our analysis thus is suppressed from the notation. It is ele¬ 
mentary that 


I I I I 

I \-^j 11 max 


“k 2)1 ||.Rj||max 






0 (u|). 


( 8 . 10 ) 


Here || • ||max represents the max-norm of a matrix. 

Recall the facts (5.10) and (5.19) 

ML{B,T) = M (£++(bi) +i__(b2)) +M^s(S,r), 

mr:(x,h) = m(-£++(xi)-£+_(x 2 )) +Mes{x,v). ( 8 . 11 ) 

In light of (5.17)-(5.19), we can also represent MK{X, V) in the following two alternative ways 

MK(X, V)=M (-i++(xi) - i++(x 2 )) + Mls{X, V) + M{K{D+, I) - K{D±,I)), (8.12) 

Mk{X, V) = M (-i+_(xi) - i+-(x 2 )) + Mis[X, V) + XI{K{D_J) - K{D±,I)). (8.13) 

We will use three representations of Mk{X, V) in (8.11), (8.12) and (8.13) for Type I’, H and HI vicinities 
respectively. In addition, we introduce the matrices 


A+ := (1 -I- a\)I -I- a\S, A_ := (1 -|- a_)/ -I- a_S. 

Then, we have the following lemma. 

Lemma 8 . 2 . With the parametrization in (8.8), we have the following approximations. 
• Let bi,b 2 G and ||bi||oo, ||b 2 ||oo = o{y/M), we have 


(8.14) 


M (£++(bi) +i__(b2)) = -b'lA+bi -G -b'2A_b2 + R^ = o( 
• Let = ± and Xi,X 2 G and ||xi||oo, ||x 2 ||oo = o{VM), we have 


'Sa=l,2 11 ^' 


alls 


y/M 


). (8.15) 


1 »/ A A , A A , r,x r,x ^7^0 = 1,2 11^0 


XI (^-£++{xi) - i+_(x 2 )) = -x'^A+xi -G -X 2 A_X 2 -G = O 

• In the Type II vicinity, we have 

XI (- 

• In the Type III vicinity, we have 


(-i++(xi) - £++(x 2 )) = -x'^A+xi -G -x' 2 A+X 2 -G R+, R\ = 0{ 


Vm 


(%)■ 


). (8.16) 


y/M 


0t 


XI ^-£+_(xi) -i+_(x 2 )) = -x'^A_xi -G -X 2 A_X 2 + Rf, Rf = o(^-^=y (8.18) 

Here Rf R^, and Rf_ are remainder terms of the Taylor expansion of the function £(a) defined in 
(5.6). 


Remark 8 . 3 . Here we stated (8.15) and (8.16) in the domains much larger than the Type I’ vicinity for 
further discussion. In addition, the restriction ||ba||oo and ||xa||oo for a = 1,2 is imposed to avoid the 
ambiguity of the definition of the logarithmic term in the function £{a). 


Proof. It follows from the Taylor expansion of the function £(a) easily. 


□ 
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Then, according to (8.11)-(8.13), what remains is to approximate M£siB,T) and Mis{X,V) in the 
vicinities. Recalling the definition in (5.6) and the parameterization in (8.2), we can rewrite 


1 


Mis{B,T) = - - SkljC 


iaj 12 


j,k 

X ^a+ — a_ + 

. (a+-a-)^ 


a+4,1 ~ o,— bj,2 

y/M 




a+ — a_ + 

gia, |2 


a+6fc,i - a-hk,2 

y/M 


j,k 


We take the above equation as the definition of R^’^. Now, we set 

Tjp := tj cos aj , Tj ^2 ■= tj sin Uj, V j = 2,..., W 
and change the variables and the measure as 

( 4 , o-j) ^ (r,-,i, T,-2), 24d4^ ^ idr,-idTj,2- 

In the Type I’ vicinity, we can do the same thing for M£s{X, R), namely. 


M£s{X,V) =: 4+ + RY, 


j,k 

where R'^^ is the remainder term. Then we set 

Vj^i := Vj cos 4 , ■= Vj sin 4 , V j = 2,..., IT 

and change the variables and measure as 

• 1 

{Vj,0j) (wyi,Uy2), 2vjdvj^ -dz;yidz;y2- 

ZTT TT 

Now, we introduce the vectors 

a — (^2,a 5 • ■ • 5 '^a — (^2,a; ■ • ■ ; '^W,a)z ^ — 1,2. 

With this notation, we can rewrite (8.19) and (8.21) as 

M£s{B,T) = -{a+-a.f Y. <S^^^ra + R^ 


ot.h 


M£s{X,V) = -{a+-a.f Y v^ + R^■ 

a=l,2 


(8.19) 


( 8 . 20 ) 


( 8 . 21 ) 


( 8 . 22 ) 


(8.23) 


According to (8.20) and (8.22), we can express (8.5) as an integral over b, x, r and w-variables. 
However, we need to specify the domains of r and w-variables in advance. Our aim is to restrict the 
integral in the domains 


Ta G Ts, Da G T 5 , 0=1,2. 

Taking t for instance, we see that 

{tj - tkf < = (r,-,i - Tk,if + (rj,2 - Tk,2f, 

which actually implies 


(8.24) 


XaGTs for a = l,2=^i eis. (8.25) 

However the reverse of (8.25) may not be true. That means, (8.24) is stronger than (t,v,cr, 0) G T 5 x 
Tg X X To show the truncation to (8.24) from {i,v,cr,6) G Tg x Tg x x 
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is harmless in the integral (8.5), we need to bound and in terms of and v'^ 5(1) fJa 

respectively. More specifically, we need the following lemma. 


Lemma 8.4. In the Type /’ vicinity, we have 
05 


\R"±"\ < 


o{^) E (- +o(^) E (- + o{:^y (8,26) 


'02 


a=l,2 a=l,2 

07 


•y/M 


a—1^2 a—1,2 

Proof. Since the proofs of these two bounds are nearly the same, we only state the details for the first 
one. By (8.4), we have 

. 02 , 


s, = 1 + O 


(m)’ = ^. 


Then it is not difficult to see that 


|s,4e“^'= - Sfe 4 e*'^" p = 


— tjC^'^^ + O 
03^ 

a=l,2 


\mJ 


oj tfcc 

= H 47.a-T-fc.a)^+0(^) ^ 4^.a - Tfe,^ | + O ( 
0 ^ 1,2 ^ " 

Now, by the fact from Definition 8.1 

||ba||oo < ||ba ||2 = 0(05) for 0=1,2, 

we have 


M 2 


( 


a+ — a_ + 


a^_4,i ~ 0'—bj,2 


y/M 


a+ — 0- + 


a+bk,i — a-bk,2 
y/M 


) -(a+-a _)2 + 0(^) 

V j 


Substituting (8.27) and (8.29) into (8.19) yields 
05 
^y/M 


a—1,2 j,k 


a—1,2 j,k 


(8.27) 

(8.28) 

(8.29) 


I'R*’''! < o(-^) -Tfc,a)^ +o(-^) -T-fc.al +0 

07 


.M^J 




a—1,2 a—1,2 

Here we used Cauchy Schwarz inequality and f. Sjk = 0(W) = 0(0) in the second step. The bound 
for can be proved analogously. Hence, we completed the proof of Lemma 8.4. □ 

Roughly speaking, by (8.23) and Lemma 8.4 we have 

Mis{B,T) = -{a+-a.+o{l)Y E + E {-+ o{1), 

a—1,2 a—1,2 

Mes{X,V)=-{a+-a.+o{l)Y E «'a*5(l)t7, + o( E ( - + o(l). 

a—1,2 a—1,2 

Then it is obvious that if one of Ti, T 2 , Vi and V 2 is not in Ts, we will get (7.26). Hence, using (8.7), we 
can discard the integral outside the vicinity, analogously to the proof of Lemma 7.7. More specifically, in 
the sequel, we can and do assume 


Ti,T2,Vi,V2 e Ts- 

Now, plugging (8.30) into (8.26) in turn yields the bound 


(8.30) 


\RY\ = 0{ 


/ 05 04 n 

WM ^ M/’ 


(8.31) 
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By the discussion above, for the Type I vicinity, we can write (8.5) as 
, M 2 ^ 

'Y'O 'Y'^ T'^ T'T' _ 


2'^I(T';, T'L, , T!i, Ts, Ts) = 


VF 


n h n 


(n!) 24 ^M ^+2 .Ay. 1^1 

» TV r ^ r 

, n L n ‘^®.-2 1 n L n ^^.-2 / n L n *^^.■2 


w 


w 


w 


w 


i=i 


i=i 


/Ts 


i=2 


X exp I - ib(A+bi - ib 2 A._b 2 - • exp | - ix^A+xi - ix 2 A_X 2 - i?±| 

X exp|(a+- a_)2 ^ • exp |(a+- a_)2 ^ 


w 


w 


X n exp - ^J,2f(hi + h 2 ? ■ A(X, B, y,T) + 0(e-®), (8.32) 

i=i i=i 

where the error term stems from the truncation of the vicinity (t,v, cr, 0 ) £ Ts x Ts x x 

to (xi,T2,'Di,D 2) G ts X ts X ts X ts. 

Now, for the Type II and III vicinities, the discussion on is{B,T) is of course the same. For £s{X, V), 
we make the following approximation. For the Type II vicinity, using the notation in (5.21), we can write 




3,k 


M 


))yVM~ y/M ^ M )) 


= ■ -Xj^2){.Xk,l -Xka) + RT- 


i,fe 


It is easy to see that 


RY = Oi 


|iilli + l|x 2 |ii' 


Vm 




(8.33) 


(8.34) 


Combining (5.15), (5.18), (8.17) and (8.33) we obtain 


M{K{X, Y) - K{D+, /)) = ix'iA+Xi + ^x' 2 A+X 2 - ^ i - x,,2){xk,i - ife. 2 ) + i?+ + 


1 V 


2 

= ii'A" X + R1 + i?"’". 


j,k 


where 


X := (x(,X 2 )', A+ := (1 + a+)/ 2 w +a+§’', 
and recall that is defined in (5.22). Analogously, for the Type III vicinity, we can write 


where 


M{K{X, V) - K{D_,I)) = ix'A" X + R^_ + 


0 “ 

r := (l + a^)J 2 iv+aiS", i?t" = o(;^). 


(8.35) 

(8.36) 

(8.37) 


Consequently, by (8.12) and (8.13) we can write (8.6) for k = +, — as 

r '1 I— 

I{TXX -, T^, T-, Ts, I^-^) = exp {m(/C(7^±, I) - K{D^, I ))} • 


w 


(W, 


w 


w 


w 


/LW-1 


n ^ n 1^ n l n n l n ^^2, 


(n !)2 sw^3W+3 

w „ w 

2 


3=2 


3=2 


'^3=1 


3 = 1 


3 = 1 


3 = 1 
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iTs 


W W 1 . o 1 „ „ 

n /. n • exp I - -b'^A+bi - -b 2 A_b 2 - 

i =2 j =2 


X exp{ --x'A^x-•exp{(a+- a_)^ ^ 

,1 + 
y/M 


w 


w 


a=l,2 
\2 


X exp|i ^^’V^-^’^ I • ]^(a;j,i - a;j, 2 )^(&i,i + fej. 2 )^ ■ A(X,B, y,T) + 0(e ®). (8.38) 

i=i X i=i 

8.2. Steepest descent paths in the vicinities. In order to estimate the integrals (8.32) and (8.38) 
properly, we need to control various remainder terms in (8.32) and (8.38) to reduce these integrals to 
Gaussian ones. The final result is collected in Proposition 8.6 at the end of this section. As a preparation, 
we shall further deform the contours of b-variables and x-variables to the steepest descent paths. We 
mainly provide the discussion for the b-variables, that for the x-variables is analogous. 

For simplicity, in this section, we assume 0 < E < \f2 — k, the case —y/2 + K < E < 0 can be discussed 
similarly. We introduce the eigendecomposition of S as 

S=US[J'. 

Note that U is an orthogonal matrix thus the entries are all real. Now, we perform the change of coordinate 

Ca — (ci,a; ■ • ■ ; ^W,a') ■— U b^. Cl — I 5 2. 


Obviously, for the differentials, we have 

w w 

j=i i=i 

In addition, for the domains, it is elementary to see 


a = 1 , 2 . 


ba e T 


Ca G T, 0 = 1,2. 


(8.39) 


Now, we introduce the notation 


7 ;- 


1 + 0 ^+ a^Aj(S) 


7, := 


1 + 0 ?.+ atAj{S) 


(8.40) 


'‘+ ^ 

and set the diagonal matrices 

]D>+ := diag( 7 )^,..., 7 ,j;», D. := diag( 7 f,..., 7 ,^). 

By the assumption 0 < E < y/2 — k and (1.4), it is not difficult to check 

17 / 1 ^ 1 , 1771 ^ 1 , arg 7 +G (-^,0], arg7“G[0,^), Vj = l,...,lF. 

With the notation introduced above, we have 

b']^A+bi = c'^D^^Ci, b2A_b2 = C2]D>I^C2. 

To simplify the following discussion, we enlarge the domain of the c-variables to 
c, G Too = Too(e):= [-05,05]’^, o = l,2. 

Obviously, T C Too- It is easy to check that (7.26) also holds when c^ G Too \ T for either a = 1 or 
2, according to (8.39), thus such a modification of the domain will only produce an error term of order 
0(exp{—0}) in the integral (8.32), by using (8.7). 

Now we do the scaling 


Cl —>■ ID+Ci, 


C2 D_C2. 
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Consequently, we have 


thus 


bi = UID)+ci, b2 = UID)_C2, 


b'lA+bi = ^ b'2A_b2 = X] S'2- 


(8.41) 

(8.42) 


Accordingly, we should adjust the change of differentials as 

w w w w 

Ild^j.i ^detro+ • Y[ dcj-i, n d 6 j ,2 —t det]D)_ • n dCj, 2 - 


i=i 


i=i 


i=i 


i=i 


III cidciitioii, tlic domS/iii of should bo chcnigod from Tfoo n]!li where 
and that of C 2 should be changed from Too to 1^=1 > where 

J 7 := ( 7 -)-M- 0 ^e^]- 


By the fact detID)+D_ = 1/ -^det A_|.A_, we can write (8.32) as 

/Vf 2 

2^I(T^, T7, , T7, Ts, Ts) = 


(n!)24’^7r4M^+2 y^det A+A_ 


w 


w 


w 


w 


Ifd^ 


'nr=i-^jT=i 
w 


id 




n /. n /. n 


i=i 


^i=i 




w r ^ r ^ r ^ ^ 

j=2 •'^S j—2 j—2 2 j=l \ M 


X exp 


®i =2 i=i 

U/A ^ 1 ^/ 


{ - - ^ 1 - - ^ ^2 - • exp I - -x;A+Xi - -x' 2A_X2 - i?±| 


X exp |(a+- a_)^ ^ r'^S'^^Va - ■ exp |(a+- a_)^ ^ 


a=l ,2 


a=l ,2 


W 


X [](x,-1 - X,.2)^(6,.1 + 6,.2)2 • A{X,B,V,T) + 0 (e-®). 
i=i 


(8.43) 


In (8.43), all the b and b-variables in the integrand should be regarded as the functions of the c-variables, 
see (8.41). 

Now, we consider the integrand as a function of c-variables on the disks, namely, 

c,,i e 0+ :={zGC: I^I < e^ld+l”^}, c ,-2 £ O" := {z G C : |z| < e^lq-r'}- 

For Cl G n]^i ^2 € nj!li ®7’ (8.40) and (8.41) we have 


|bi| 


|b 2 ||oo < 0{Q). 


(8.44) 


Here we used the elementary fact ||C/a||oo < v^tHHaHoo for any a £ and and unitary matrix U. Then, 
we deform the contour of from jF to 

(-e+)ul;ue; 

for each j = 1,..., W, where 

L+ := K n 0+, E+ = {z £ 50+ : 0 < argz < - argy+j . 







It is not difficult to see that 
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ReCj 1 > 0 , for cj^i € (-S^) U 
by using (8.40). Consequently, by (8.42), we have 


w 

exp { - it.; A+bi} = exp { - i ^ ^} < 0( 

i=i 


.-e\ 


Then using (8.7), we can get rid of the integral over and analogously to the discussion in Section 

7. Similarly, we can perform the same argument for C 2 . Consequently, we can restrict the integral in 
(8.43) to the domain 

w w 


J =1 

tW tt- 


i=i 


So we can assume that n 7 =ij;and n7=iJ- are replaced with n7=i L 7 and L^- respectively in 
(8.43). 

By (8.15), (8.44) and the fact ||a ||3 < ||a||oo||a||| for any vector a, we see that 


|i?7 < C- 


|bi||oo + ||b 


2 Moo 


|bi 


VM VM-M2 . 1102112 ; < 

for some positive constant C, where in the last step we also used the fact that ||ba ||2 = 0 (||ca|| 2 ) for 
0 = 1 , 2 , which is implied by (8.41) and (8.40). Consequently, we have 

exp{ - i||ci ||2 - ^I|C 2||2 -^'’1 = exp I - (i + o(l))||ci ||2 - (i + o(l)) Ucallaj. 

This allows us to take a step further to truncate Ci and C 2 according to their 2-norm, namely 

Ci,C2eT. (8.46) 


0 


|ci ||2 + IIC2II2) 


(8.45) 


Similarly to the discussion in the proof of Lemma 7.7, such a truncation will only produce an error of 
order exp{—0} to the integral, by using (8.7). 

Now, analogously to (8.41), we can change x-variables to d-variables, defined by 


di = (dpi,.. ■,dw,i) ■■= D_^^U'xi, 

Thus accordingly, we change the differentials 
w w 

nd^pi^detD+.J]^ ddyi, 

1=1 1=1 

In addition, like (8.46), we deform the domain to 


d2 = (di,2, ■ • ■,dw,2) := D_^U'x2. 


w w 

nd^p2^detD_ ddj,2- 
1=1 1=1 


di,d2 G t. 


Finally, from (8.43), we arrive at the representation 


2^I(t7, t 7, T7, T)i, Ts, Ts) 
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X exp|(a+- a_)^ ^ • exp |(a+- a_)^ ^ 


a=l,2 


a=l,2 


W 


X J](x,- 1 - x,,2)^(6,.1 + b,,2f ■ A(X,B,V,T) + 0 (e-®), 
i=i 


(8.47) 


in which x and x-variables should be regarded as functions of the d-variables, as well, b and b-variables 
should be regarded as functions of the c-variables. 

Now, in the Type II and III vicinities, we only do the change of coordinates for the b-variables, which 
is enough for our purpose. Consequently, we have 

I(T'+, T'_, T^, rz, Ts, I^-^) = exp {M{K{D^, I) - K{D^, /))} • 


(n!)2 s'^t:^w+3 


J detA+A_ Jiw-i -H- 2tt Jjw-i -*y' 

V J — Z 1 — 




w 


w 


w 


n L n dcyi / n 


i=2 


i=i 
w 


i=i 


w . w r ^ r ^ ^ °_L° 

n <ii„i / n / n ■i-w n n “p 

j — 1 i—\ ^ * 0—2 ^ 0—2 o—\ * 


i=i 


i=2 


1=1 


xexp{-i||di||2-i||d2||2-i?'’}-exp|(a+-a_)2 


'{ 


1 V. 


w 


X exp ■{ — -X A^x — RZ — R 


<V^X 


} - ^^■. 2 )'(&,.i + bj, 2 f ■ A(l, B, I/, T) + 0(e-®). 

i=i 


By (8.40) and (8.41), it is easy to see 


Similarly, we have 


Ci,C2 e T =4 p bi,b2 € T 


di,d 2 e T =» Xi,X2 e T. 


(8.48) 


(8.49) 


(8.50) 


We keep the terminology “Type T, II and III vicinities” for the slightly modified domains defined in terms 
of c, d, r and 7 ^-variables. More specifically, we redefine the vicinities as follows. 


Definition 8.5. We slightly modify Definition 8.1 as follows. 

• Type r vicinity: Ci,C 2 ,di,d 2 £ T, ti,' r 2 , Wi ,^2 £ Ts- 

• Type II vicinity: Ci, C 2 , Xi, X 2 £ T, Ti,'r 2 £Ts, Vj £ U{2) for all j = 2,... ,W, 

where 'k-variables are defined in (8.8) with x = +. 

• Type III vicinity: Ci, C 2 , Xi, X 2 £ T, ti,T 2 £Ts, Vj £ U{2) for all j = 2,... ,W , 

where x-variables are defined in (8.8) with x = —. 

Now, recall the remainder terms i?^, RZ, -R+ and Rf_ in Lemma 8.2, and RZ^ in (8.23), Rf)^ 
in (8.33) and in (8.36). In light of (8.49) and (8.50), the bounds on these remainder terms are the 
same as those obtained in Section 8.1. For the convenience of the reader, we collect them as the following 
proposition. 

Proposition 8.6. Under Assumptions 1.1 and 1.13, we have the following estimate, in the vicinities. 
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[ 11 ] 


w = o 



- of 

Ri - of 

R^ - of 

V\/M/ 





z 0 .y 

J^V,X _q1 \ 




(m): RY=0 


Proof. Note that, (i) can be obtained from (8.26), and (ii) follows from Lemma 8.2, and (iii) is implied 
by (8.34) and (8.37). Hence, we completed the proof. □ 

Analogously, in the vicinities, ||bi|||, ||b2|||, ||xi|||, ||x2|||, ||t||oo and ||v||oo are still bounded by 0. 

9. Integral over the Type I vicinities 

With (8.47), we estimate the integral over the Type I vicinity in this section. At first, in the Type 
F vicinity, we have ||xa||oo = 0{&i) and ||ba||oo = 0(0^) for a = 1,2. Consequently, according to the 
parametrization in (8.2), we have 


j,i - Xj ,2 - a+ - a- +o(^-^=j, 


/ 02 \ 

bj,i + 6 j ,2 = a+ - a- + Oy--j= J , 


which implies 


w „| 




(9.1) 


(9.2) 


Hence, what remains is to estimate the function A(X, B, H, T). We have the following lemma. 

Lemma 9.1. Suppose that the assumptions in Theorem 1.14 hold. In the Type I’ vicinity, for any given 
positive integer n, there is Nq = NqIji), such that for all N > Nq we have 

02 p(yCo 


\A{X,B,V,T)\ < 


■ I det A+p • det(5'^^^)^. 


M{Nr])^+^ 

for some positive constant Co and some integer I = 0 ( 1 ), both of which are independent of n. 

With (8.47), (9.2) and Lemma 9.1, we can prove Lemma 5.8. 

Proof of Lemma 5.8. Using (8.47), (9.2), Lemma 9.1, Proposition 8.6 with (5.38), the fact detA_|_ = 
det A_ and the trivial estimate 

for sufficiently large constant Co, we have 


2 ‘^|I(T%T^,T^,T^,Ts,Ts)| < 
„ w „ w 


(Xg) 

NCo 


1 


(Ng)'^ (27r2)2iv 

w „ w „ w 


• det(5'^^^)^ • (a+ — 

w r ^ 

RdcYi y ndcY2 y iid^i y 11^42 y ndu.i y ndu.2 

i=i •^^1=1 •^^®i=2 '^^^ 3=2 

w IV r . 

/. exp<^--(llcill^ + IIC2II2 + lldill^ + ||d2||2) [ 

X exp |(a+ — a_ f + 45(^42 + + 45(1)1^2) } + 0 (e-®). 

Then, by elementary Gaussian integral we obtain (5.40). Hence, we completed the proof of Lemma 
5.8. □ 


/Ts 
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The remaining part of this section will be dedicated to the proof of Lemma 9.1. Recall the definitions 
of the functions A(-), Q(-), P(-) and F(-) in (3.31), (4.3), (4.4) and (4.5). Using the strategy in Section 6 
again, we ignore the irrelevant factor Q(-) at the beginning. Hence, we bound P(-) and F(-) at first, and 
modify the bounding procedure slightly to take Q(-) into account in the end, resulting a proof of Lemma 

9.1. 

9.1. P(X,i3, U, T) in the Type I’ vicinity. As mentioned above, we should always regard b or fa- 
variables as functions of c-variables, regard x or x-variables as functions of d-variables. Our aim, in this 
section, is to prove the following lemma. 

Lemma 9.2. Suppose that the assumptions in Theorem 1.14 hold. In the Type I’ vicinity, we have 

w2+7C)2 

P{X,B,V,T) < ^ |detA+|^det(,gW)^ (9.3) 

Before commencing the formal proof, we introduce more notation below. In the sequel, we will use the 
notation 

kj = kjiX, B, V, T) := |xyi| -F + l^j.2| + I'^jl + = 0(0), (9-4) 

where the bound holds in the Type T vicinity, according to the facts ||xa||oo = 0(0^), ||ba||oo = 0(05) 
for 0 = 1,2, ||t||oo = 0(0) and ||v||oo = 0(0). 

Recalling (6.17) with zuj defined in (6.13) and in (6.15). Now, we write 

= exp { - Mlogdet (l -F } 

where 

A,., := Tr{v;X-W,n,T-^B-%E,Y. (9.6) 

The second step of (9.5) follows from the Taylor expansion of the logarithmic function. The expansion 
terminates at the 4th order term since Aij is a homogeneous polynomial of flj and -variables with 
degree 2£, regarding all the complex variables as fixed parameters. Now, we expand the first factor of 
(9.5) around the Type T saddle point, namely 

exp { - TrV;X-^V,n,T-^B-%E ,} =: exp { - TrDg^n,Dg% } exp { - ^ A, }. (9.7) 

We take (9.7) as the definition of Aj, which is of the form 

4 

Aj = Pj.a./3 • 

1 

for some function pj,a,i 3 of x, fa, v and t-variables, satisfying 

Pj,a,p = 0{kj), V a,/? = 1,... ,4. (9.8) 

One can check (9.8) easily by using (8.8)-(8.10). Analogously, we can also write 

4 i 

= 'y ^ Pt,j,a,/3 J_ ) Ck ■= {o^l, ■ ■ ■ t Oil), /3 .= (/3i, . . . ,/3f) , (9.9) 

ai,...,ai, j—1 

where 

= 0(1), V = 2,... ,4;ai,... ,a^,/3i,... ,/3£ = 1,... ,4. (9.10) 



=: exp { - TrV;X-W,n,T-^B-%E,} exp { 
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The bound on in (9.10) follows from the fact that all the 1^-, Tj, ^ and ^-entries are 

bounded in the Type F vicinity, uniformly in j. Consequently, we can write for j ^ p,q 


exp ■ 


V-'” f,=2 e=l 

9 i....,/ 3 f = l 


In a similar manner, we can also write for k = 


exp« 


I ^ V I'^fc = Po(‘) ( ^ ^ ^ ^ ^ V WA:,ai^fc,/3i j 5 (9-12) 

£—2 ^ ^—1 ai,...,ai, i—i / 


ai 


where po(‘) = detXfc/detilfe, which is introduced in (6.15), and ^ej^a^p is some function of X, B, V and 
T-variables, satisfying the bound 

qo.«,/3 = O((l + %)0, V£=l,...,4, j = l,...,VF. (9.13) 

Obviously, we have po(‘) = 0(1) in Type F vicinity. 

Now, in order to distinguish i, ot and /3 for different j, we index them as (.j, otj and /3j, where 

In addition, we introduce the vector 

i:= {£i,...,£w), a = a{£):={ai,...,aw), 0 = ^{£) := {13^,..., (3^,). 

Let ||.^|i = ^0 be the 1-norm of £. Note that a and $ are ||£||i-dimensional. With these notations, 

using (6.17), (9.5), (9.7), (9.11) and (9.12) we have the representation 


w 


1 - 1 ^ 


'P{n,E,X,B, V,T) = po{Xp,Bp)fo{Xq,Bq) ■ exp | - y^SjkTrnjEk - 

j,k j=l 

W W 

1 + M-^ nqc.F«,-./ 3 j • 

FgI 0,41'^, 5,dG[1.4]M4li i=l 1=1 i=l 

s.t.||d|i>l 

where we made the convention 

0 

90 , 1 , 0,0 = 1; 11 ^l,“j,>?l,/^ 3 ,i “ J = 1,...,1F. 

i=l 

According to (9.13) and (9.15), we have 


(9.14) 


w 


w 




1=1 


1=1 


In addition, we can decompose the sum 


4W 


It is easy to see 


Moreover, it is obvious that 


E =E E ■ 

Fg[0,4]“', Fg[0,4]1^, 

s. t.||i||i>l s. t.|| 7 ||i=m 


H{/g lO,4r:|Kl|i=m}< 


tl{d,/9 G [i,4iib1ii} = leii^lib 


(9.15) 


(9.16) 


(9.17) 


(9.18) 

(9.19) 
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Therefore, it suffices to investigate the integral 

» w w t-i 

•“ / drid^exp I — ''^^SjkTrflj^ik — IT 

hk j=l i=li=l 

for each combination {i,a,(3), and then sum it up for {i,a,P) to get the estimate of P{X, B,V,T). 
Specifically, we have the following lemma. 


Lemma 9.3. With the notation above, we have 

= IKlli =0 or 1 . (9.20) 

Moreover, we have 

|detA+|2det(.SW)2(||^l|i-l)!(2lT^)(ll''ll-i), if \\i\U>2. (9.21) 

We postpone the proof of Lemma 9.3 and prove Lemma 9.2 at first. 


Proof of Lemma 9.2. By (4.4), (9.14) and (9.20) and the fact that po(‘) = 0(1); we have 

w 

IP(X,B,V,T)I<C E (9.22) 

^eIO,4lW, a,/3e[l,4]ll4li 1=1 

s. t.||d|l>2 

Substituting the bounds (9.16), (9.19) and (9.21) into (9.22) yields 
IP(X,B,V,T)I < |detA+pdet(S'(i))2 

w 

X • (ll^li - l)!(2WT')0^1li-i) • ]^(1+ Kj)^A (9.23) 

^ 610 , 41 ^, 1=1 

s. t.||d|l>2 

Now, from (9.4) we have 


w 

]^(1+ Kj)^i < 011^11% (9.24) 

1=1 

which can absorb the irrelevant factor Using (9.17), (9.18) and (9.24), we have 

w 

• A/-^ • (ll^lli - l)!(21UT')0^1li-i) • J|(l + Kj)^i 
^6 [ 0,41 w , 1=1 

s.t.||d|i>2 

41V /,TTr\ 

-E( J • e"' 

m=2 V ^ 1 

< • 0"' • = o(- — —j, (9.25) 

m—2 

where the last step follows from (5.33) and (5.38). Now, substituting (9.25) into (9.23), we can complete 
the proof of Lemma 9.2. □ 


Hence, what remains is to prove Lemma 9.3. We will need the following technical lemma whose proof 
is postponed. 


Lemma 9.4. For any index sets I, J C {1,..., W} with |l| = |J| = m > 1, we have the following bounds 
for the determinants of the submatrices of S, A+ and A_ defined in (8.14). 
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• For and we have 

|det(A+)('IJ)| ^ ^ |det(A_)('IJ)| ^ 

I det A+1 ~ ’ I det A_ | “ 

• AorS'ClJ), we have 


I det 5(1) I 


< (m 


l)!( 2 W^i')(’"-i). 


Proof of Lemma 9.3. Recall the definition in (6.20). Furthermore, we introduce the matrix 

H = (a+^A+) © 5 0 5 © (aI^A_). 

Using the fact a+a_ = — 1, we can write 

w 

-'^SjkTrQjEk -^TrD^^rtjDg^Ej = -OHS'. 
j,k j=l 

Now, by the Gaussian integral of the Grassmann variables in (3.2), we see that 

|%3_^| = |detH('IJ)| 

for some index sets I, J C {!,..., 41U} determined by a and /3 such that 


|l| = |J| = IKlli- 


(9.26) 


(9.27) 


(9.28) 


(9.29) 


Here we mention that (9.29) may fail when at least two components in a.j coincide for some j. But 
/3 “ ill i'iii® case because of = 0 for any Grassmann variable y. 

Now, obviously, there exists index sets U, Jq, C {!,..., W} for a = 1,..., 4 such that 

HCII) = (a;2A+)('d-ii) © ^Cdi^) © 5('db) 0 (al2A_)('d-'4)^ ^ ||^| = ^ |j^| = ||£j|^. 

a OL 

It suffices to consider the case 


|U| = |J„|, Va = l,2,3,4. 

Otherwise, detElCIi) is obviously 0, in light of the block structure of H, see the definition (9.28). Now, 
note that, since det 5 = 0, we have 

detHCH) = 0, if ||£^|i=0,l. 

In addition, if ||^||i = 2, by using (9.30) below, one has 

r |detA+A_|det(5(i))2, if lb] = [Ul = IJ 2 I = IJ 3 I = 1, 

|detH('l-')| = 

[ 0, otherwise. 

For more general by Lemma 9.4, we have 

IdetHClDl < |detA+A_|det(5(i))2(||£l|i - l)!(2IFi')(lldd-i). 

Then, by the fact | det A_|.A_| = | det A+p, we can conclude the proof of Lemma 9.3. □ 

To prove Lemma 9.4, we will need the following lemma. 

Lemma 9.5. For the weighted Laplacian S, we have 

det5(*d) = (-I)J-Met5(d, V i, j = 1,..., IF (9.30) 


Remark 9.6. A direct consequence of (9.30) is 

det 5 ( 1 ) = ... = detS‘'^\ 


(9.31) 
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Proof of Lemma 9.5. Without loss of generality, we assume j > i'm. the sequel. We introduce the matrices 

Ij-i-l 


Pij ■— Ii—1 


It is not difficult to check 


1 


) Iw-j, E,-.= I- 2e,e* - ^ e^e*. 


sm = 

Then, by the fact det PijEj = (—we can get the conclusion. 


(9.32) 

□ 


Proof of Lemma Q.f. At first, by the definition in (8.14), (1.4) and the fact Rea^ = Rea^ > 0, it is easy 
to see that the singular values of A+ and A_ are all larger than 1. With the aid of the rectangular 
matrix (A_|_)*^'l®^ as an intermediate matrix, we can use Cauchy interlacing property twice to see that 
the fc-th largest singular value of (A_|_)^'l'') is always smaller than the fc-th largest singular value of A_|.. 
Consequently, we have the first inequality of (9.26). In the same manner, we can get the second inequality 
of (9.26) 

Now, we prove (9.27). At first, we address the case that I fi J 0. In light of (9.31), without loss 
of generality, we assume that 1 € I fl J. Then is a submatrix of Therefore, we can find two 
permutation matrices P and Q, such that 

A B 
C D 




where D = Now, by Schur complement, we know that 


|det5'('l-')| 




= Idet(A-BD-^C) 


I dots'll)! 

Moreover, (A — BD~^C)“^ is the (|l| — 1) by (|l| — 1) upper-left corner of 

That means det / det is the determinant of a sub matrix of (5'^^^)“^ (with dimension |l| — 1), up 
to a sign. Then, by Assumption 1.1 (hi), we can easily get 

|det5'('l-'V|det5'(^)| < (|l| - l)!W(l'l-i)'>'. 

Now, for the case I fl J = 0, we can fix one i € I and j G J. Due to (9.30), it suffices to consider 

det S'CI-') 

By similar discussion, one can see that (9.33) is the determinant of a sub matrix of (S'^*d))-i with 
dimension |l| — 1. Hence, it suffices to investigate the bound of the entries of From (9.32) we 

have 




(9.34) 


Observe that 




Ij-i-l 


Iw-j , 


E-^=E,. 


Then, it is elementary to see that the entries of (5^1-?)) ^ are bounded by 21F^, in light of (9.34) and 
Assumption 1.1 (iii). Consequently, we have 

< (||| - l)i(21T'>')l'l"^ 

|det5(®h-)| - ^ ’ 

which implies (9.27). Hence, we completed the proof of Lemma 9.4. □ 
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9.2. P{X, B ,V,T) in the Type I’ vicinity. Neglecting the yt^l and -variables in Q(-) at first, 
we investigate the integral F(X,i3,y,T) in the Type F vicinity in this section. We have the following 
lemma. 


Lemma 9.7. Suppose that the assumptions in Theorem 1.14 hold. In the Type I’ vicinity, we have 

1 




(9.35) 


.(iV7?)"+2, 

Recalling the functions G{B,T) and F(X,fo) defined in (6.23) and (6.24), we further introduce 

<G{B, T) = exp {(a+ — a-)Nr]}G{B, T), F(W, V) = exp { — (a+ — a_)fV77}F(X, V). (9.36) 

Then, we have the decomposition 

f(X,B,V,T) = G{B,T)¥{X,V). (9.37) 

Hence, we can estimate IF(X,fo) and G(R,T) separately in the sequel. 

9.2.1. Estimate of'¥(X,V). We have the following lemma. 


Lemma 9.8. Suppose that the assumptions in Theorem 1.14 hold. In the Type /’ vicinity, we have 


(9.38) 


Proof. Using (8.8) and (8.9), we can write 


X, = pfv;x,v,p^ = pfD±p^ + o(^), 


(9.39) 


where the remainder term represents a 2 x 2 matrix whose max-norm is bounded by Q/y/M. Using (9.39) 
and recalling N = MW yields 

w - 'QNr]'^ 


exp [MrjJ^TrX.j] = exp |iV?7rrP*D±Pi j|(l -)- 


i=i 


(9.40) 


Substituting (9.40) into (3.19) and (6.24), we can write 

F(X,U)= J d^(Pi)dXW exp{fV77TrP;L>±Pij} 


n 


ktp,q det^(wW) 


exp {iTrX™ JZ - j} 




k,i—p.,q 

Recalling the parametrization of Pi in (3.24), we have 

TrPfD±PiJ = (1 - 2v^)ia+ - a_). 

Consequently, we have 

F(X,U)= f dJfdl f vdv f — exp I — 2(a+— a_)iV7yz;^]- P - 

J 7 J tt I J det^C 


det^(xW) 

P exp [iTrXl^bz - j}(l + o(l)). 


k—p,q 
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Obviously, by the fact that -variables are all bounded and | det \ = 1 ior k = p,q, it is easy to see 
that 

/■i ^ 

|]F'(X, V)\ — ^ J ~ 2(a+ — a-)Nr]v^'^ = ■ 

Therefore, we completed the proof. □ 

9.2.2. Estimate of <G{B,T). Recall the definition of G{B,T) from (9.36), (6.23) and (3.20). In this 
section, we will prove the following lemma. 

Lemma 9.9. Suppose that the assumptions in Theorem 1.14 hold. In the Type /’ vicinity, we have 

1 




,(iV7?)"+iy 

Note that y^p \ y^l^ and t in the parametrization of Qi (see (3.24) ) are not bounded, we shall trun¬ 
cate them with some appropriate bounds at first, whereby we can neglect some irrelevant terms in the 
integrand, in order to simplify the integral. More specifically, we will do the truncations 

1 


t < 


and 


Accordingly, we set 


< {Nr])i. 


(9.42) 

(9.43) 


G{B,T) ;= [ — 

Jh 277, 

ANriy^ 


uWuWduWduW 






AN-q)-B 




[ 1 ] 


2Mt [ dcrWdo-W 5(Qi,T,i3,yW,wW), 

.70 

where we have used the parameterization of wdl in (3.15). We will prove the following lemma. 


(9.44) 


Lemma 9.10. Suppose that the assumptions in Theorem I.I 4 hold. In the Type I’ vicinity, we have 

<G{B, T) = G(B, T) + 0(e-^') 

for some positive constant e. 

Proof. At first, by (6.29)-(6.32), we have for any j, 


ReTrBjYl^b> . 

(s + t )2 


2 /W min{Re6j- 1 , Re6j, 2 } ^ 


k=p,q, 


(9.45) 


{Sj + tj)^ 

for some positive constant c, where the last step follows from the facts that Rebj^i, Rebj 2 = Rea+ -I- o(l) 
and tj = 0 ( 1 ) in the Type T vicinity. In addition, it is not difficult to get 

TrR,J= (a+-a_+0(^)) (1 + 2^2), Vj = l,...,W, 


which implies that 


o o./v 

MpY^TrBjJ = (a+ - a-)Nr] + 2(a+ - a-+ Npt^ + o(^-j=y 




Tm- 


(9.46) 


Note that the second and third factors in the definition of g{-) in (3.20) can be bounded by 1, according 
to (6.26). Then, as a consequence of (9.45) and (9.46), we have 

[ 1 ] I [ 1 ] 

g(a+-a.)7v^|^(.)| < c(yWy[i])"+3exp{-c'7V77t2}exp{ - 


(9.47) 
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for some positive constants C, c and d. By integrating and out at first, we can easily see that 
the truncation (9.42) only produces an error of order 0(exp{—iV®}) to the integral G{B,T), for some 
positive constant e: = £(£2) by the assumption 77 > in (1.15). Then one can substitute the bound 

(9.42) to the last factor of the r.h.s. of (9.47), thus 

""P{ ■ - ""P • 

We can also do the truncation (9.43) in the integral G{I3,T), up to an error of order 0(exp{—A^®}), for 
some positive constant e. Therefore, we completed the proof of Lemma 9.10. □ 

With the aid of Lemma 9.10, it suffices to work on G{B,T) in the sequel. We have the following 
lemma. 


Lemma 9.11. We have 




Proof of Lemma 9.11. Recall the parameterization of in (3.15) again. To simplify the notation, we 
set 

.„[i] [1] [1] 1 

tt’fe =<<> k=p,q. 

Similarly to (9.39), using t = o(l) from (9.42), we have the expansion 

Qf^T-^B,T,Q, = Q-^D^Q,+o[^). 

Consequently, we have 

/0A^r/N 

— Mr] ''^^TrBjJ = —Nr]{a+ — a_)(l + 2t^) + 0( 
i=i 

In addition, for k = p,q, using the fact ^ Sjk = 1, we have 


Bj = 


Vm - 


(9.48) 


J = TrQf^D±QiYpj 


0 


Tri?fcK 


[ 1 ] 


= Vk^ ((“+ - - a-(4^')^)) 

+ 4^'((a+ + -^TrRkY^ 

where is a 2 x 2 matrix independent of , satisfying 

||Rfe||max=0(l). 


[ 1 ] 
k ’ 


(9.49) 


(9.50) 


Observe that the term in (9.48) is obviously independent of wbl-variables. In addition, for k = p or q, 
we have 


iTryW JZ ={-v + iE{l - 2(4'l)2))yW, 

and for k,£ = p or q, we have 


Moreover, we have 


( KkwW)*)^^ (wW(wW)*)^J" = 




(9.51) 

(9.52) 

(9.53) 
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Substituting (9.48), (9.49) and (9.51)-(9.53) to the definition of g{-) in (3.20) and reordering the factors 
properly, we can write the integrand in (9.44) as 

exp{(a+ - a_)Nr]}g{-) 

= exp{m(crW - cr^)} • exp | - (a+ - a-)st ^ | 

k=p,q 

X exp { - TrRkYl'^^'^ ■ exp | | 

k=p,q 

X n ■ n exp { - ^^keVk^Ve^ (»fc’tt’™) } ■ exp { - 2Ng{a+ - a-)t'^} 

k—p,q k,£=p,q 

X n exp { - Vk'' ( (a+(Ufc')^ - a-(4^')^) + (a+ - a_)f + g - iE (l - 2(?;W)2^ ) | 

k—p,q 

(9.54) 

where the last factor is independent of the wPl-variables. Here, we put the factors containing and 
(Tg^^ together, namely, the first two lines on the r.h.s. of (9.54). 

For further discussion, we write for k = p,q 

TrRkYl^^ = 2/W (r+e^^ + + r^), (9.55) 

where r^, and tfc are all polynomials of and with bounded degree and bounded coefficients, 
in light of (9.50), the definition of Y^^^ in (3.14) and the parametrization in (3.15). 

Now, we start to estimate the integral (9.44) by using (9.54). We deal with the integral over and 
at first. These variables are collected in the integral of the form 

Ia{hj2) ■■= dcr^dcrW exp{i(n+ £i)crW} ®xp { - i(n+ £2)0-^} ®xp I -^ 

^ ^ k—p-,q 


X exp I — (a+— a_)st ^ | 


k=p,q 


exp I - ^SpgyWy^tt’p'tt’™ | 


with integers £1 and £2 independent of n. Note that according to (9.54), it suffices to consider 1^(0,0) 
for the proof of (9.41). We study Xc{i\,l- 2 ) for general £1 and £2 here, which will be used later. 

Now, we set 


Cp.g 


Cfe,i := -(a+ - a_)s£2/WroWe - -^yt\ . k=p,q, 
Ck ,2 := -(a+ -a_)s£2/WluWe''" - ■^2/Wr+, k=p,q. 


In addition, we introduce 

dp^q := dk := (iT k=p,q. 

Obviously, when (9.42) is satisfied, we have 

Cp,q = 0{dp^q), Ck,i = 0{dk), Ck,2 = 0{dk), k=p,q. 


(9.56) 


(9.57) 


(9.58) 
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With the aid of the notation defined in (9.55) and (9.56), we can write 

2i(7(^i,'^2) = exp I --^(yWrp + yWrg)| J dapdaq exp |i(n + £i)crW | exp |-i(n + £2)^^ | 

X exp|cfe,ie““^fc’+Cfc,2e’‘^^’|exp| - (9.59) 

k—p,q 

We have the following lemma. 


ii + £2 


(^3 : = 


Lemma 9.12. Under the truncation (9.42) and (9.43), we have 

for some positive constant C. 

Proof. At first, by Taylor expansion, we have 
exp |i(n + £i)crW| exp |-i(n + £2)crW | exp | 


= E 


(-1)' 


y(^) exp Ii(n + £1 + ni - 712)0-^1 exp I-i(n +£2 + ni - 712)0-^1 ■ ( 9 . 60 ) 


(711)1(772) 


'£^1 ^^2—0 

Now, for any mi, m 2 G Z, we denote 


ia{mi,m2) •.= / dfjWdfjW explimicrW} exp{-i7?720-W} TT exp Cfc, 

■7L2 


• [1] • [1] ' 


= ItT^ 1(773 + 7771 > 0) 


n3=0 


7731(773 + 777l)! 


k=p,q 
00 

1(774 + 7772 > 0 )- 


(Cp,i)"3+™7(Cp,2)"^ ^ ^„^(c,.l)”"(c,.2)”^+'"= 


114 = 0 


7741(774 + 7772)! 


Setting 


ni+n 2 


mi := 77 + £1 + 77l — 772, 7772 := 77 + £2 + 77l — 772, 

and using (9.60), we can rewrite (9.59) as 

711,712=0 

For simplicity, we employ the notation 

7773 := 7773 (£ 1 , 77i, 772, 773 ) = 7771 + 773 , 7774 := 7774(£2, 77i, 772, 774) = 7772 + 774. 

Consequently, by (9.63) and (9.61) we obtain 


(9.61) 


(9.62) 


Xcr (7771, 7772 ). (9.63) 


(9.64) 


0 °° 1 
|X.(£i,£ 2)| < 477^ exp { - + yl\)} E_^ {ni)l{n 2 )l 


^p,q 


M 


ni+n 2 


ic„ir=^ic„oi"3 “ 


E i(-3 > o)^:^7ir^^ 

^ nslms'. 


n 3=0 

< C max 

ni ,712 ,713 ,774,7773 ,7774 CN 

< C max 

rii ,712,713 ,’^'4 57713,7714 GN 

for some positive constant C, where in the last 


V i(„., > 


714=0 


( r \ ^ 1+^2 

-j^) {cp,ir^{cp,2r{cq,ir{cq,2r^ 

( r \ 711+712 

+,i)’"Mc,, 2)’"‘ 
step we used the fact 


(9.65) 


|cfc,i| < 1, |cfe,2| < 1, k = p,q 


(9.66) 
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which can be seen directly from the definition in (9.58), the truncations (9.42) and (9.43) and the as¬ 
sumption rj < Analogously, we also have 

Cri 


M 


< 1 . 


(9.67) 


According to the definitions (9.62) and (9.64), we have 

2(ni -I- 77,2) -\- + mA> 2n + + £ 2 - 

Hence, by using (9.66) and (9.67), we have the trivial bound 


max 

ni ,712 ,713 ,7?l4>0 




Therefore, we completed the proof by using (9.58). 


□ 


Now, we return to the proof of Lemma 9.11. Using (9.54) and Lemma 9.12 with £1 = £2 = 0 to (9.44), 
and integrating the bounded variables and a out, we can get 


|G(H,T)| < C 


r{Nr,)-B 


dyp 


[ 1 ] 


f{NT])S „ 

/ dj/W / 2tdt ■ (dp^q) 

Jo Jo 

X ((^) + dp” + d^”) exp { - 2Nr){a+ - a-)t^} 

V4-A2 
2 

k—p,q 


f'{Nr}) 4 


n+3 


X exp < — 


(1+0(1)) 


where the last two factors come from the facts 


exp. 

1-E4''( 

k=p,q 

exp. 

!-E^f'( 


= exp < - 


{ 


V4- A2 


E 4 "}. 


(9.68) 


k—p,q 


k—p,q 

X! 2^fc'((“+““-)^^+^“l^(l“2(ufe')^))| =1 + 0(1). 


In (9.68) we used the fact (u|3^)^ + (^^i^^)^ = 1- Now, we integrate yp^ and out. Consequently, by the 
definition in (9.57), we have 


|G(H,T)| < C 


^(A^?7) 4 


2 tdt 


(+ (;!?)'” + *'”) 1 ■ ■ “->‘d = e(?«;^) 


{NyY 


where in the last step we have used the assumption rj < M in (1.15), Assumption 1.13, the definition 

of 0 in (5.33) and the fact N = MW. Hence, we completed the proof of Lemma 9.11. □ 


Finally, we can prove Lemma 9.9, and further prove Lemma 9.7. 


Proof of Lemma 9.9 . This is a direct consequence of Lemmas 9.10 and 9.11. □ 

Proof of Lemma 9.7. This is a direct consequence of (9.37), Lemma 9.8 and Lemma 9.9. □ 

9.3. Summing up: Proof of Lemma 9.1. In this section, we slightly modify the discussions in Sections 
9.1 and 9.2 to prove Lemma 9.1. The combination of Lemmas 9.2 and 9.7 would directly imply Lemma 
9.1 if the Q(-) factor were not present in the definition of A(-). Now we should take Q(-) into account. 
This argument is similar to the corresponding discussion in Section 6.4. 




















75 


Proof of Lemma 9.1. At first, we observe that ki, K 2 and K 3 in (6.5) are obviously independent of n. 
Then, by the fact ki = it suffices to consider one monomial of the form 


(j/W) \(2/W) 



a.b=l,...,4 


where the degrees of pi(-), p 2 (-) and q(-) are all 0(1), and independent of n, in light of the fact K 3 = 0(1) 
in (6.5). Especially, the order of (?/p^^)“^and ( 2 / 5 ^^)“^ are not larger than 2, which can be easily seen from 
the definition of Q{-) in (3.18). 

Now, we reuse the notation Pq(A, B, V, T) and Fp(A, B, V, T) in (6.40), by redefining them as 
P,(l,S,E,r) := [ dOdS7’(.)-q({^^^},,,^i,...,^), 

F,{X,B.V,T):= j dXHdyI'lclwI'ld^(Pi)di/(Qi)P(.) 

xp.(<, s, (4")-'. ({e'-", J. 


It is easy to check P(-)q(-) also has an expansion of the form in (9.14). Hence, Pq(-) can be bounded in 
the same way as P(-)j thus we have 


|Pq(X,H,H,r)| < ^ |detA+pdet(^W)^. 

For Fp(-), the main modification is to use Lemma 9.12 with general £1 and £2 independent of n, owing 
to the function p 2 (-)- addition, by the truncations in (9.42) and (9.43), we can bound pi(-) by some 
constant C. Hence, it suffices to replace n by n + ^3 in the proof of Lemma 9.11. Finally, we can get 

r,ix.B.v,T) = o[j^^y 

with some finite integer £3 independent of n. Consequently, we completed the proof of Lemma 9.1. □ 


10. Integral over the Type II and III vicinities 


In this section, we prove Lemma 5.9. We only present the discussion forI(T()_, T(L, , T^, 

i.e. integral over the Type H vicinity . The discussion on I(T()_, T(!_, , Ti!i, is analogous. 

We start from (8.48). In this section, we will regard H-variables as fixed parameters, and consider the 
integrand as a function of all the other variables. 

Similarly, we shall provide an estimate for the integrand. At first, under the parameterization (8.2) 
with x: = +, we see that 


Xj,i-Xj,2 = -i;j-2) +o(^), bj,i+bj^2 = a+ -a_ +o(-^y 


y/M 


02 

y/M' 


Consequently, we have 

'V 12 \W r\- 

n(a;j,i - Xj^2f{bj,i + bj,2f = (»+ “ a_)2'^(l + 


1=1 


w 

n 


1=1 




( 10 . 1 ) 


Now, analogously to the case of Type I vicinity, what remains is to estimate A{X, B,V,T). Our aim 
is to prove the following lemma. 
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Lemma 10.1. Suppose that the assumptions in Theorem I.I 4 hold. In the Type II vicinity, we have 

\A{X,B,V,T)\ < e-^^’'|detA+|2det(S'(i))2 (10.2) 

for some positive constant c. 

With the aid of (10.1) and Lemma 10.1, we can prove Lemma 5.9. 

Proof of Lemma 5.9. Recall (8.48). At first, by the definition of Alj_ in (8.35), (5.23) and the fact Rea^ > 
0, we can see that 


Re(x'A!(_x) > ||x| 


(10.3) 


for all {Vj}Y =2 £ (^(2))^ Substituting (5.20), (10.1), (10.2), (10.3) and the estimates in Proposition 
8.6 into (8.48) yields 

, Ts, 

/ \2w r ^ ^ 

s ^ I **s'*’M *tA.,1. n ^ n 


i=2 


i=2 


/RW-i 


W 

n 


w 


w 


w 


w 


w 


7=2 


'RW-i 


/ Y[dxj,i / 


7=2 


7=1 


7=1 






X exp{(a+ - a_)^TiS'(^Vi}exp{(a+ - a_)^T25(^V2} exp{-i||ci ||2 - i||c 2 || 2 } 


X expt — 


1 1 w 

p{-2pill2- 2ll*2||l}n (%1 -%2 + o( 

7 = 1 


0 


where we absorbed several factors by exp{— cA^t;}. We also enlarged the domains to the full ones. Then, 
using the trivial facts 


w 


w 


n ^ L . n 2^'7di^7=1 


/lw-1 277 


/fW-l 


7=2 


and performing the Gaussian integral for the remaining variables, we can get 


/ / R \ \ 

|I(T^+,T^_,T^,T^,Ts,I^-i)| < C|det^W|- |detA+| • (^l + o(^^Jj 
Observe that 


(10.4) 


(10.5) 


|detA+| < |l + a;r < 2"^. 

Moreover, by Assumption 1.1 (ii), we see that |Sii| < (1 — co)/2 for some small positive constant cq. 
Consequently, since 5*^^^ is negative definite, we have 

'l-co\W^ 


i*tswi<ni...i<(i^) 


by Hadamard’s inequality. Substituting (10.5) and (10.6) into (10.4) yields 


(10.6) 


(10.7) 


|I(T']_,T'L,T^,T“,Ts,l'^-i)| = 0(e-“'^) 

for some positive constant S. Hence, we proved the first part of Lemma 5.9. The second part can be 
proved analogously. □ 


In the sequel, we prove Lemma 10.1. We also ignore the factor Q(-) from the discussion at first. 
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10.1. P(X,il,y,T) in the Type II vicinity. Our aim, is to prove the following lemma. 


Lemma 10.2. Suppose that the assumptions in Theorem 1.14 hold. In the Type II vicinity, we have 

w2+7C)2 

P{X,B,V,T) < ^ |detA+pdet(.gW)^. (10.8) 

Proof. We will follow the strategy in Section 9.1. We regard all 1^-variables as fixed parameters. Now, 
we define the function 

1= ij{X,B,T) := + \xj^2\ + + \bj,2\ + \ij\. 

Then, we recall the representation (9.5) and the definition of in (9.6). We still adopt the represen¬ 
tation (9.9). It is easy to see that in the Type II vicinity, we also have the bound (9.10) for pej.a./s- The 
main difference is the first factor of the r.h.s. of (9.5). We expand it around the saddle point as 

exp { - TrV;X-^V,n,T-^B-^T,E,} =: exp { - TrD4^n,D^%] exp { - ^A,}. 

We take the formula above as the definition of A^-, which is of the form 

4 

Aj = Pj,a,p ' 

where Pj,a,p is a function of X, B, V and T-variables, satisfying 

Pj,a,p = 0{i). 


Let 

H = (a+^A+) © 5 © (a;^A+) © S. 

Recalling the notation in (6.20), we can write 

w 

— SjkTrPtjAk — ’^^TrDff^TljD^'E.j = — 
j,k j=l 

Now, via replacing Aij- by Ai^-, kj by Ij, El by HI in the proof of Lemma 9.2, we can perform the proof 
of Lemma 10.2 in the same way. We leave the details to the reader. □ 


10.2. f{X,B,V,T) in the Type II vicinity. In this section, we will prove the following lemma. 


Lemma 10.3. Suppose that the assumptions in Theorem 1.14 hold. In the Type II vicinity, we have 

'exp{—(a+ — a-) Nr]} 


F{X,B,V,T) = o(^- 




(iV7?)"+i A 

Proof. Recall the decomposition (9.37). Note that Lemma 9.9 is still applicable. Hence, it suffices to 
estimate F(A, V). Now, note that in the Type II vicinity, it s obvious to see that 

+ I A,2 I 


TrXjJ = 0( 


Vm 


, Vj = I,...,W, 


which implies that 


Consequently, we have 




1=1 




exp 


^Xlrj'^TrXjj'^ = exp{0{&'/MT])} = 1 + o(I) 


1=1 
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by our assumption on ry. From (3.19) we can also see that all the other factors of /(Pi, V, X, are 
0(1). Hence, by the definition (9.36), we have 

F(X, V) = 0(exp{-(a+ - a-)Nr)}), 

which together with Lemma 9.9 yields the conclusion. □ 

10.3. Summing up: Proof of Lemma 10.1. Analogously, we shall slightly modify the proofs of Lemma 
10.2 and Lemma 10.3, in order to take Q(-) into account. The proof can then be performed in the same 
manner as Lemma 9.1. We omit the details here. 

11. Proof of Theorem 1.14 

The conclusion for Case 1 is a direct consequence of the discussions in Sections 3.5-10. The proofs of 
Case 2 and Case 3 can be performed analogously, with slight modihcations, which will be stated below. 

• (Case 1) 

In this case, by using Lemmas 5.1, 5.6, 5.8 and 5.9, we can get (1.19) immediately. 

• (Case 2) 

In this case, we shall slightly modify the discussions in Sections 3.5-10 for Case 1, according to the 
decomposition of supermatrices in (3.9). Now, at first, in (3.12) and (3.13), for A = 5, X, T, H or S, 
we replace A'^'^ and by and Aq respectively, and replace A^^ by Ap^ In addition, in the last 

three lines of (3.13), we shall also replace Sjq by Sjp, and replace Spq and Sqp by Spp, and in the first line, 
we replace ^i,g,i^i,p,i^ 2 .p,i((' 2 , 9 ,i by ^i,p, 2 (('i.p,i<^ 2 ,p,i?!> 2 ,p, 2 - Then, in (3.14) and (3.15), for A = X,Y, Vt, 
S, u), w, y, u, V or a, we replace A,^^ by Ap^. With these modifications, it is easy to check the proof 
in Sections 3.5-10 applies to Case 2 as well. The main point is we can still gain the factor 1/(A^77)"+^ 
from integral of y(-) defined in (3.20) (with and replaced by and wj,^^). Heuristically, we 
can go back to (4.6), and replace by a^p^ therein. It is then quite clear the same estimate holds. 
Consequently, Lemmas 5.1, 5.6, 5.8 and 5.9 still hold under the replacement of the variables described 
above. Hence, (1.19) holds in Case 2. 

• (Case 3) 

Analogously, in this case, we can also mimic the discussions for Case 1 with slight modifications. We also 
start from (3.12) and (3.13). For A = S, X, Y, H, S, uj and we replace Aq^^ by Aq, and replace a[^^ by 
0. In addition, in the first line of (3.13), we replace 4 >i,q,i4>i,p,i4>2,p,i(t>2,q,i by 4 >i,p,i4>i,p,i4>2,p,i4'2,p,i- Con¬ 
sequently, after using superbosonization formula, we will get the factor (yp^^|(wp^(wp^)*)i 2 |)^" instead 
of (yp^^y 5 ^^(wq^^(w 5 ^^)*)i 2 (wp^(wp^)*) 2 i)" in (3.16). Then, for the superdeterminant terms 

-pr detjXk - nk{Yk)-^Ek) -pr yk^ 

IL ’ kiL det2(xW) 

we shall only keep the factors with k = p and delete those with k = q. Moreover, we shall also replace 
a\^^ by 0 for A = X, Y, H, S, a;, w, y, u, V or a in (3.16). In addition, dA^^l shall be redefined 
as the differential of Ap ^-variables only, for A = X, y, w, w and One can check step by step that 
such a modification does not require any essential change of our discussions for Case 1. Especially, 
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note that our modification has nothing to do with the saddle point analysis on the Gaussian measure 
exp{—M{K{X, V) + L{B, T))}. Moreover, the term V{-) in (3.28) can be redefined by deleting the factor 
with k = q in the last term therein. Such a modification does not change our analysis In addition, 

the irrelevant term Q(-) can also be defined accordingly. Specifically, we shall delete the factor with k = q 
in the last term of (3.29) and replace by 0 for ^ = fl, S, u;, w, y. It is routine to check that 
Lemma 6.3 still holds under such a modification. Analogously, we can redefine the functions /(•) 
and g{-) in (3.18)-(3.20). Now, the main difference between Case 3 and Case 1 or 2 is that the factor 
(?/p'|(wW(wW)*)i 2|)2” does not produce oscillation in the integral of g{-) any more. Heuristically, the 
counterpart of (4.6) in Case 3 reads 

g(a+-a_)iv^ J dy[i]dwWdi^(Qi) • g{B, T, Qi, yW, w^) 

/*00 n [1] -I 

~ / 2tdt / daW • e-cArr,t^+cie—p t -^ 

7o k ^ Ng 

Hence, (1.19) holds for Case 3. 

Therefore, we completed the proof of Theorem 1.14. 


12. Further comments 

In this section, we make some comments on possible further improvements on our results. 


• (Comment on how to remove the prefactor in (1.19)) 


As mentioned in Remark 1.15, we have used N^° to replace /{NgY. However, the latter 

is also artificial. It can be improved to some n-dependent constant Cn via a more delicate analysis on 
A(-), i.e. the integral of Such an improvement stems from the cancellation in the Gaussian 

integral. At first, a finer analysis will show that the factor Q(-) can really be ignored, in the sense that 
it does not play any role in the estimate of the order of ]E|Gij(z)p". Hence, for simplicity, we just focus 
on the product P(-)F(-) instead of A(-). Then, we go back to Lemma 9.2 and Lemma 9.7. Recall the 
decomposition (9.37). A more careful analysis on F(-) leads us to the following expansion, up to the 
subleading order terms of the factors G(-) and F(-)! 


w 

FH = «»*■<■> ~ (jv5^^ S 

V ^ j-l 


( 12 . 1 ) 


where lj(-)’s and lj(’)'® are some linear combinations of the arguments. Analogously, we shall write down 
the leading order term of P(-) in terms of x, b, t and v explicitly, instead of bounding it crudely by 
using (9.24). Then it can be seen that the leading order term of P(-) is a linear combination of Xj^iXk, 2 , 
bj,ibk, 2 , Xj^abk,i 3 , Vj,aTk,p for j, fc = 1,..., IF and a, /3 = 1, 2, in which all the coefficients are of order 1/M. 
Observe that the Gaussian integral in (8.47) will kill the linear terms. Consequently, in the expansion 
(12.1), the first term that survives after the Gaussian integral is actually 


1 

(_/V7y)’T-+2 


Mg 

7^ 


w 

1=1 


Mg 


w 




j=i 


fo.2, 




( 12 . 2 ) 
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Replacing A(-) by the product of the leading order term of P(-) and (12.2) in the integral (8.47) and tak¬ 
ing the Gaussian integral over c, d, t and D-variables yield the true order l/{Nr])'^, without additional 
iV-dependent prefactors. 

• (Comment on the restriction \E\ < \/2 — k) 

This restriction is used in several places, we mention the two most important ones. 

The most critical issue is the term is{B, T) dehned in (5.6). A direct consequence of \E\ < ^/2—Hi is that 
Re{bj^i + bj^ 2 ){bk,i +bk, 2 ) > 0 for bi G T^ and b 2 G f^, thus Re£s{B,T) > 0. However, once \E\ > v^, 
Re£s{B,T) can be negative for bi G T'^ and b 2 G f’’^. Consequently, the measure exp{—M£s{B,T)} 
is not well defined, considering the domain of the t-variables is not compact. Actually, such a problem 
is unavoidable if we independently deform the contours of &yi’s (resp. bj^s) from R+ to any contour 
passing through the saddle point a+ (resp. a_), starting from the 0, since one can always choose and 
bk,i to be a+ while bj ^2 and bk ,2 to be very close to 0, such that Re(6jp -I- &y 2 )(&fc.i + bk, 2 ) < 0 in case 
\E\ > v^. A possible way to solve this problem is to change the variables (5j,i, &j, 2 ) to (r^p, rj_ 2 ), defined 
by Tju = (5j,i -I- 6 j, 2)/2 and rj ^2 = (&j,i — ^j,2)/2, and perform the saddle point analysis with respect to 
the latter. It is not difficult to calculate that the saddle point of L{B, T) restricted on these r-variables is 
1 = i/4 — and rj _2 = iA. To guarantee the positivity of £s{B,T), it suffices to choose a contour for 
passing through staying in the sector K := {w G C : | argu;| < 7r/4}. However, the explicit 

way to choose the contours for r^-^i and rj ^2 and the subsequent analysis in terms of r-variables instead of 
b-variables would be more involved and we leave it to future work. 

Another point where \E\ < n is used is the proof of Lemma 5.9 in Section 10. In (10.3), we used 
ReAlj. > /, which is a consequence of Rea^ = (2 —A^)/2 > 0 and (5.23), see the definition of A!(_ in (8.35). 
If |A| > v^, the bound in (10.3) should be weakened to Re(x'A!(.x) > 0(711) ||x||f for some A-dependent 
constant c{E) < 1, which is uniform on Vj G U{2) for all j = 2,..., W. This modification, however, is not 
critical. An easy way to remedy this situation is to impose an additional condition det < {c{E) — c)^ 
for arbitrarily small constant c, to guarantee the estimate (10.7). 
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